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ABSTRACT 


Interference  effects  between  a  highly  underexpanded,  sonic  or  supersonic 
jet  in  a  subsonic  or  supersonic  crossflow,  and  the  surface  from  which 
the  jet  exhausts  are  examined.  For  subsonic  free  stream  Mach  numbers, 
existing  data  is  examined  and  correlated.  Various  semi-empirical 
models  to  represent  the  interference  pressure  distribution  on  flat 
plates  are  then  developed.  For  supersonic  freestream  Mach  numbers, 
a  computer  program  for  calculating  jet  interference  effects  on  axisym- 
metric  bodies  at  angle  of  attack  is  described.  Interference  effects 
between  the  jet  plume  and  control  fins  on  a  cruciform  missile  are 
analyzed.  A  semi-empirical  model  of  the  jet  in  a  crossflow,  valid  at 
large  distances  from  the  nozzle  is  developed.  The  results  of  this 
model  are  then  used  to  compute  interference  forces  and  moments  on 
fins  located  aft  of  the  nozzle,  both  for  subsonic  and  for  supersonic 
freestream  Mach  numbers. 
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FOREWORD 


This  report  describes  results  of  an  analysis  conducted  by  the 
McDonnell  Douglas  Astronautics  Company--Western  Division 
(MDAC-WD),  under  United  States  Army  Contract  DAAHO 1 -68-C- 1 91 9. 
The  contract  was  initiated  under  DA  Project  No.  1M2623XXA206  and 
AMC  Management  Structure  Code  No.  522C,  11.  148.  The  technical 
effort  was  conducted  between  1  June  1968  and  31  August  1969.  The 
project  was  administered  under  the  direction  of  the  Aerodyna.mics 
Branch,  Advanced  Systems  Laboratory,  U.  S.  Army  Missile  Com¬ 
mand,  Redstone  Arsenal,  Alabama.  The  Army  technical  monitors 
for  the  study  were  Mr.  D.  J.  Spring  and  Mr.  T.  A.  Street. 

In  addition  co  the  authors,  Mr.  j.  G.  Davis  and  Dr.  R.  Rosen  of 
MDAC-WD  made  significant  contributions  to  the  boundary  layer  sepa¬ 
ration  analysis  and  subsonic  flow  modeling,  respectively. 
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Section  1 
INTRODUCTION 


The  problem  of  analyzing  the  effectiveness  of  reaction  jet  control 
systems  on  flight  vehicles  operating  in  the  atmosphere  has  received 
considerable  interest  in  the  past  several  years  (e.  g. ,  References  1 
to  11).  It  is  well  known  that  the  interaction  between  a  reaction  con¬ 
trol  jet  and  flow  around  the  jet  on  a  vehicle  surface  generates  a 
force  on  the  vehicle  which  is  often  larger  than  the  jet  thrust.  Con¬ 
sequently  the  problem  of  analyzing  reaction  jet  control  effectiveness 
usually  reduces  to  one  of  determining  the  magnitude  and  behavior  of 
the  force  due  to  this  interaction.  The  presence,  and  frequent  domi¬ 
nance,  of  the  force  due  to  interaction  leads  to  the  term  "jet  inter¬ 
action"  (JI)  control  which  is  usually  applied  to  endoatmospheric 
reaction  control  systems. 

Reports  of  both  experimental  and  analytical  studies  of  JI  are  common 
in  the  literature.  These  studies  are  generally  classified  according  to 
mainstream  flow  conditions  and  the  jet  configuration.  The  case  of 
widest  interest  (e.  g.  ,  References  5,  12,  13,  14,  15)  has  been  that  of 
two-dimensional  interaction  between  the  jet  from  an  infinite  sonic 
slot  and  a  uniform  supersonic  stream  over  a  flat  plate.  The  inter¬ 
action  is  considered  to  be  three-dimensional  whenever  a  velocity 
component  exists  normal  to  the  plane  of  intersection  of  the  mainstream 
velocity  and  the  jet  centerline  in  the  interaction  flow  region.  The 
three-dimensional  interaction  between  a  jet  from  an  orifice  or  nozzle 
in  a  flat  plate  and  a  uniform  supersonic  stream  has  received  attention 
from  various  investigators  (e.  g.  ,  References  2,  16,  17,  18).  More 
complex  JI  problems  involving  variously  configxired  sonic  and  super¬ 
sonic  jets  exhausting  transverse  to  axisymmetric  or  three-dimensional 
supersonic  flows  are  more  infrequently  discussed  in  the  unclassified 
literature  (e.  g.  ,  References  1,  2,  3,  6,  8).  Finally,  studies  of 
under  expanded  sonic  or  supersonic  jets  interacting  with  subsonic 
mainstreams  are  comparatively  rare  for  any  flow  geometry  (e.  g.  , 
References  8,  9,  10,  19). 

Some  success  in  scaling  experimental  data  for  JI  control  forces  has 
been  demonstrated  in  the  literature,  notably  for  two-dimensional 
flows.  However,  JI  control  force  prediction  techniques  for  complex 
flow  geometries  typical  of  those  encountered  in  application  are  rare 
and  generally  very  limited  in  range  of  applicability.  The  blast  wave 
analogy  methods  proposed  by  various  investigators,  as  discussed  in 
Reference  19,  have  been  most  commonly  applied. 
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The  objective  of  the  study  reported  here  has  been  to  develop 
approximate  JI  control  force  prediction  techniques  applicable  to  JI  in 
three-dimensional  ..subsonic  and  supersonic  flows.  In  the  initial  phase 
of  the  study  (Reference  19),  emphasis  was  placed  on  devising  a 

"  technique  applicable  to  JI  controls  on  axisymmetric  missiles  in 
supersonic  flight,  and  a  lesser  effort  was  devoted  to  the  subsonic 
flight  problem.  In  the  second  phase  of  the  study,  emphasis  has 
been  placed  on  the  development  of  a  technique  for  analyzing  the  JI 
problem  when  the  mainstream  is  subsonic.  A  secondary  effort  in  the 
second  phase  has  been  devoted  to  expanding  the  range  of  applicability 
of  the  supersonic  mainstream  analysis  technique  developed  in  the 
initial  phase. 

4  The  general  complexity  of  three-dimensional  JI  flows  is  well  known, 
particularly  with  regard  to  their  boundary  layer  separation  and  jet 
plume  aspects.  Consequently,  approaches  to  interaction  force 
prediction  methods  are  usually  through  analogies  to  the  interaction 
flow  field  rather  than  descriptions  of  it.  This  has  been  the  case  in 
the  present  study.  JBased  on  experimental  data  available  in  the 
literature,  or  provided  by  the  U.  S.  Army  Missile  Command 
(AMICOM)  from  recent  experiments,  the  significant  governing  aspects 
of  the  flow  field  have  been  identified.  Then  inviscid  flow  analogies 
have  been  developed  to  represent  the  various  aspects  of  the  flow  field. 
This  approach  has  been  taken  with  both  supersonic  and  subsonic 
mainstream  JI  with  various  degrees  of  success  as  discussed  in 
Reference  19  and  this  report. 

Since  the  second  phase  of  the  study  has  been  strongly  oriented  toward 
JI  in  a  subsonic  mainstream,  the  bulk  of  this  report  is  devoted  to 
this  subject.  The  general  nature  of  the  JI  flow  field  in  a  subsonic 
mainstream  is  discussed  along  with  empirical  models  of  its  behavior 
in  Section  2.  Various  incompressible,  potential  flow  models  of  the 
subsonic  mainstream  JI  flow  are  discussed  in  Section  3.  The  expanded 
equivalent  body  analogy  JI  effectiveness  prediction  program  for  three-  i 
dimensional  jets  on  axisymmetric  missiles  in  supersonic  flight  is 
described  in  Section  4.  Methods  of  calculating  the  effects  of  jet-to-fin 
interference  when  fins  are  located  aft  of  JI  controls  in  subsonic  or 
supersonic  flight  are  discussed  in  Section  5. 
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Section  2 


SUBSONIC  MAINSTREAM  JET  INTERACTION 
DESCRIPTION  AND  SCALING 


The  development  of  mathematical  models  of  JI  in  a  subsonic  mainstream 
requires  a  basic  understanding  of  the  fluid  mechanics  involved  which 
can  only  be  derived  from  detailed  experiments.  In  this  section  the 
physical  aspects  of  the  interaction  flow  will  be  described,  based  on 
available  experimental  data..  Empirical  scaling  of  interference  pres¬ 
sures  in  the  interaction  region  will  also  be  described. 

Reaction  control  jet  systems  typically  employ  very  high  chamber 
pressures,  so  that  downstream  of  the  nozzle,  the  jet  exhibits  the  inter¬ 
nal  shock  system  characteristic  of  highly  underexpanded  plumes  (see 
Figure  1).  It  will  be  shown  later  that  a  characteristic  dimension 
of  this  shock  system  is  an  important  scaling  length  for  the  interfer¬ 
ence  pressure  distribution. 

No  experiments  were  conducted  by  MDAC-WD  during  the  study.  However, 
previously  unpublished  experimental  data  taken  by  AMICOM  was  made 
available  for  use  in  the  study.  To  the  knowledge  of  the  authors,  'the 
experimental  data  concerning  an  underexpanded  jet  in  a  subsonic  main¬ 
stream  which  is  being  provided  by  AMICOM  is  the  only  reasonably 
detailed  data  in  existence  (References  10  and  20  contain  limited  data). 
The  more  detailed  experiments  by  AMICOM  are  not  yet  complete; 
consequently,  experimental  data  regarding  a  subsonic  jet  in  a  crosswind 
have  been  relied  on  heavily  during  the  present  study.  Existing  data 
may  be  classified  into  two  categories:  (1)  Flow  field  surveys  in  the 
vicinity  of  the  jet  plume,  along  the  jet  trajectory  and  (2)  pressure  dis¬ 
tribution  on  the  surface  from  which  the  jet  exhausts  or  forces  and 
moments  on  the  body  from  which  the  jet  exhausts. 

2.  1  DESCRIPTION  OF  THE  JET 

The  flow  field  of  interest  is  illustrated  in  the  schlieren  photograph  in 
Figure  1.  Reaction  control  jets  are  typically  highly  unde rexpanded, 
causing  a  shock  engulfed  plume  at  the  jet  exit  similar  to  that  evident 
in  the  figure.  The  existence  of  this  plume  differentiates  the  subsonic 
mainstream  JI  flow  field  from  that  often  studied  in  reference  to  the 
vertical  take-off  and  landing  (VTOL)  aircraft  transitional  flight  prob¬ 
lem.  VTOL  related  studies  typically  deal  with  low  subsonic  ;et 
velocities . 

To  the  authors'  knowledge,  no  surveys  of  the  flow  in  a  highly  under¬ 
expanded  rocket  plume  exhausting  normal  to  a  subsonic  free  stream 
have  been  conducted.  Examination  of  flow  visualization  photographs 
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(b) 

FjqUre  i.  Jet  in  Subsonic  Cross  Flow  with  Barrel  Shock  System 


such  as  Figure  1  indicates  that t  the  plume  shape  is  not  influenced  too 
strongly  by  the  crossflow,  indeed,  since  the  flow  within  the  barrel 
shock  system  is  highly  supersonic,  most  of  the  flow  field  in  this 
portion  of  the  jet  coi’e  would  be  unaffected  by  any  changes  in  the 
plume  boundary.  Near  the  nozzle  exit,  a  highly  underexpanded 
plume  in  a  .subsonic  crossflow  might  therefore  be  expected  to 
behave  approximately  as  a  jet  plume  exhausting  into  still  air. 

The  jet  plume,  however,  has  a  blockage  and  entrainment  effect  on 
the  mainstream,  in  the  neighborhood  of  the  nozzle  exit.  Some 
aspects  of  this  behavior  can  be  surmised.  A  shear  layer  forms 
around  the  jet  plume  and  mixing  between  the  jet  and  the  mainstream 
begins  in  this  region  where  the  jet  axial  velocity  is  relatively  low. 

The  mainstream  flow  past  the  plume  and  shear  layer  is  believed  to 
behave  in  a  manner  similar  to  flow  over  a  cylinder  at  low  Reynolds 
number  with  separation  and  the  formation  of  vortices  on  the  leeward 
side.  The  behavior  of  the  mixing  region  is  more  evident  in  Figure  2 
where  the  underexpanded  jet  is  supersonic  at  the  exit  and  the  plume 
exhibits  the  familiar  diamond  shock  pattern.  Since  the  supersonic 
jet  shown  has  an  expansion  ratio  of  4.  0,  it  is  less  underexpanded 
than  the  sonic  jet  in  Figure  1,  and  the  shock  bounding  the  plume  is 
not  as  strong. 

At  some  distance  from  the  nozzle  exit,  the  jet  eventually  becomes 
subsonic,  probably  before  significant  bending  occurs.  From  this 
point,  more  definitive  descriptions  of  the  jet  are  possible  based  on 
the  more  detailed  experimental  data  for  a  subsonic  jet  in  a  crosswind. 

Several  excellent  observations  of  the  qualitative  behavior  of  subsonic 
jets  in  a  subsonic  crossflow  exist  in  the  literature,  such  as  those  in 
References  21  and  22.  The  remarks  in  this  paragraph  are  based 
mainly  on  these  references.  Beginning  at  the  nozzle  and  progressing 
along  the  jet  axis,  the  first  region  encountered  consists  of  a  potential 
core  surrounded  by  a  turbulent  mixing  layer.'  The  low  momentum 
flow  in  the  mixing  layer  is  deflected  downstream  by  the  crossflow, 
causing  a  deformation  of  the  jet  cross-section  into  a  kidney  shape. 

The  flow  separates  near  the  edges  cf  the  jet  and  two  counter¬ 
rotating  vortices  are  formed  on  the  leeward  side,  as  in  low  Reynolds 
number  flow  about  a  circular  cylinder.  The  potential  core  is  con¬ 
sumed  in  a  shorter  distance  than  if  the  jet  were  exhausting  into  still 
air,  but  its  centerline  remains  undeflected,  except  in  the  case  of 
relatively  low  jet  velocities  compared  to  free  stream  velocities. 

After  the  jet  has  become  fully  turbulent,  vigorous  mixing  with  the 
free  stream  occurs.  The  vortices  on  the  leeward  side  apparently 
enhance  the  entrainment  of  external  air.  This  mixing  causes  the 
jet  axial  velocity  to  decay  much  more  rapidly  than  for  a  jet  in  still 
air-,  as  illustrated  by  Figure  1 1  of  Reference  23.  The  entrainment 
of  free  stream  momentum  causes  the  jet  to  bend  quickly  to  a  direc¬ 
tion  which  is  almost  aligned  with  the  free  stream.  It  is  evident 
from  experimental  data  such  as  that  in  Reference  23  that  mixing  or 
momentum  entrainment  and  not  a  "cross-flow  pressure  drag"  causes 
most  of  the  jet  bending.  Some  analogies  between  the  cross  flow  pres¬ 
sure  drag  on  a  solid  obstacle  and  the  jet  bending  have  been  drawn, 
but  experiments  seem  to  indicate  that  the  pressure  in  the  jet  soon 
adjusts  1o  the  free  stream  value  and  any  pressure  drag  influence  is 
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cGiiliiicu  to  i'ne  immediate  neighborhood  of  the  orifice.  For  example, 
Jordiaon  states  in  Reference  21  that  surveys  show  that  the  pressure 
in  the  jet  is  everywhere  equal  to  the  free  stream  pressure  a  few 
diameters  from  the  orifice.  Some  more  evidence  of  the  rapid  decay 
to  free  stream  pressure  is  shown  in  Reference  24,.  In  the  zone  of 
maximum  bending,  the  leeward  vortices  have  been  absorbed  into  the 
jet  and  continue  to  grow  in  strength.  The  final  state  in  the  jet  appears 
<o  be  a  long  region  where  the  jet  direction  differs  little  from  the  free 
stream,  and  where  the  axial  jet  velocity  is  practically  the  same  as  the 
free  stream  velocity.  This  last  region  is  called  the  "vortex  zone"  in 
Referent:  22  because  the  counter-rotating  vortices  still  persist, 
although  with  diminishing  strength. 

....  -  ( _t  *1 

2.  2  UNDER  EXPANDED  JET  INTERACTION  PRESSURE 
i  DISTRIBUTIONS 

Several  experiments  have  been  conducted,  by  the  Advanced  Systems 
Laboratory,  ..  Research  and  Engineering. ^Directorate,  AMICOM, 
whtlre  pressure  distributions  due  to.  interaction  between  an  under- 
expkttded  .jet  and  ,a  subsonic  stream,  were  measured.  Except  for  the 
data  presented  by  Spring  and  Street  in  Reference  8,  results  of  these 
(kxperiment8.are  not  yet  generally. available.  In  the  present  study, 
data  from  'wind  tunnel  tests 'conducted  by  AMICOM  at  Cornell 
Aeronautical  Laboratory  (CAL.)  and  Arnold  Engineering  Development 
Center  (AEDC)  were  tnad«e  available  to  the  authors  by  Spring  and 
Street.  These  data,  some  of  yvhich  will  be  presented  in  this  section, 
are  contained  in  References  25  to  28,  which  are  not  generally  avail¬ 
able  except  through  AMICOM.  Schematic  drawings  of  the  test  models 
are  shown  in  Figures  3  anfl  4. 

■/  '• 

References  25  and  26  have  limited  pres  sure  data  for  a  highly  under¬ 
expanded  jet  exhausting  from  a  flat  plate.  The  tests  of  Reference  25 
were  conducted  with  several  different  nozzle  configurations  of  slightly 
different  exit  Mach  numbers.  Different  jet  gases  were  also  used  in 
order  to  investigate  the  effect  of  changing  the  jet  specific  heat  ratio, 

Yj.  The  free  stream, Mftch  number  was  varied  from  0,  6  to  1.  2,  and 
the  jet  chamber, -to-free-stream  pressure  was  varied.  The  experi¬ 
ment  reported  in  Reference  26  was  conducted  with  essentially  the 
same  model,  but  included  tests  of  hot  gas  effects. 

References  27  and  28  contain  extensive  interference  pressure  data 
for  an  underexpanded  jet  exhausting  just  forward  of  the  nose  juncture 
from  an  ogive- cylinder  missile  configuration.  The  ogive- cylinder 
model  of  Reference  27  had  interchangeable  circular  and  slot  nozzles. 
The  two  circular  nozzles  had  different  exit  diameters  but  the  same 
exit  Mach  number.  The  model  of  Reference  28  was  tested  with  one 
circular  nozzle  and  one  slot  nozzle,  both  of  which  were  sonic. 

During  these  tests,  the  free  stream  Mach  number  varied  from  0.  20 
to  1.  25  and  the  jet  chamber-to-free-stream  pressure  ratio  ranged 
from  0  to  120.  The  pressure  distribution  in  the  neighborhood  of 
the  nozzle  was  measured  with  a  large  number  of  taps  distributed  on 
the  surface.  Total  forces  and  moments  on  the  model  were  also 
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Figure  4.  Army  Missila  Command  Modal 


measured  through  a  sting  balance  system.  During  some  of  the 
experiments  rectangular,  cruciform  stabilizing  fins  were  placed 
on  the  model.  Two  of  the  fins  were  instrumented  with  pressure  taps, 
and  the  other  two  had  their  own  internal  balances  to  measure  forces 
and  moments. 

Generally,  the  pressure  distributions  measured  in  References  25  to 
28  exhibit  the  same  behavior  as  those  observed  for, a  supersonic  main¬ 
stream.  There  is  a  limited  region  of  positive  pressure  coefficient  on 
the  windward  side  of  the  jet  and  a  larger  negative  Cp  region  on  the 
leeward  side.  It  has  been  found  that  it  is  possible  to  correlate  some 
of  the  experimental  results  by  properly  choosing  the  scaling  length. 

The  ogive -cylinder  data  of  References  27  and  28  will  be  discussed  first, 
since  the  pressure  distributions  were  measured  in  greater  detail. 

Using  the  coordinate  system  shown  in  Figure  5,  pressure  coefficients 
along  the  plane  of  symmetry  S  =  0  are  plotted  in  Figure  6  for  one  value 
of  the  free  stream  Mach  number  and  various  pressure  ratios.  The 
pressure  coefficient  has  been  defined  in  the  conventional  manner: 


C 
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Figure  5.  Reference  Coordinate  Syrtem  for  Interference  Pressure  on  an  Ogive-Cylinder 
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Although  the  nozzle  was  not  quite  located  on  the  cylindrical  portion 
of  the  model,  the  effects  of  curvature  have  been  neglected  and  x  has 
been  assumed  to  be  equal  to  the  distance  along  the  model  axis  in 
inches.  Also,  the  jet-off  pressure  distribution  has  not  been  sub¬ 
tracted  out,  so  that  the  Cp  shown  in  Figure  6  is  not  strictly  an 
interference  pressure  coefficient.  The  data  shows,  however,  that 
the  jet-off  Cp  is  very  small  and  has  a  negligible  effect  on  the  curves 
of  Figure  6.  These  curves  exhibit  the  characteristic  positive 
pressure  coefficients  on  the  windward  (x<0)  side,  with  large  negative 
pressure  coefficients  on  the  leeward  (x>0)  side. 

The  jets  used  in  References  27  and  28  are  highly  underexpanded  and 
possess  the  internal  shock  structure  characteristics  shown  in  Fig¬ 
ure  2.  It  has  been  found  that  a  significant  dimension  of  this  shock 
structure  provides  a  reasonably  valid  scaling  length  for  pressure 
distributions  in  the  neighborhood  of  the  jet  exit.  Based  on  flow 
visualization  data  suchas  shown  inFigures  1  and  2itcanbe  assumed 
that,  insofar  as  internal  shocks  are  concerned,  jets  penetrating  a 
subsonic  cross  flow  behave  as  if  they  were  exhausting  into  still  air. 
Let  h  denote  the  distance  from  the  nozzle  exit  to  the  Mach  disk  or 
to  the  first  intersection  ol  the  "diamond"  shock  pattern  when  this 
configuration  exists.  The  data  of  Love  el  al  in  Reference  29  indi¬ 
cate  that,  for  both  shock  configurations,  the  ratio  of  h  to  the  jet  exit 
diameter  can  be  correlated  as  a  function  of  the  jet  exit  to  free  stream 
pressure  ratio  at  a  fixed  value  of  the  jet  exit  Mach  number.  For 
sonic  nozzles,  Crist,  Sherman,  ahd  Glass  have  correlated  experimental 
values  of  h  for  a  wide  range  of  conditions  by  plotting  (h/dj)  vs  the  jet 
stagnation  to  ambient  static  pressure  ratio  P,  as  shown  in  Refer¬ 
ence  30.  Based  on  experimental  data,  they  obtain  the  empirical 
equation: 

J  =  (0.645)  V"P"  (1) 

j 

Although  this  equation  does  not  strictly  apply  to  a  diamond  shock 
pattern  and  a  supersonic  nozzle,  it  will  be  used  to  calculate  h  for 
scaling  purposes  with  the  nozzle  throat  diameter  dj.  substituted  for 
dj.  A  length,  h,  thus  calculated,  permits  correlation  of  interference 
pressure  data  for  different  pressure  ratios  and  nozzle  diameters  for 
an  underexpanded  jet  in  a  subsonic  crossflow. 

Figure  7  shows  the  same  data  as  Figure  6,  but  with  x  scaled  by  the 
shock  intersection  height  h.  Evidently,  the  data  for  all  pressure 
ratios  fall  on  a  s;ngle  curve.  Data  for  larger  diameter  nozzle 
(dt  =  0.  33  in.  )  are  also  included  in  the  figure,  and  the  points  corre¬ 
late  well  with  dt  =  0.  22  in.  data.  Figure  8  shows  the  same  data  cor¬ 
relation  for  cases  when  the  free  stream  Mach  number  is  0.  20.  Note 
that  in  this  case  Cp  is  negative  upstream  as  well  as  downstream  of 
the  nozzle.  The  data  for  M^  =  0.  20  has  also  been  correlated  along 
the  line  x  =  0,  as  shown  in  Figure  9.  The  abscissa  represents  the 
arc  length: 

S  =  Rtp 

normalized  by  h. 
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Figure  6.  Pressure  Distribution  Along  Plane  of  Symmetry  of  an  Ogive 


IRE  COEFFICIENT.  Cp 


n 


4 


Flfurt  9.  Comlattd  Circumfaianttol  Fnaura  Distribution  at  tho  Jtt  Station,  ■  Oi 


t 


15 


v  10.  Rtftnnca  Coordinate  Syatam  for  DaacHbing  Intarfannct  Praaura  Distribution  on  a  Flat  Plata 


When  the  flat  plate  interference  pressure  measurements  of  Refer¬ 
ence  25  are  evaluated,  it  is  found  that  the  scale,  h,  correlates  this 
data  also.  The  coordinate  system  used  is  shown  schematically  in 
Figure  10.  Figure  11  is  a  plot  of  correlated  data  along  the  centerline 
of  the  plate,  and  Figure  12  shows  a  limited  amount  of  data  at  inter¬ 
mediate  values  of  0.  It  must  be  pointed  out  that  the  results  shown  in 
Figures  11  and  12  are  not  as  conclusive  as  those  on  the  previous 
graphs  because  the  pressure  ratios  are  not  very  different.  Further¬ 
more,  the  jet  exit  Mach  numbers  in  the  cases  shown  arc  not  exactly 
the  same,  and  it  appears  that  the  correlation  only  holds  for  a  fixed 
exit  Mach  number,  as  d'seussed  below. 

In  sum,  it  has  been  found  that  a  scale  which  is  characteristic  of  the 
internal  shock  structure  in  a  highly  underexpanded  jet  and  varies 
directly  as  the  square  root  of  the  pressure  ratio  correlates  the  Inter¬ 
ference  pressure  distribution  for  fixed  and  Mj.  A  striking  feature 
of  the  correlated  data  is  the  limited  extent  of  the  “disturbance  in  the 
flat  plate  as  well  as  the  ogive-cylinder  cases.  It  appears  that  the 
induced  pressures  reduce  to  zero  within  four  to  five  lengths,  h,  from 
the  nozzle. 

Pressure  coefficients  from  Reference  28  for  tests  with  a  sonic  nozzle 
were  compared  to  the  Mj  =  2.  94  data  at  equal  values  of  M^.  The 
curves  did  not  agree.  Consideration  was  given  to  improving  the  cal¬ 
culation  of  Mach  disk  height  to  include  the  dependence  on  nozzle  Mach 
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Figure  11.  ComtoMd  Praam  0»trib«tini  Akreg  PUm  of  SyrenMiy  of  •  Ftat  Pktt 
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Figure  12.  Correlation  of  Flat-Plate  D.tta  off  the  Centerline 
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number,  M:,  in  order  to  extend  the  correlation  to  other  nozzle  Mach 
numbers.  A  method  for  computing  h  in  a  highly  underexpanded 
supersonic  plume  exhausting  into  still  air  was  developed,  as  discussed 
in  Appendix  A.  However,  this  stitl  failed  to  correlate  the  effects  of 
nozzle  exit  Mach  number. 

The  interaction  forces  and  moments  produced  by  the  interference 
pressures  due  to  the  underexpanded  jet  in  a  subsonic  stream  were 
measured  as  noted  above  in  the  tests  described  in  References  27  and 
28.  For  the  ogive- cylinder  configuration  shown  in  Figure  4,  the  jet 
thrust  amplification  faccor,  defined 

T.  _  Thrust  f  Interaction  Force 

A  "  Thrust 

is  shown  in  Figure  13  as  a  function  of  the  momentum  ratio  parameter 
described  in  Reference  19.  It  is  evident  that  the  amplification  factor 
scales  reasonably  well  with  the  momentum  ratio  parameter,  even 
when  jet  throat  area  and  free  stream  Mach  number  are  varied.  In 
fact,  for  these  low  subsonic  Mach  numbers  the  correlation  is  better 
than  the  same  correlation  for  a  similar  configuration  at  the  tran¬ 
sonic  and  supersonic  Mach  numbers  shown  in  Reference  31.  The 
effect  of  exit  Mach  number  on  amplification  factor  is  not  correlated 
by  the  momentum  ratio,  just  as  its  effect  on  pressure  distribution 
was  not  correlated  by  the  plume  dimension,  as  described  above. 

2.3  EMPIRICAL  MODELING  OF  THE  INTERFERENCE 

PRESSURE  DISTRIBUTION 

Since  the  objective  of  the  study  reported  here  has  been  to  expedi¬ 
tiously  develop  reliable  engineering  analysis  techniques,  the  first 
analysis  models  developed  were  empirical.  The  empirical  models 
are  not  only  of  direct  use,  they  are  also  of  considerable  value  in  the 
development  of  more  analytical  analysis  models,  as  will  be  described 
in  Section  3.  In  this  section,  experimental  interference  pressure 
distributions  on  a  flat  plate  are  represented  by  means  of  a  Fourier 
series  in  the  azimuthal  angle,  0,  defined  in  Figure  10. 

The  principal  difficulty  encountered  in  the  development  of  empirical 
models  of  the  interference  pressure  distribution  was  the  lack  of 
sufficiently  detailed  experimental  data.  The  data  in  References  27 
and  28  are  representative  of  the  level  of  detail  required;  however, 
the  three-dimensional  nature  of  the  undisturbed  flow  leads  to  a 
special  case  in  the  interference  pressure  distribution.  The  data  in 
References  25  and  26  for  the  jet  exhausting  into  a  uniform  stream 
are  ideal  for  differentiating  the  effects  of  interaction;  however,  the 
pressure  distributions  were  not  measured  in  sufficient  detail  to 
provide  data  for  empirical  modeling.  Experiments  are  presently 
being  conducted  by  AMICOM  to  provide  more  detailed  data  for  the 
behavior  of  an  underexpanded  jet  exhausting  from  a  fl'"  plate  lr 
uniform  subsonic  crossflow.  In  the  interim,  data  for  interierence 
pressures  due  to  a  subsonic  jet  as  reported  in  Reference  32  by  Vogler 
have  been  used.  Comparison  of  the  limited  data  Reference  25 

with  the  data  from  Reference  32  indicates  that  the  general  shape  and 
levels  of  the  interference  pressure  distribution  are  approximately  the 
same  in  both  cases. 
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For  subsonic  jets,  References  32  to  34  indicate  that  the  significant 
correlating  parameter  is  the  ratio  of  free  stream  to  jet  velocity 
(Uoo/Uj  ).  Geometrical  similarity  also  exists,  so  that  data  for  the 
same  velocity  ratio,  normalized  by  the  nozzle  diameter,  falls 
approximately  on  a  single  curve.  Plots  of  Vogler’s  data  at  a  fixed 
value  of  r  (where  r  has  been  normalized  by  the  jet  diameter) 
reveal  that  for  a  fixed  value  of  (Uoo/Uj)  the  data  varies  regularly 
with  0.  It  appears  therefore  that  a  truncated  Fourier  series  of  the 
form 

m  /  IT  \ 

Cp(r.  9;  U./U.)  =  cn  ^r,  ^  cos  =8  i2> 

should  represent  the  data  quite  well  with  relatively  few  terms.  The 
coordinates  used  in  Equation  (1)  are  depicted  in  Figure  10.  Due  to 
the  symmetry  of  the  flow  about  the  lines  0=0  and  0  =  rr ,  the  series 
will  not  contain  any  sine  terms.  The  coefficients,  cn,  may  be 
evaluated  at  a  fixed  value  of  r  by  .numerically  integrating  the  data 
and  using  the  orthogonality  of  the  cosine  function.  That  is,  the 
expressions 

and  ^ 

c  ( r ;  U  /U.)=  —  [  C  cos  n0  d0 
nl  “  J/  *4  P 

fcr  n>0 
can  be  used. 

Some  difficulty  was  encountered  in  using  Vogler's  data,  principally 
because  the  measured  pressure  coefficient  does  not  decay  to  zero 
as  it  should.  As  mentioned  in  Reference  32,  this  fact  is  probably 
caused  by  misalignment  of  the  plate  with  respect  to  the  free  stream 
or  possibly  by  warping  of  the  plate  under  the  loads  induced  by  the  jet. 
The  data  was  therefore  adjusted  along  each  ray  9  =  constant  so  that 
the  pressure  coefficient  would  be  zero  at  r  =  10.  The  results  of  plot¬ 
ting  the  data  for  U^/Uj  =  0.  3  as  a  function  of  9  for  various  values 
of  the  normalized  radius  r  are  shown  in  Figures  14a-  14d.  Using 
Equations  (3)  and  (4),  the  data  has  been  integrated  numerically  to 
obtain  the  coefficients  cn.  Two-  and  three-term  series 
representations 

Cp  =  CQ(r)  +  ^(r)cos  9 

and 

C  =  c  (r)  +  c,(r)cos  9  +  c-(r)cos  20 
p  o  1  2 

are  also  plotted  in  Figures  14a- i4d.  Evidently,  three  terms  in  the 
series  are  sufficient  to  represent  the  pre.ssure  distribution  quite  » 
well . 
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(3) 

(4) 


AZIMUTHAL  ANGLE,  0  (DEGREES) 


Figurg  14b.  Comparison  of  Voglsr's  Data  and  Fourw  Sarin,  r  =  2.0 
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Figvra  14c  Comparison  of  Voglar's  Data  and  Fouriar  Sariaj,  r  *  3.0 


AZIMUTHAL  ANGLE.  0  (DEGREES) 


Figure  14d.  Comparkcn  of  Voglar's  Data  and  Fouriar  Sarias,  r  =  4.0 


There  is  one  other  distinct  advantage  to  this  truncated  Fourier 
series  technique.  The  force  or  moment  on  the  plate  due  to  the 
interaction  can  be  calculated  from  the  equation 


F  or  M 


C  cos  n0  r^n+^  dr  dG 
P 


(5) 


where  the  force  is  obtained  if  n  equals  0  and  the  moment  results 
if  n  equals  1.  If  the  pressure  coefficient  is  written  in  the  form  of 
Equation  (2),  the  integral  becomes 


2.  C» 


n  s  m 


cos  m0  r^+n^dr  d9 


rj  0 


n=  0 


or 


m  , 

^  /*/  n  cn  C0S  n®  COS 


r(1+n)  dr  dG 


*—<  -v  -'n 
n  =  0  rj  0 


Due  to  the  orthogonality  of  the  cosine,  this  now  becomes 

; 

and 

n2  it  2  2 

c .  cos  0r  dr  dG 

J  1 


F  =  q  f«c  r  dr  dG 


Thus,  the  total  contribution  to  the  forces  comes  from  the  first  term 
of  the  Fourier  series  (truncated  or  complete)  and  the  total  contribu¬ 
tion  to  the  moment  comes  from  the  second.  Matching  the  first 
couple  of  terms  of  a  Fourier  series  takes  on  a  new  significance  in 
the  light  of  this  result.  If  analytical  models  could  be  found  that 
would  have  very  close  agreement  in  these  terms  for  all  values  of  r, 
then  the  two  most  important  quantities  could  be  predicted  quite  well. 

Figures  15  and  16  show  the  first  two  Fourier  coefficients  obtained 
from  all  of  Vogler's  data  as  functions  of  r  for  various  values  of 
U^/Uj.  The  variation  in  the  third  coefficient,  c,,  is  much  more 
irregular,  as  shown  in  Figure  17. 


From  the  above  results,  it  is  evident  that  Fourier  series  techniques 
will  yield  a  simple  and  efficient  empirical  description  of  interference 
pressure  distributions  on  a  flat  plate.  In  addition,  this  means  of 
data  analysis  has  been  helpful  in  the  development  of  semi- empirical 
flowfield  models,  as  described  in  the  next  section. 
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Figure  1 


Section  3 

FLOWFIELD  MODELS  FOR  THE  INTERACTION  REGION 


The  subsonic  mainstream  JI  flowfield  is  characterized  by  complex 
phenomena,  including  mixing  of  turbulent  and  laminar,  compressible 
and  comparatively  incompressible  flows,  as  well  as  three-dimensional 
flow- separation  phenomena.  Consequently,  detailed  mathematical 
modeling  of  the  flow  represents  a  task  of  considerably  higher  magni¬ 
tude  than  that  intended  in  the  study  reported  here.  In  this  study,  the 
approach  taken  to  mathematical  modeling  was  indirect.  Analogies  to 
the  actual  flowfield  were  postulated  that  could  be  expected  to  yield 
pressure  distributions  in  the  region  of  the  jet  exit  which  would 
behave  as  those  derived  from  experiment.  The  model  flows  were 
assumed  to  be  incompressible  and  inviscid,  and  mathematical  models 
were  developed  to  allow  maximum  use  of  empirical  data.  The 
amenability  of  pressure  distribution  data  to  empirical  description,  as 
described  in  Section  2,  provides  some  degree  of  flexibility  in  combin¬ 
ing  analytical  and  empirical  methods. 

The  models  developed  may  be  divided  into  two  general  categories: 
phenomenological  and  pressure.  The  former  introduce  in  some  form 
the  gross  physical  effects  that  the  jet  may  be  expected  to  induce  on  the 
surrounding  stream,  such  as  blockage  and  entrainment.  The  latter 
postulate  a  flow  field  intended  only  to  give  the  proper  qualitative  pres¬ 
sure  distribution.  All  the  models  contain  arbitrary  constants  which 
are  adjusted  by  matching  the  resultant  pressures  to  experimental 
data,  in  most  cases  by  means  of  the  Fourier  series  representation 
described  in  Section  2.  Due  to  the  lack  of  sufficient  detailed  inter¬ 
ference  pressure  distributions  on  flat  plates  with  highly  under¬ 
expanded  jets,  the  data  of  Vogler  has  been  used  throughout.  It  is 
expected  that  the  techniques  developed  may  also  be  used  to  describe 
interference  pressures  when  the  jet  is  highly  underexpanded. 

3. 1  MODELS  PROPOSED  IN  THE  LITERATURE 

Various  attempts  have  been  made  at  theoretically  predicting  the 
interference  pressure  distribution  on  a  flat  plate  from  which  a 
transverse  jet  exhausts.  The  main,  concern  has  been  with  the  VTOL 
problem,  so  that  the  discussion  applies  to  jets  of  relatively  low 
velocity  whose  exit  pressure  is  roughly  equal  to  the  free  stream 
pressure. 

Numerous  investigators  (References  19,  32,  33  and  35)  have  attempted 
to  represent  the  interference  pressure  on  the  plate  by  the  inviscid  flow 
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about  an  infinite  solid  circular  cylinder  normal  to  the  cross  flow. 

While  for  some  velocity  ratios  the  agreement  is  not  unreasonable  on 
the  windward  side  of  the  jet,  it  is  not  good  on  the  downstream  portion 
of  the  plate. 

In  Reference  33,  Bradbury  and  Wood  have  represented  the  interference 
effects  produced  by  the  jet  by  means  of  an  entrainment  and  a  blockage 
contribution  to  the  velocity  on  the  plate  surface.  It  is  assumed  that 
the  entrainment  contribution  is  axisymmetric  about  the  nozzle  center- 
line,  whereas  the  blockage  term  is  not.  Entrainment  is  assumed  to 
vary  with  the  ratio  (U^/Uj),  whereas  the  blockage  term  is  taken  to  be 
independent  of  this  ratio.  Bradbury  and  Wood  show  that  along  the 
centerline  of  the  plate,  the  pressure  coefficient  should  then  have  the 
form 


Y  1  -  CP  =TuT7u“)  +  blockage  term 

t 

where  the  blockage  term  is  independent  of  (U^/Uj).  Consequently,  a 
plot  of  (1  -  Cr>)  vs  (Uj/Uco  )  should  be  a  straight  line  for  large  values 
of  (Uj/Uco).  Bradbury  and  Wood  show  that  this  is  the  case.  The 
entrainment  function  f(r)  is  calculated  by  postulating  a  sink  distribution 
along  the  axis  of  the  jet  which  will  yield  the  same  entrainment  as 
calculated  from  a  turbulent  mixing  analysis.  However,  it  is  shown 
that  the  contribution  of  this  entrainment  function  to  the  overall  pres¬ 
sure  coefficient  is  extremely  small  so  that  most  of  the  observed 
pressure  coefficient  would  have  to  come  from  the  blockage  term.  In 
Reference  33,  Bradbury  and  Wood  indicated  that  they  had  been  unable 
to  develop  such  a  blockage  term. 

Wooler,  et  al,  describe,  in  Reference  36,  a  very  complete  model  which 
they  have  formulated  for  predicting  interference  pressures  on  wings 
with  lift  fans.  Their  method  includes  a  scheme  for  predicting  jet 
trajectories,  provided  that  some  constants  can  be  adjusted  empirically. 
These  constants  fix  the  entrainment  rate  and  the  growth  of  the  jet  in 
cross-section,  and  they  are  evaluated  by  matching  theoretical  and 
experimental  jet  trajectories.  With  the  entrainment  and  blockage  of 
the  jet  thus  obtained,  Wooler,  et  al,  represent  the  effect  of  the  jet  on 
the  surrounding  fluid  by  distributing  sinks  and  doublets  along  the  jet 
trajectory.  In  particular,  the  upwash  on  the  plane  of  the  wing  is  cal¬ 
culated,  from  which  the  interference  pressure  is  then  found.  The 
authors  compare  their  theoretical  results  to  measured  pressure 
coefficients  on  a  rectangular  wing  which  they  have  tested.  They  show 
satisfactory  agreement  at  fairly  large  distances  from  the  jet  orifice, 
but  the  agreement  near  the  nozzle  is  not  good.  Figure  18  is  a 
comparison  made  from  the  results  of  Reference  36.  The  pressure 
coefficient  is  plotted  in  the  coordinates  of  Figure  10  for  6  =  90°,  by 
cross-plotting  the  chordwise  Cp  distribution  given  in  Reference  36 
for  the  midpoint  on  the  chord. 
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INTERFERENCE  PRESSURE  COEFFICIENT,  C 


Figuie  18.  Comparison  of  Data  and  Theory  of  Reference  36 


Another  model  for  predicting  interference  pressures  is  proposed  by 
Kuiper  in  Reference  9.  He  represents  the  blockage  produced  by  the 
jet  by  the  potential  flow  about  a  three-dimensional  source  in  a  uniform 
stream.  The  source  is  located  at  some  distance  from  the  exit  plane  of 
the  nozzle,  and  thisv. distance  is  adjusted  to  match  the  data  of  Ref¬ 
erence/32.  .  In  Reference  9,  however,  the  jets  are  considered  to  be 
■located  on  the  aft  end  of  the  vehicle,  so  that  Kuiper  is  only  concerned 
with  the  0  range 

90°  £  0  <  27o°. 

3.  2  PHENOMENOLOGICAL  MODELS 

Two  models  will  be  considered  under  this  heading.  Referring  to 
Figure  10,  it  is  assumed  that  the  flowfield  is  two  dimensional  in 
planes  z  =  constant,  and  that  the  free  stream  is  incompressible, 
inviscid,  and  irrotational  so  that  flow  in  the  (x,  y)  plane  obeys  the 
two-dimensional  Laplace's  equation.  The  first  model,  called  the 
doublet  model,  consists  of  the  superposition  of  a  free  stream,  source 
or  sink,  and  a  doublet.  The  doublet  attempts  to  account  for  the  block¬ 
age  effects  of  the  jet,  while  the  sink  would  represent  the  entrainment 
caused  by  turbulent  mixing.  The  second  flow  model,  calletf  the  vortex 
model,  is  the  same  as  the  first,  except  for  the  addition  of  two  counter¬ 
rotating  vortices  downstream  of  the  origin.  It  is  known  experimentally 
that  these  vortices  exist,  and  that  the  flow  bears  some  similarity  to 
separated  flow  behind  a  cylinder. 

It  has  been  shown  in  Section  2  that  a  truncated  Fourier  series  in  G 
succeeds  in  representing  the  data  with  reasonable  accuracy.  Con¬ 
sequently,  the  following  procedure  has  been  adopted  for  completely 
determining  the  models.  Once  the  number  and  types  of  singularities 
have  been  chosen,  the  complex  potential  for  the  flowfield  is  written  in 
terms  of  the  singularity  strengths.  From  this,  an  expression  for  the 
pressure  coefficient  is  obtained  and  put  in  the  form  of  a  Fourier 
series  in  0,  as  in  Equation  (2).  By  setting  each  coefficient  at  a  given 
value  of  r  equal  to  tbs  corresponding  one  obtained  from  the  data, 
equations  are  obtained  which  may  then  be  solved  simultaneously  for 
the  unknown  singularity  strengths.  This  may  be  done  for  several  values 
of  the  velocity  ratio,  using  the  data  of  Reference  32,  and  the  singularity 
strengths  obtained  as  functions  of  this  velocity  ratio.  The  benefits 
of  this  technique  are  twofold.  First,  its  use  will  allow  the  extension 
of  limited  amounts  of  data  to  other  velocity  ratios.  Second,  by 
determining  which  singularities  are  strongest,  the  dominant  factors 
in  the  flowfield  can  be  inferred.  This  information  could  be  extremely 
important  in  trying  to  relate  the  results  obtained  from  flat  plate  data 
to  an  axisymmetric  body.  The  model  will  yield  results  at  one  value 
of  r  which  are  as  close  to  the  data  as  the  Fourier  series  representa¬ 
tion  is.  The  accuracy  of  results  at  other  values  of  the  radius  will 
depend  upon  how  closely  the  model  depicts  the  actual  situation  and 
will  determine  whether  the  model  is  valid.  Also,  as  discussed  in 
Section  2,  any  model  that  accurately  reproduces  the  r  dependence  of 
the  first  two  terms  in  the  Fourier  series  will  accurately  predict 
integrated  forces  and  moments. 


3.  2.  1  Doublet  Model 

The  doublet  model  is  a  two-dimensional,  potential  flow  model  derived 
from  the  superposition  of  complex  potentials  for  a  uniform  flow,  a 
doublet,  and  a  source  (or  sink)  which  are  described  in  Reference  37. 

Let  4  denote  the  complex  variable 

C,  =  x  +  iy 

where  x  and  y  have  been  normalized  by  the  jet  diameter,  and  let  w(  £) 
denote  the  complex  potential.  Then  the  complex  potential  for  flow 
about  a  doublet  and  a  source  may  be  described  by  the  Equation 

A* 

w*  =  A*  ^  +  A*  log  £  +  -~.  (6) 


where  the  coefficients  An  are  purely  real.  If  the  model  and  actual 
flows  are  to  have  the  same  velocity  at  large  distances,  Ag-must  be 
chosen  equal  to  U^,,  and  the  pressure  coefficient 


(7) 


where  a  bar  denotes  the  complex  conjugate,  as  described  in  Ref¬ 
erence  37.  The  above,  choice  of  Aq,  however,  creates  a  problem  in 
matching  the  model  and  the  experimental  pressure  coefficients.  It 
may  be  shown  that  the  Equations  (6)  and  (7),  when  combined  and 
transformed  to  real  variables,  lead  to  a  three-term  Fourier  series 
for  the  pressure  coefficient  (up  to  and  including  a  cos  20  term). 
Equating  coefficients  in  that  series  term  by  term  to  the  experimentally 
determined  coefficients,  cn  (r),  would  lead  to  three  equations  for  the 
two  unknowns  Aj  and  This  over  specification  can  be  avoided  by 

leaving  Aq  unspecified  and  obtaining  its  value  from  a  matching  of  the 
coefficients  for  all  three  terms.  Since  the  resulting  value  of  Aq  will 
differ  from  U®  ,  the  uniform  stream  specified  in  the  model  will  not 
have  the  same  velocity  as  the  actual  free  stream  flow.  Since  the 
objective  of  the  model  is  to  match  pressures,  it  is  of  course  desirable 
that  the  model  and  actual  flows  have  the  same  free  stre~m  static 
pressure.  Consequently,  differing  free  stream  velocities  require  that 
their  stagnation  pressures  be  different.  By  allowing  different 
stagnation  pressures  in  the  flow  and  the  model,  the  pressure  at  any 
stagnation  points  which  arise  in  the  model  flow  may  be  adjusted  for 
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better  agreement  with  data.  This  may  prove  tu  ut:  nacssary,  since 
the  data  of  Vogler  shows  that,  even  on  the  windward  side,  the'  meas¬ 
ured  pressures  never  reach  the  value  of  the  free  stream  stagnation 
pressure. 

Bernoulli's  equation  for  the  model  flow  can  be  written 


,2 
U  * 


where 


|  qj'j  -magnitude  of  the  dimensional  velocity  vector 

i 

U«  =  uniform  flow  velocity  in  the  model 

The  pressure  coefficient  can  be  written  In  terms  of  the  actual  free 
stream  velocity  as 


(8) 


where  it  has  been  assumed  that  the  model  and  actual  flow  densities  are 
the  same.  Now,  in  Equation  (6)  A’,:  -  Uj.  ,  so  that  normalizing  Equa¬ 
tion  (6)  by  Uw  yields  ° 


w  =  Aq£  +  Aj  log  ?,  +  ~ 
where 


A 


o 


U 

U 


» 

CO 

00 


and  therefore  the  pressure  coefficient  is 


Substituting  Equation  (9)  into  Equation  (10)  gives 


C  = 

P 


-  A  A. 
o  1 


(9) 


(10) 


34 


This  expression  for  the  pressure  coefficient  is  now  to  be  written  as 
a  Fourier  series  in  the  azimuthal  angle  6.  It  can  be  done  directly  in 
this  case  by  the  substitution: 


4  = 


ie 

re 


4  =  re 


ie 


or  for  more  complex  models  by  the  integral  definitions,  Equa¬ 
tions  (3,  4).  After  performing  the  substitution  and  equating  the 
resulting  expression  term  by  term  to  the  series 


2 

C  (r,  Q)  s  \  c  (r)  cos  n0 

P  Z-  n 

n  =  0 


the  following  three  equations  result: 


co(r)  = 


A2  A2 
A1  .  A2 

1  +  “ T 

r  r 


c^(r)  =  - 


(11) 

(12) 


c 


2(r)  =  2 


A  A, 
o  2 

2 

r 


(13) 


These  expressions  for  the  coefficients,  cn,  are  to  be  set  equal  to 
experimentally  determined  values  at  a  fixed  value  of  r. 

The  value  of  r  chosen  for  matching  data  is  to  some  extent  arbitrary. 
The  circle  r  =  1  was  selected  because  this  is  the  region  with  the 
highest  values  of  Cp  for  which  data  were  consistently  available. 
Equations  ( 1 1 ,  12,  13)  then  become 


Aj+A^.co(l> 

(14) 

2A,(Ao  .  A2)  .  .  Cjdl 

(15) 

2Ao  A2  =  c2* 1 ' 

(16) 

These  may  be  solved  as  follows: 

Adding  Equations  (14)  and  (15)  yields 

(A1  *  A2>2  +  2A1  Ao  =  -  (co  +  clh 
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and  subtracting  Equation  (15)  from  Equation  (14)  yields 


(A  +  A_)  -  2A  A  •  -  (c  -  c,  ) 

12  1  o  o  1 


Now,  subtracting  Equation  ( 1 6 )  from  each  of  these  leads  to 

(A.  -  A_)2  +  2A  (A.  -  A_)  +  (c  +  c,  +  c-)  0 

1  Z  o  1  2  o  1  2 


(A  +  A,)  -  2A  (A.  -I-  A_)  +  (c  -  c,  +  c«)  -  0 

1  Z  012  O  1  2 

The  quadratics  in  Equations  (17,  18)  may  now  be  solved  for 
(Ai  -  A2)  and(A]  +  A2)  in  terms  of  A0  and  the  coefficients,  cn. 
Then  the  resulting  linear  equations  can  be  solved  for  Aj  and  A 2- 
The  results  are 


A  -  ±4  <  A2  -  (c  -  c 
1  2)o  o  1 


+  c2)  2+  A2  -  <co  +cj  +c2)  2| 


AZ  Ao*i  I  Ao-(co-C)  +c2}  2*  Ao'  (co  +  Cl  +c2>  2(20; 


Finally,  substitution  of  (20)  into  (16)  yields  a  single  equation  for 
A  in  the  form 


2A  ± 
o 


(co  “  C1  +  C2 


(co  hcl  +c2 


Equation  (21)  must  be  solved  numerically  for  AQ,  and  once  An  is 
known,  A]  and  A 2  tan  be  determined  by  substituting  into  Equa¬ 
tions  (19,  20).  As  the  equations  indicate,  it  is  possible  in  principle 
to  have  more  than  one  solution.  This  is  not  surprising  as  it  is  caused 
by  the  fact  that  in  matching  pressures,  a  nonlinearity  is  introduced 
through  the  use  of  Bernoulli's  equation.  For  a  specific  case,  however, 
there  has  been  no  difficulty  in  choosing  the  solution  which  is  physically 
significant. 
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For  a  velocity  ratio  (TJoo/Uj)  =  0.4.  the  coefficients  cn  calculated  from 
‘v'ogier's  data  have  the  following  values:* 


c0m  - 

-0.  619 

(22) 

Cjd)  = 

-0.  843 

(23) 

c2(l>  - 

0.  287 

'24) 

A  graphical  solution  of  Equation  (21)  is  shown  in  Figure  19.  With  the 
numerical  values  given  in  Equations  (22,  23,  24),  the  left  hand  side 
(L.  H.  S.  ),  the  right-hand  side  for  positive  sign  (R.  H.  S.  +),  and  right- 
hand  side  for  negative  sign  (R.  H,  S.  " )  of  Equation  (21)  are  shown  in 
Figure  19.  This  figure  contains  the  restriction  that  AQ  be  real,  and 
also  that  it  be  positive,  so  that  actual  and  model  free  stream 
velocities  will  be  in  the  same  direction.  Evidently,  only  one  solution 
is  possible,  and  it  lies  in  the  neighborhood  enclosed  by  the  circle 
labeled  "solution"  in  the  figure.  Further  iterations  in  this 
neighborhood  yield  the  value 

Aq  =  0.  745,  (25) 

and  substitution  in  Equations  (19,  20)  (using  the  +  sign)  then  gives 

A^  -  0.  762  (26) 

.  A2  =  0.  192  (27) 

These  results  indicate  that  the  model  free  stream  velocity  amounts  to 
approximately  3/4  of  the  actual  velocity,  so  that  the  model  stagnation 
pressure  is  Bmaller  than  that  in  the  actual  stream.  Furthermore, 
since  A]  is  positive,  the  source  is  indeed  a  source,  and  not  a  sink. 

It  is  possible  to  show  that  the  radial  velocity  at  r  »  1/2,  (which 
corresponds  to  the  rim  of  the  jet),  is  positive.  This  is  contrary  to 
what  would  be  expected  physically  since  the  jet  entrains  free  stream 
air  and  the  net  effect  should  be  that  of  a  mass  sink. 

The  real  test  of  the  model  is  the  agreement  with  data  at  values  of  r 
other  than  unity.  Comparisons  of  pressures  predicted  by  the  model 
with  data  for  seven  rays  (0  =  const.  )  are  shown  in  Figures  20  through 
26.  Examination  of  these  figures  shows  disagreement  around  the 


It  has  been  mentioned  in  Section  2  that  due  to  difficulties  with 
Vogler's  data  at  large  distances  from  the  orifice,  this  data  has  in 
some  cases  been  adjusted  so  that  C_  will  decay  to  zero.  This 
adjustment  has  not  been  made  for  tne  data  used  with  the 
phenomenological  models. 
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Fifiurs  19.  Diagram  for  Solution  of  Equation  (21),  UH/Uj  -  0.4 
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upstream  and  downstream  rays  and  good  agreement  only  near  the 
ray  9  =  "/ 2.  The  accuracy  of  other  results  obtained  from  this  model 
at  different  velocity  ratios  are  similar  to  those  described  above.  Thus, 
it  can  be  concluded  that  this  choice  of  singularities  will  not  yield  very 
accurate  pressure  distributions.  Also,  the  velocities  obtained  do  not 
behave  in  the  way  which  would  be  expected  physically. 

3.  2.  2  Vortex  Model 


The  second  phenomenological  model  studied  consists  of  a  free  stream, 
a  source  or  sink,  a  doublet,  and,  in  addition,  two  vortices  of  equal 
strength  but  opposite  sign  located  symmetrically  in  the  leeward  quad¬ 
rants.  The  latter  are  included  to  represent  the  vortex  motion  that  is 
known  to  exist  in  the  flow.  The  location  of  the  singularities  is  shown 
in  Figure  27. 


In  this  case,  model  and  free  stream  velocities  are  left  equal  since 
two  additional  unknowns  are  introduced  by  the  vortex  locations. 
Normalizing  coordinates  by  the  nozzle  diameter  dj,  and  velocities 
by  the  free  stream  value  U*  ,  the  complex  potential  for  the  flow  is 


*(4> 


A2 

t,  +  Aj  log  4  +  —  +  iA3  log 


(28) 


where  4q.  40  are  the  complex  vortex  position  vector  and  its  conjugate 


i6 


-ie 


4  =  r  e 

”o  o 


4  =  r  e 

o 


The  pressure  coefficient  is  in  this  case  given  by 


/  dw  \  /  dw  \ 
‘  \d4/\  d4  / 


(29) 


Differentiating  (28),  taking  its  complex  conjugate,  and  substituting  in 
(29)  gives  the  result 
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PRESSURE  COEFFICIENT,  Cp 


Figure  20.  Results  of  Doublet-Vortex  and  Doublet  Model  Compared  With  Data  ( 0  =  0°) 
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Figure  21 .  Results  of  Doublet-Vortex  and  Doublet  Models  Compared  with  Data  (6  -  30°) 
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Figure  22.  Rasults  of  Doublet-Vortax  and  Doublet  Modal  Compared  With  Data  (0  =  60°) 
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Figure  24.  Results  of  Doublet-Vortex  end  doublet  Models  Compared  With  Data  ( 9  =  120°) 


This  expression  contains  the  five  unknowns,  Aj,  ,  A^,  rD,  and  0O. 
These  unknowns  can  be  determined  in  several  ways.  One  would  be 
to  express  the  pressure  coefficient  as  a  five-term  cosine  series  in 
0  and  equate  the  coefficients  to  a  five-term  series  for  the  data. 

Another  way  would  be  to  use  a  three-term  series  and  equate  Ihe 
coefficients  at  two  values  of  r  (i.  e.  ,  equate  three  coefficients  at  one 
r  and  two  at  another).  A  third  option  is  to  use  physical  considera¬ 
tions  to  provide  two  equations,  and  match  a  three-term  series  to 
determine  the  other  three.  This  third  option  was  taken  as  being  more 
compatible  with  the  physical  reasoning  which  motivated  selection  of 
the  model.  The  physical  consideration  used  is  that  the  vortices  should 
remain  stationary  in  the  x,  y  plane,  and  consequently  that  the  vector 
sum  of  the  velocities  induced  at  the  location  of  a  vortex  by  all  other 
singularities  be  zero.  For  instance,  removing  the  vortex  at  l»0  from 
the  complex  potential  in  Equation  (28)  and  differentiating  the  result 
yields,  at  t,  -  C,Q, 
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Equating  real  diui  imaginary  parts  oi  this  expression  to  zero  yields 
two  relations  between  the  five  unknowns,  which  may  be  written  in 
the  following  form 


i 


t 


(31) 

(32) 


The  next  step  is  to  write  Equation  (30)  as  a  Fourier  series  in  0.  Due 
to  the  presence  of  singularities  which  are  not  at  the  origin,  the  Fourier 
series  will  in  this  case  have  an  infinite  number  of  terms,  instead  of 
terminating  as  it  did  for  the  doublet  model.  Of  these  terms,  only  the 
first  three  are  matched  to  experimental  values  of  the  pressure  coeffi¬ 
cient  at  r  -  1.  This  procedure  leads  to  three  nonlinear  equations 
which  are  to  be  solved  simultaneously  with  Equations  (31)  and  (32)  for 
the  five  unknowns  A| ,  A 3,  A3,  rQ,  and  0Q. 

The  derivation  of  the  Fourier  series  for  Equation  (30)  and  the  solution 
of  the  set  of  five  nonlinear  algebraic  equations  are  quite  complicated, 
as  described  in  Appendix  B.  There  are  a  great  number  of  possible 
solutions,  but  one  is  again  chosen  on  physical  grounds.  Calculations 
were  carried  out  using  the  coefficients,  cn,  obtained  from  Vogler's 
data  at  (U^/Uj)  =0.4  (as  for  the  doublet  model,  the  data  was  not 
adjusted).  The  numerical  results  obtained  for  this  case  are  only 
approximate,  and  further  iterations  would  have  been  necessary  to 
obtain  more  exact  numbers.  Nevertheless,  it  was  felt  that  they  were 
sufficiently  accurate  for  purposes  of  comparing  the  model  to  data  in 
order  to  determine  its  validity.  The  numerical  values  found  are  listed 
below: 


A,  =  -0.  121 

a 


(33) 


The  values  of  pi  assure  coefficient  determined  from  the  three- te  rm 
Fourier  server,  representation  of  Fvquat  ion  (30),  using  the  numerical 
values  listed  in  Equations  { i i)-  {3  () ,  are  compared  to  data  in  Fig¬ 
ures  20-26.  Evidently,  the  agreement  is  not  good.  For  much  o  f  the 
0  range,  the  agreement  is  worse  than  for  the  simpler  doublet  mode). 
Thus,  even  though  the  singularities  representing  the  jet  have  approxi¬ 
mately  the  correct  behavior  expected  on  physical  grounds,  the 
predicted  pressures  are  not  realistic. 

The  failure  of  ine  vortex  model  along  the  leeward  ray  is  especially 
significant,  because  it  points  to  a  fundamental  difficulty  of  all  inviscid 
models  which  attempt  to  simulate  the  observed  velocities  near  the 
plate  surface,  including  three-dimensional  models.  It  is  known  from 
experiment  (Reference  21),  that  immediately  behind  the  jet  the  fluid 
velocity  is  inward  toward  the  orifice.  The  presence  of  vortices  in 
the  vortex  model  was  supposed  to  account  for  this  fact,  because  of 
symmetry  and  the  condition  of  no  flow  through  the  plate,  the  velocity 
vector  in  the  plane  of  the  plate,  along  the  x  axis,  must  be  aligned 
with  the  x  axis  (see  Figure  10).  Consequently,  if  far  downstream 
perturbationu  are  to  decay  and  the  velocity  along  the  x  axis  is  to 
become  equal  to  U,,,,  the  flow  must  reverse  direction  and  have  a 
stagnation  point.  If  no  account  is  taken  of  viscous  dissipation,  the 
pressure  at  chis  stagnation  point  will  be  equal  to  the  free  stream 
stagnation  pressure,  and  the  pressure  coefficient  will  be  unity.  This 
difficulty  will  he  encountered  with  any  inviscid  model  that  attempts  to 
represent  the  inward  velocity  observed  experimentally,  no  matter  how 
complicated.  As  a  matter  of  fact,  the  failure  of  the  doublet  model  to 
approximate  the  experimental  velocities  on  the  surface  may  be  traced 
to  the  same  sources.  As  the  doublet  model  was  originally  envisioned, 
the  blockage  of  the  jet  would  be  accounted  for  by  the  doublet,  and  the 
entrainment  by  the  sink.  At  it  turned  out,  the  negative  pressure 
coefficient  on  the  leeward  side  of  the  jet  so  influenced  the  sink  strength 
that  the  sink  became  a  source. 

3.  3  PRESSURE  MODELS 

Since  the  general  shape  of  the  interference  pressure  on  a  flat  plate  is 
known,  model  flows  are  constructed  which  will  yield,  approximately, 
this  pressure  distribution.  At  this  point,  free  parameters  are  adjusted 
to  obtain  best  agreement  with  data.  No  attempt  is  made  to  qualitatively 
reproduce  the  velocities  observed  on  the  surface.  With  this  viewpoint, 
it  is  clear  that  the  results  obtained  from  the  model  should  have  a  high 
pressure  region  on  the  windward  side  of  the  jet  and  a  low  pressure 
region  on  the  leeward  side.  These  both  asymptotically  decay  to  free 
stream  pressure  at  infinity.  To  avoid  the  difficulties  encountered  with 
inviscid  models  in  the  previous  section,  the  velocity  induced  by  what¬ 
ever  represents  the  jet  must  be  in  the  same  direction  as  the  free 
stream  on  the  leeward  side. 

The  simplest  model  that  will  meet  these  criteria  is  a  source  in  a  free 
stream.  Two  models  of  this  type  have  been  developed.  The  first  uses 
the  same  assumption  of  two  dimensional  flow  in  the  x,  y  plane  (Fig¬ 
ure  10)  and  consists  of  a  two-dimensional  source  and  a  free  stream, 
adjusted  by  the  same  methods  employed  for  the  phenomenological 
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models.  II  ia  referred  to  as  the  "source  model."  The  second  model 
is  fully  three  dimensional,  and  it  consists  of  a  free  stream  super¬ 
imposed  on  axisymmetric  flow  through  an  orifice,  It  is  referred  to 
as  the  "orifice  model." 

3.  3.  I  Source  M  odel 

The  source  model  is  derived  from  the  potential  for  a  two  dimensional 
source  in  a  uniform  stream.  This  model  is  similar  to  the  doublet 
model  discussed  in  Section  3,2. 1.  It  is  simpler  than  the  doublet  model, 
but  it  yields  pressures  which  agree  better  with  data  than  those  cal¬ 
culated  from  the  doublet  model.  The  complex  potential  for  a  source 
in  a  uniform  stream  is  simply 


8  V  +A1 


log  i 


As  in  Section  3.2.  1,  the  model  and  actual  flows  are  here  assumed  to 
have  different  stagnation  pressure  and  free  stream  velocities.  The 
above  potential  has  again  been  normalized  by  the  actual  free  stream 
velocity  and  bv  the  nozzle  diameter.  Using  Equations  (10)  and  (38), 
the  pressure  coefficient  can  be  written  in  real  variables  in  the  form 


2 A  A, 
o  1 


cos  6 


Note  that  in  this  case  the  Fourier  series  contains  only  two  terms.  As 
before,  A0  and  A]  are  calculated  by  equating  (39)  term  by  term  to  a 
two-term  Fourier  series  representation  of  the  data  at  r  =  1.  For  this 
case,  however,  the  data  of  Reference  32  haB  been  adjusted  at  each 
value  of  6  so  that  Cp  will  be  zero  at  r  =  10.  For  U^/Uj  =  0.4,  the 
first  two  Fourier  coefficients  thus  obtained  are: 


-Al  =  c0(1)  = 


■  0.  725 


- 2 A  A, 
o  1 


■  0.  823 


These  yield: 


A  =  0.483 
o 


0.  851 


Comparison  of  this  model  with  experimental  data  is  shown  on 
Figures  28-34.  As  before,  these  are  plots  of  the  pressure  coefficient 
as  a  function  of  the  radial  distance  from  the  center  of  the  nozzle  for 
several  values  of  0.  As  can  be  seen,  this  model  gives  quite  good 
agreement  in  the  vicinity  of  8  =  90°.  The  agreement  near  0  =  0°  or 
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leailtJ  of  fount  and  Orifico  Modih  Compand  With  Data  (6  =  30°) 


Figure  30.  Route  of  Sourca  and  Orif  ica  Models  Compared  Wth  Oats  (0  =  60°) 


Figure  31.  Route  of  Sown  and  Orifkt  Modab  Compand  With  Data  (0 
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Figure  32.  Results  of  Sown  and  Orifice  Models  Compered  With  Data  (6  =  120°) 


•jure  34.  Results  of  Source  and  Orifice  Models  Compared  With  Data  (8  =  180°) 
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6  -  180“  Is  not  as  goou,  but  it  is  better  than  that  obtained  with  the 
doublet  model.  The  curves  are  typical  of  results  obtained  for  other 
velocity  ratios.  Curves  of  the  values  of  A0  and  Aj  as  functions  of  the 
velocity  ratio  are  shown  in  Figure  35. 

3.3.2  Orifice  Model 

The  orifice  model  consists  of  the  three-dimensional  potential  flow  derived 
by  superimposing  a  transverse,  uniform  flow  on  axisymmetric  flow 
through  an  orifice.  First  the  velocity  potential  for  flow  through  a 
circular  orifice  is  derived.  Then  the  compononts  of  the  velocity  in 
the  plane  of  the  orifice  are  determined.  A  uniform  stream  normal  to 
the  orifice  axis  of  symmetry  is  then  superimposed  on  the  orifice  flow 
by  adding  velocity  components.  Bernoulli's  equation  is  then  employed 
to  calculate  pressure  distributions  in  the  plane  of  the  orifice. 

A  derivation  of  the  potential  for  flow  through  an  orifice  is  described 
by  Lamb  in  Reference  37.  The  derivation  of  the  orifice  flow  model 
follows  basically  from  Lamb's  solution.  The  coordinate  aystom  of 
interest  is  given  by  ' 

7  =  k  cos  6  sinhq  (40) 

7  =  k  sin  4  cosh q  (41) 

where  the  variables  ~z  and  7  correspond  to  dimensional  coordinates 
as  defined  in  Figure  10.  The  constant  k  is  an  as  yet  unspecified  scale 
length.  Squaring  and  combining  Equations  (40),  (41)  yields 


-2 

r 


.2  .  2  * 

k  sin  ; 


_2 

z 


,  2  2e 
k  cos  4 


(42) 


Equation  (42)  represents  a  family  of  hyperboloids  of  revolution  (since 
it  is  independent  of  8).  The  hyperboloids  have  foci  on  the  circle 
T  =  k,  z  =  0.  The  variable  4  i®  *  parameter  that  varies  from  hyper¬ 
boloid  to  hyperboloid.  The  value  4  =  0  corresponds  to  the  line  7=0. 
The  value  |  =  n/Z  corresponds  to  the  plane  z  =  0  with  the  circular 
region  7  <  k  removed.  Values  of  |  between  0  and  v/2  correspond 
to  hyperboloids  between  these  two  limiting  cases.  Negative  values 
of  4  are  not  considered  since  7  is  always  positive.  Combining 
Equations  (40)  and  (41)  in  a  different  manner  yields 


_  2 

_ r _ 

,Z  .2 
k  cos  h  q 


+ 


.2 

_ z _ 

.  2  .  .2 

k  sin  h  r| 


(43) 


so  that  lines  q  =  constant  correspond  to  confocal  ellipsoids  of  revolution 
with  foci  on  the  circle  7  =  k.  The  coordinates  are  shown  in  Figure  36. 
From  this  figure  it  can  be  seen  that  a  solution  of  Laplace's  Equation  in 
the  (4,  q)  system  in  which  lines  4  =  constant  correspond  to  the  stream¬ 
lines  will  transform  ur.der  Equations  (40)  and  (41)  to  the  flow  out 
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of  (or  into)  an  orifice.  In  terms  cf  (£,  q)  coordinates,  the  equation  of 
continuity  becomes: 


where: 


p  =  cos  £ 
v  =  sinhq 

®  =  velocity  potential 

Any  dependence  on  0  drops  out  since  the  flow  is  axiaymmetric  about 
the  z  -  axis.  Since  the*  coordinate  system  is  orthogonal,  and  so  are 
the  streamlin2s  and  equipotential  lines,  it  follows  that  if  £  -  constant 
is  to  be  a  streamline,  then  q  -  const  should  correspond  to  an  equi- 
potential  line.  Therefore  ®  must  be  only  a  function  of  q--and 
consequently  only  a  function  of  v.  With  ®  a  function  of  v  only, 
Equation  (44)  may  be  integrated  once  to  yield: 

(  +  l'Z)  dv"  =  con8t-  r  -A  (45) 

which  can  be  integrated  again  to  yield: 
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<&  -  A  cot  '  v  (4b ) 

The  orifice  flow  velocity  components  parallel  to  the  plane  of  the  orifice 
are  given  i >y 


Equations  (40),  (41),  (46)  are  combined  with  these  equations  to  yield 
the  velocity  components 

—  (A/k)  sin  £  tan  hq  ... 

vr  r  .  2  t  .  ,2  ,  2  ,  ,  2  ' 

sin  9  sinh  q  +  cos  §  cosh  q 


(A/k)  cos 


z  sin^  £  3inh^  q  +  cos^  |  cosh^  q 


It  may  be  shown  that  these  velocities  satisfy  the  boundary  conditions 
of  the  problem.  It  may  also  be  shown  that  for  very  large  distances 
from  the  aperture  the  velocity  decays  like  the  flow  from  a  three 
dimensional  source,  as  would  be  expected. 

Although  it  is  evident  from  Equation  (42)  that  the  constant  k 
corresponds  to  the  radius  of  the  orifice,  the  constant  A  in 
Equations  (47),  (48)  is  still  to  be  determined.  It  can  be  fixed  in 
terms  of  the  volume  flow  rate  out  of  the  orifice.  Defining  V  as 
the  volume  flow  rate 


=  2  IT  J  v_  (“*,  0  ) 


and  substituting  from  Equations  (40),  (41),  (48)  gives 

a  ■  ,50) 

The  effects  of  the  free  stream  can  be  introduced  by  adding  appropriate 
components  of  Uw  vectcrially  to  the  velocity  components  on  the  plate 
surface.  Since  as  shown  in  Figure  36,  the  plate  surface  corresponds 
to  £  =  ir/2,  the  radial  velocity  component  due  to  the  orifice  flow  is, 
according  to  Equation  (47) 


-  2&l£A 


(V/2irk2) 

v  _(•?,  1 )  = - - - 

r  2  BXnhr|  coshr| 

The  transformation  Equations  (40),  (41)  become,  at  £  =  ir/2, 

F  =  k  cosh  q 

so  that  the  above  may  be  written 

(A/2trk2) 
v 

where  r'  s  r/k. 

Referring  to  the  polar  coordinates  of  Figure  10,  the  free  stream 
velocity  vector  has  the  following  components  along  the  F  and  0 
directions,  respectively: 

v  _  =  UB  cos  9 

r  or,  ” 

v0of  =  -U„  sin  e 

The  velocities  can  be  substituted  in  Bernoulli's  equation, 

-2  ^  -2 
Vr  +  Vfi 

C  =  1  -  -s 
P  (p»U^/2) 

to  yield  the  interference  pressure  coefficient  in  the  form 


(V/2iTk2U  )2  ( V / 7t k 2 U  )  cos  e 

CO  co 


In  Equation  (51),  the  pressure  coefficient  for  the  orifice  model  con¬ 
tains  two  undetermined  parameters.  These  are  the  volume  flow  V  and 
the  scale  factor  k.  In  considering  ways  to  match  Equation  (51)  to  data 
for  obtaining  V  and  k  it  was  founds  in  this  instance,  that  it  was  best 
to  assume  that  the  model  and  actual  jets  have  the  same  volume  flow 
given  by 


V  =  irk2U. 

J 


Equation  (51)  then  becomes, 
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c 

p 


i  (Ui/UJ 

•  J 

4 


r'2(r’2-l) 


(u./u  ) cos  e 

J  a 


(52) 


The  scale  factor  k,  however,  should  not  necessarily  be  set  equal  to  the 
radius  of  the  actual  jet  exit,  but  left  to  be  determined  as  an  effective 
orifice  radius. 


The  data  of  Vogler  in  Reference  32  show  positive  pressure  coefficients 
in  the  region  upstream  of  the  orifice,  which  indicate  the  predominance 
of  blockage  effects  for  velocity  ratios  U^/Uj  greater  than  0.  2.  For 
velocity  ratios  close  to  unity,  where  the  effect  of  the  jet  is  primarily 
blockage,  and  in  the  underexpanded  jet  case  of  interest  there  is  prob¬ 
ably  a  large  separated  region  behind  the  jet.  Thus,  the  jet  appears  to 
the  subsonic  free  stream  as  an  obstacle  considerably  larger  than  the 
size  of  the  orifice.  As  the  velocity  ratio  decreases  the  entrainment 
increases  and  the  apparent  obstacle  size  of  the  jet  decreases.  It  is  to 
be  expected  then  that  (k/dj)  will  decrease  as  the  velocity  ratio  (Ua,/Uj) 
decreases.  The  numerical  values  for  k  are  found  by  comparing  Equa¬ 
tion  (52)  to  the  data  of  Vogler  along  the  rays  0=0",  and  0  =  180",  At 
a  fixed  value  of  (U^/Uj),  the  scale  was  adjusted  for  an  approximate 
best  fit,  and  a  value  for  (k/dj)  deduced.  The  relationship  between 
(k/dj)  and  (Uo/Uj)  which  yields  consistently  good  results  in  the  appli¬ 
cation  of  Equation  (52)  is  shown  in  Figure  35.  The  curve  may  also  be 
represented  empirically  by  the  relation 

+  12,  153/77^  -  10.  813 

\Ui/ 


~  =  0.  #76  -  1.  874 


© 


A  comparison  of  Equation  (52)  with  the  data  of  Vogler  (Reference  6) 
for  the  specific  case  U^/Uj  =  0.  4  is  made  in  Figures  28-34.  It  is 

recalled  that  the  data  has  been  adjusted  at  each  value  of  0.  | 

Agreement  is  good  near  the  windward  and  leewai'd  planes  of  sym¬ 
metry,  but  not  too  good  near  0  =  ir/2.  In  this  neighborhood  the  decay 
predicted  by  Equation  (52)  is  too  fast  to  properly  represent  the  data. 

3.  4  PRESSURE  MODEL  COMPUTER  PROGRAMS 

A  computer  program  has  been  developed  which  calculates  pressure 

distribution,  interaction  forces,  and  interaction  moments  on  flat 

plates  and  cylindrical  shapes,  using  the  pressure  models  described  * 

in  Section  3.  3.  The  program  is  coded  in  the  FORTRAN  IV  pro-  ! 

gramming  language  for  use  on  the  IBM  7094  computer.  j 

i 

1 

» 
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3.  4.  1  Flat  Plate  Interaction  Forces 

The  data  analysis  and  calculations  of  pressures  by  the  computer  pro¬ 
gram  are  carried  out  for  the  plane  of  the  jet  exit.  The  program 
requires  either  a  set  of  pressure  data  in  that  plane,  or  the  values  for 
the  empirical  data  fit  Fourier  coefficients  and  singularity  strengths 
for  the  source  model. 

Three  different  methods  of  calculating  the  pressure  distributions  can 
be  accomplished  by  the  program.  The  first  is  either  an  empirical  fit 
of  the  data  which  was  input,  based  on  a  five-term  Fourier  series,  or 
the  pressure  distribution  calculated  by  the  same  resulting  equation 

2 

CP  =  2  Cn  C0S  (n6)  (53) 

n  =  o 

based  on  the  input  Fourier  coefficients,  cn.  The  second  alternative 
is  to  calculate  the  pressure  distribution  by  the  source  and  uniform 
stream  model,  as  given  by  Equation  (39).  The  third  alternative  is  to 
calculate  the  pressure  distribution  based  on  the  orifice  flow  in  a  uni¬ 
form  stream  model  according  to  Equation  (52). 

After  the  pressure  distribution  is  calculated,  it  is  integrated  in  the 
plane  of  the  jet  exit  to  yield  interaction  force  and  moment  coefficients. 
The  integration  of  the  input  pressure  distributions  is  accomplished 
numerically  in  the  coordinate  system  shown  in  Figure  10.  In  that 
coordinate  system,  the  numerical  integration  scheme  is  given  by 


% 


CN  -  S7  II  Vri’  V  •  *Ai! 

R  i  i 


(54) 


and 


CM.  =  I  2  riCp'V  V  c°< ! V 


(55) 


i  l 


where 

AA. 


(56) 


The  pressure  integration  for  the  source  model  can  be  accomplished 
analytically  by  integrating  Equation  (39h  The  orifice  flow  model 
interference  pressures  are  integrated  numerically  to  smooth  effects 
of  the  singularity  at  r'  =  1. 
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3.4.2  Interaction  Forces  on  a  Cylinder 


'T'V,  «  1  -.  4  **„<*.,.  «-. ^  T  T  •>.«.*« 
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vide  pressure  distributions  only  in  the  plane  of  the  jet  exit.  However, 
an  approximation  and  coordinate  transformation  have  been  introduced 
in  order  to  calculate  pressures  on  a  cylindrical  body  and  integrate 
them.  The  pressure  distribution  on  the  cylinder  is  approximated  by 
wrapping  the  plane  of  the  jet  exit  into  a  cylinder.  This  transformation 
is  made  in  such  a  manner  as  to  maintain  constant  distance  on  the  sur¬ 
face  between  the  jet  exit  and  the  point  (S,  6)  in  the  plane  located  by  the 
polar  angle  0  as  shown  in  Figure  37.  The  transformation  is  given  by 
the  equations 


x-x.  =  n  sin  8 
J  1 

y  =  q  cos  0 


2  •  2q 

q  sin  6 


where  the  distance  q  is  determined  by  numerically  evaluating  the 
integral 


'  > 


rr) 

«•  =  / 
Jo 


rn  r_2  2  2  2C  1  J^2 

I  k  -  n  sin  a  cos  0  ,  , 

- 2  - 2 - 2 -  dT1 

'0  R  -  n'  sin  0 


to  determine  the  upper  limit, 

The  pressure  coefficient  at  the  point  (S,  0)  on  the  cylinder  surface  is 
given  by  the  pressure  coefficient  in  the  plane  of  the  jet  exit,  as  the 
average  value 

cPi,  •  i  [CP  "W  *  W  +  c,"w  'M> 

+  V6i+1'  r<]  (“ 


The  associated  area  increment,  AA..,  is  the  projection  in  the 
x,  y  plane  shown  in  Figure  37. 

In  the  program,  the  pressure  coefficients  given  by  Equation  (61)  are 
integrated  numerically  to  yield  the  interaction  force  and  moment 
coefficients 


:»4  I  1  c?u AA„ 


CM  '  SR  LR  2  S  CPil  AAi»  |xj  ‘  Xcg  +  xl»*  <(,3) 

where  the  moment  is  referred  to  the  point  x 

eg. 

The  output  of  the  program  includes  the  coefficients  Cft  and  Cm  as  well 
as  the  pressure  distribution  on  the  cylinder. 
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Section  4 


CONTROL  EFFECTIVENESS  PREDICTION 
FOR  SUPERSONIC  FLIGHT 


One  major  objective  of  the  study  reported  here  has  been  to  develop  a 
computer  program  for  predicting  JI  control  effectiveness  for 
axisymmetric  missiles  in  supersonic  flight.  The  analysis  method 
employed  in  the  computer  program  which  has  been  developed  is  based 
on  the  equivalent  solid  obstacle  in  inviscid  flow  analogy.  The  analogy 
and  its  basis  are  described  in  detail  in  Reference  19.  In  this  section, 
improvements  and  increases  in  capability  and  flexibility  of  the  program 
are  described.  A  description  of  the  program  and  instructions  concern¬ 
ing  its  use  are  contained  in  Appendix  C. 

4.  1  GENERAL  DESCRIPTION  OF  THE  PROGRAM 

The  final  version  of  the  equivalent  solid  obstacle  program  calculates 

*  static  stability  derivatives  and  JI  amplification  factors  for  axisym- 
metric  missiles  with  circular  lateral  jets.  The  jet  location  on  the 
vehicle  is  arbitrary  and  all  parameters  are  calculated  as  functions  of 
angle  of  attack. 

The  equivalent  solid-obstacle  analogy  was  described  in  Reference  19 
and  is  based  on  a  momentum  balance  criterion  that  is  independent  of 
viscous  or  geometric  effects.  It  is  required  that  the  free  stream  exert 
a  drag  force  on  the  jet  plume  as  it  accelerates  the  jet  fluid  downstream. 
The  equivalent  circular  cross-sectional  area  of  the  plume  is  calculated 
by  the  method  of  Reference  19  and  the  jet  plume  is  replaced  by  a  solid 
obstacle,  a  hemisphere -cylinder  having  the  same  frontal  area. 

The  hey  assumption  involved  in  the  equivalent-body  analogy  is  that  the 
jet  plume  can  be  replaced  by  a  solid  obstacle.  It  is  further  assumed 
that  the  shock- wave  pattern  caused  by  the  equivalent  body  alone  is 
unaltered  by  the  presence  of  the  vehicle  surface,  and  the  pressures  on 
the  vehicle  surface  are  altered  by  a  factor  equal  to  the  pressure  ratio 
at  the  corresponding  point  in  the  equivalent-body  flow  field.  The 
analysis  of  Reference  19  showed  satisfactory  agreement  of  shock 
shapes  caused  by  a  jet  and  an  equivalent  hemisphere- cylinder  aligned 

*  with  the  free  stream  when  these  analyses  were  restricted  to  jets 
exhausting  from  flat  plates.  The  method  has  now  been  extended  to 
include  arbitrary  bodies  of  revolution  at  angle  of  attack.  It  is  not  to  be 
expected  that  the  shock  shape  or  details  of  die  pressure  field  behind  the 

*  shock  wave  will  match  experimental  data,  but  the  integrated  value  of 
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the  Interaction  force  should  be  related  to  the  wave  drag  of  the  equiva¬ 
lent  body  if  the  flow  field  is  indeed  shock  dominated.  Empirical  data 
can  be  introduced  to  shape  and  scale  the  geometry  of  the  disturbed 
region. 

Although  the  analysis  described  herein  is  based  entirely  on  nonviscous 
aerodynamics,  an  effort  has  been  made  to  simulate  the  effects  of 
boundary-layer  separation  resulting  from  the  impingement  on  the 
vehicle  surface  of  the  bow  shock  wave  caused  by  the  jet  plume.  The 
resulting  pressure  distributions  produce  more  realistic  control 
moment  increments  for  a  given  interaction  force.  The  secondary 
viscous  effects  altering  the  pressure  distribution  in  regions  of  high 
velocity  gradients  have  been  neglected. 

The  equivalent  solid  obstacle  analogy  has  been  combined  with  an 
existing  surface  pressure  integration  technique  currently  operational 
at  MDAC-WD.  The  integration  technique,  described  by  Timmer  and 
Stokes  in  Reference  38,  is  based  on  local  inclination  pressure  laws 
and  is  used  to  predict  and  integrate  surface  pressures  on  bodies  of 
revolution  at  angle  of  attack.  On  that  portion  of  the  vehicle  in  the 
region  of  influence  of  the  jet,  the  surface  pressures  are  multiplied  by 
pressure  ratios  determined  independently  by  a  method  of  character¬ 
istics  analysis  of  the  axisymmetrlc  equivalent  body  flow  field. 

The  complete  analysis  method  has  been  automated  and  is  currently 
available  in  the  form  of  a  MDAC-WD  FORTRAN  IV  computer  program. 
Equivalent  body  flow  fields  (i.  e.  ,  pressure  distributions  and  shock 
shapes)  for  a  unit  hemisphere-cylinder  have  been  calculated.  These 
flow  fields  were  calculated  for  the  local  undisturbed  (jet  off)  Mach 
numbers  {Mj}  at  the  jet  location  on  the  vehicle  surface.  The  jet 
penetration  height  is  then  calculated  and  the  equivalent  body  flow  field 
is  scaled  accordingly.  Finally,  the  pressures  are  integrated  over  the 
vehicle  surface  taking  into  account  the  presence  of  the  jet.  In  this 
manner,  the  angle  of  attack  variation  of  force  and  moment  amplifica¬ 
tion  factors  and  vehicle  aerodynamic  coefficients  can  be  determined 
for  any  combination  of  jet  location  and  jet  pressure  ratio. 

In  this  section  the  application  of  the  equivalent  body  analogy  to 
simulate  the  presence  of  control  jets  on  bodies  of  revolution  at  angle 
of  attack  is  described.  Appendix  C  contains  detailed  flow  charts  and 
specific  instructions  fox  using  the  computer  program, 

4.  2  VEHICLE  GEOMETRY 

The  vehicle  surface  is  described  with  respect  to  body-fixed  axes 
(x,  y,  z)  with  the  origin  at  the  nose  and  the  positive  x-axis  as  the  axis  of 
symmetry.  The  vehicle  may  be  made  up  of  one  to  eight  components 
which  are  described  in  the  x-z  plane  as  straight  lines,  circular  arcs, 
or  arbitrary  curves.  Each  component  is  then  subdivided  into  eight 
patches,  each  subtending  a  45*  angle  on  the  surface.  A  16  rectangle- 
per-patch  integration  mesh  is  constructed  on  all  patches  upstream 


■  •,  ~  *  -‘'-‘AV- 


of  the  jet.  Aft  of  and  including  the  patch  on  which  the  jet  is  located 
the  mesh  fineness  is  increased  to  64  rectangles  per  patch.  As  indi¬ 
cated  in  Figure  38,  the  free  stream  velocity  vector  is  specified  to  be 
in  the  x-z  plane,  so  consideration  of  the  half-space  y  <  0  (and  there¬ 
fore  only  4  patches  per  component)  is  sufficient  for  vehicle  geometry 
considerations. 

4.  3  ANGLK-OF1* ATTACK  DETERMINATION 

The  operational  method  of  the  program  requires  that  calculations  only 
be  made  for  vehicle  angles  of  attack  which  correspond  to  specific  local 
Mach  numbers  at  the  jet  location,  with  the  jet  off.  These  local  Mach 
numbers  correspond  to  those  for  which  equivalent  solid  obstacle  pres¬ 
sure  distributions,  as  obtained  from  method  of  characteristic  calcula¬ 
tions,  are  stored  on  magnetic  tapes. 

The  determination  of  the  vehicle  angles  of  attack  <*•  which  produce 
the  specified  local  Mach  numbers  M{  at  the  jet  location  requires 
a  numerical  solution  of  the  isentropic  flow  relations.  Different 
methods  of  calculation  are  used  depending  on  whether  relatively 
sharp  or  blunt  nosed  vehicles  are  considered. 

4.  3.  1  Blunt-Nosed  Vehicles 

The  basic  assumption  involved  in  determining  the  angle  of  attack  for 
a  specific  local  Mach  number  on  blunt  vehicles  is  that  th  fluid  wetting 
the  vehicle,  at  the  jet  location,  passed  through  a  normal  shock  at  the 
nose  and  expanded  isentropically  to  the  local  Mach  number  Mj.  The 
local  static -to- stagnation  pressure  is  given  by  tho  isentropic  relation 


-i.  ( 
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Combining  this  with  the  Rayleigh  pitot  formula. 
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the  static  pressure  ratio 
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Figure  38.  Coordinate  Syitem  For  Vahicla  Qaomatiy 


can  be  determined.  The  vehicle  angle  of  attack  that  produces  this 
pressure  is  determined  using  one  of  the  following  local  inclination 
laws: 

Tangent  Cone  (windward  surfaces) 


F=  1+^-K2 


1  + 

lx± 

(  Y  - 

1  )  K*  +  2 

1  )  k}  +  2 

* 

Prandtl-Meyer  (lee  side  hypersonic  small  disturbance) 


(67) 
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P  3  t1  +XtlK}  (68) 
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In  the  above  expressions, 
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K  =  v Mat  -  1  sin  a, 


where  »T  is  the  local  angle  of  attack  and  is  givenbythe  scalar  product 

Lj 


sin  aL  =  ^ 


where  Is  the  free-stream  velocity  unit  vector  and  n^  is  the 

unit  outer  normal  at  the  jet  location,  ~~ 


n  =  n  i+n  j+n  k 
x  — y  z  — * 


If  the  jet  is  on  the  lee  side,  the  Prandtl -Meyer  equation  is  used  and 
can  be  solved  for  K  explicitly.  Then,  since 


Ya  =  _^sin  Q  +  k  cos  o 


the  appropriate  angle'  of  attack  is  the  root  of  the  equation 


T  ( a  )  5  n  sin  a  +  n  cos  a  + 
X  z 


-  1 
QO 


The  root  can  be  found  by  numerical  solution  in  the  interval 


ir  rr 

T  c  a  <  T 


If,  however,  the  jet  is  on  the  windward  side,  the  tangent  cone  equa¬ 
tion  must  be  used.  Since  this  cannot  be  solved  explicitly  for  K,  an 
iterative  method  is  used.  In  the  program,  this  method  is  used  to 
find  the  root  of  the  function: 


G  (k)  S  I  •  P  +  TK  1  + 
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where  k  s  K^.  Using  the  value  of  K  thus  obtained,  the  angle  of 
attack  is  found  using  the  same  method  as  described  above  for  the 
lee  aide  case. 


4.  3.  2  Sharp-Nosed  Vehicles 

The  problem  of  determining  the  angle  of  attack  for  a  specific  local 

Mach  number  on  sharp  nosed  vehicles  is  slightly  more  complicated 

because  of  the  presence  of  the  attached  shock  on  the  nose  of  the 

vehicle.  In  the  case  of  a  blunt  nosed  vehicle  (Subsection  4.  3.  1), 

calculation  of  the  pressure  ratio  P  was  straightforward  and  independent 

of  the  angle  of  attack.  For  the  sharp  nosed  vehicle,  the  total 

pressure  on  the  vehicle  changes  as  the  angle  of  attack  (  and,  hence, 

the  strength  of  the  attached  shock)  changes.  It  is  assumed  that  the 

fluid  wetting  the  vehicle  surface  is  that  fluid  which  passed  through  the 

oblique  shock  at  the  nose  on  the  windward  ray.  Thus,  the  total 

pressure  at  the  jet  location  is  determined  by  the  strength  of  the  nose 

shock  on  the  windward  ray.  It  is  further  assumed  that  the  flow  at  the 

nose  is  conical.  When  the  vehicle  is  confined  to  small  angles  of  attack 

and  the  flow  near  the  sharp  nose  (half  angle  =  6  )  is  assumed  to  be 

conical,  then  the  shock  strength  on  the  windward  ray  is  approximately 

equal  to  the  strength  of  the  shock  produced  by  a  cone  of  half  angle 

6=6+  |a|  at  zero  angle  of  attack, 
c 


An  iterative  procedure  is  required  to  determine  the  angle  of  attack  in 
this  case.  It  is  reduced  in  the  current  analysis  to  a  binary-chop 
method  of  finding  the  root  of  the  equation 


H  (  or  )  lPj(«).P2(»)i  -  0.  3  <  a  <  0.  3 


(74) 


where  | or  I  =  0.  3  radians  was  arbitrarily  chosen  to  be  the  upper 
limit  on  "small  angles  of  attack".  The  functions  Pj(a)  and  F*2(a)  are 
the  static  pressure  ratios  determined  by  shock- expansion  theory  and 
a  local  inclination  law.  Ha/ing  guessed  an  initial  value  of  a,  the 
functions  Pi  (a)  and  P2(a)  are  calculated  and  the  function  H(a)  is 
evaluated.  It  can  be  shown  that  H  is  a  monotonic  function  of  a,  so 
If  the  function  H  has  the  same  sign  at  a  ±  0.  3  radians,  there  is  no 
root  in  this  interval.  In  that  case,  no  pressure  integration  is  carried 
out  for  the  angle  of  attack  corresponding  to  the  local  Mach  number 
Mi  at  the  jet  location. 

To  calculate  the  function  Pi  (a)  for  a  particular  angle  of  attack,  the 
windward-ray  shock  angle  must  first  be  determined.  An  effective 
cone  half  angle  0C  =  6  +  |a|  is  calculated,  and  the  shock  angle  08 
is  given  by  the  following  relation  from  Reference  39 


'  V 


flin  2  ®8  =  772  +  ~T  S2  +  81  ^  2  8C 


(g2  -  gj  Bln  2  6c)  2 


'  (S3  -  f?l  )  sin20c 


where: 


V  +  1 

gi  =-r- 


g2Sl‘^ 


g3 s  Y  [1  +  nr 

Moo 


Using  this  shock  angle,  the  static  pressure  ratio  Pj  (  a  )  corresponding 
lo  the  local  Mach  number  is  given  by 


Pj  (  Of  )  *  1  +  jj —  Mi 


2  Y  M2  sin2  0  -Y  +  1 
oe  s 


y?T 


(Y+  1)M*  sin20s[(Y  -  1)M*  +  2] 
2  [  (  Y  -  D  Mj  sin2  0s  +  2  j 


The  function  P,  {  «  )  is  now  obtained  from  a  local  inclination  law.  On 
the  windward  side,  the  pressure  is  obtained  directly  from  the  tangent 
cone  law 
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(77) 


where  K  is  known  for  a  given  Of  from  the  relation 


K  =  - 


(  n  sin  rv  +  n  cos  a  ) 
'  x  z  ' 


(78) 


On  the  leeward  side,  Prandtl -Meyer  hypersonic  small  disturbance 
theory  is  used  to  yield 


P2 ( «  )  s 


2  V 
Y-  1 


(79) 


The  above  analysis  was  carried  out  for  sharp  nosed  vehicles  assuming 
an  attached  shock.  To  as-sure  consistency,  a  check  is  performed  to 
verify  that  the  shock  is  attached.  The  new  equivalent  cone  angle  is 

e  =  6  +  |o| 

c 


For  a  given  free- stream  Mach  number  Mn  ,  the  shock  is  attached  on 
cones  of  half  angles  8  <  8...  The  angle  8^.  is  given  in  Reference  39 

by  the  equation 


(80) 


If  the  shock  is  detached,  the  windward  ray  streamline  is  considered  to 
have  passed  thr^  -gh  a  normal  shock. 
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4.  4  EQUIVALENT  BODY  SIZING 

The  procedure  used  in  sizing  the  equivalent  body  is  exactly  the  same  as 
that  presented  in  Reference  19.  That  is,  the  radius  of  the  hemirphere- 
cylinder  is  assumed  to  be  scaled  by  one-half  the  jet-penetration 
height,  h  .  Tho  height  hj  is  given  by  the  expression 


4 


i 


Q 

where  Xj  .  the  equivalent  body  drag  coefficient  at  the  local  Mach 
number  Mj,  is  given  by 


(82) 


and  the  jet  exit  velocity  ratio  is  given  by 


4 


’  (V  +  l)Mj2 
2  +  (  V  -  1)  Mj2 


(83) 


The  diameter  of  the  equivalent  obstacle,  2  S,  is  scaled  to  be  approxi¬ 
mately  equal  to  the  penetration  height  when  the  penetration  height  is 
five  times  larger  than  the  vehicle  diameter  at  the  jet.  When  the  jet 
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penetration  height  is  very  small  compared  to  the  radius  of  curvature  of 
the  vehicle  surface  (i.  e. ,  approaching  the  flat-plate  case),  the  drag  of 
the  equivalent  obstacle  is  assumed  to  be  associated  with  half  of  an 
axisymmetric  shock.  Then  the  equivalent  obstacle  radius  ie  adjusted 
eo  that  the  cross- sectional  area  of  the  obstacle  is  half  that  for  the 
diameter  equal  to  the  penetration  height.  Between  these  two  extremes, 
the  ratio  between  the  cross-sectional  area  of  the  equivalent  obstacle 
and  the  vehicle  cross-sectional  area  at  the  jet,  ,  is  scaled 
exponentially.  The  equation  employed  lor  the  adjustment  is 


where  A*  =  xli  /4  and  4 


0. 00736. 


4.  5  VEHICLE  SURFACE  PRESSURE  INTEGRATION 

As  local  inclination  pressure  laws  are  used,  the  pressure  on  a  partic¬ 
ular  mesh  element  outside  the  equivalent  obstacle  Bhock,  depends  only 
on  its  orientation  with  respect  to  the  free-stream  velocity  vector. 
However,  if  the  element  lies  within  the  interaction  region,  its  pres¬ 
sure  is  multiplied  by  a  static  pressure  ratio  associated  with  the 
corresponding  point  ir.  the  equivalent  body  flow  field.  To  check  whether 
a  mesh  element  with  coordinates  (x,  y,  z)  lies  within  the  interaction 
region,  the  point  must  first  be  transformed  into  the  coordinate  system 
of  the  equivalent  body. 

4.5.1  Equivalent  Obstacle  Coordinate  System 

In  the  equivalent  body  analogy,  the  hemisphere -cylinde r  is  assumed  to 
lie  parallel  to  the  local  flow  velocity  vector.  Since  the  equivalent 
body  flow  field  data  are  specified  in  an  axisymmetric  coordinate  sys¬ 
tem  X’=R'  with  the  origin  at  the  nose,  it  is  appropriate  to  place  the 
origin  of  this  coordinate  system  at  a  point  (x0>y0>z0)  out  a  distance 
=1/2  hj  along  the  unit  outer  normal  (nx.  n y,  nz)  from  the  jet  loca¬ 
tion  (xj,  yj,  zj),  as  shown  in  Figure  39.  The  X'  axis  is  aligned  with  the 
local  flow  velocity  vector.  Since  the  equivalent  body  flow  field  corres¬ 
ponding  to  a  local  Mach  number  Mj  is  based  on  a  hemisphere -cylinder 
of  unit  radius,  all  coordinates  (x,  y,  z)  are  divided  by  the  scale  factor 
Sj  so  the  X'-R'  coordinate  system  is  compatible  with  the  equivalent 
body  coordinate  system. 
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fhe  appropriate  transformation  is  derived  by  first  expressing  X'  and 
R'  in  terms  of  the  point  (xj,  y j ,  z±)  on  the  X*  axis  nearest  the  desired 
point  (x,  y,  z) 


Figure  39.  Coordinate  System  For  Equfcaitnt  Body  FlowfWd 
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and  second,  that  the  line  from  (x,  y,  z)  to  (xj,  yj,  zj)  is  normal  to 
the  X1  axis 


[(x-y,).  (y-y,),  (z-zj)]  .  [v  uy,  uz]  =  o 


(87) 


Solving  these  equations  simultaneously  gives 

u  u  (x-x  Ha  y  (y-y  \  +( u  2  +  u.j  +  u  2  z 
_  .  x  z  V  o/  y  7  z  v  1  °f  \  *  y  /  * 

*i  r~ — r~  2 

u  +  u  +  u 

X  y  Z 

uy  /  \ 

*1  =  yo+  -T  (z  -  zo) 

X  =  x  +  —  (z  -  Z  ) 

1  o  uz  \  0/ 


(88) 


(89) 


(90) 


4.5.2  inviscid  Interaction  Region 

The  interaction  region  is  defined  as  that  part  of  the  vehicle  where  the 
surface  pressure  is  affected  by  the  presence  of  the  jet.  A  purely 
inviscid  analysis  would  locate  this  region  downstream  of  the  line  where 
the  equivalent  body  shock  intersects  the  vehicle  surface.  In  that  case, 
a  considerable  portion  of  the  interaction  force  would  be  concentrated 
along  this  line  due  to  the  spike  in  the  resultant  surface  pressure  profile. 
Experimental  data  presented  in  References  16  and  17  show  a  gradual 
rise  to  a  peak  pressure  which  is  much  less  than  the  predicted  inviscid 
value  and  then  a  gradual  decay  to  the  undisturbed  pressuic.  The 
effect  of  viscosity  therefore  is  to  "smea  "  tb  spike  in  the  inviscid 
pressure  profile  by  raising  the  pressu  ..  rough  shocks  asso¬ 
ciated  with  boundary -layer  interacti'  nee  "pstream  of  the 

inviscid  shock  impingement  line.  T  id  version  of  the 

analysis,  an  option  is  included  whic  >.  ue  viscous  effects  of 

reducing  the  peak  pressure  and  enlai  ^eractlon  region  into 

an  area  upstream  of  the  shock  imping t  m. 


* 


In  the  inviscid  case,  the  interaction  region  is  defined  by  requiring  that 
the  inequalities 


be  satisfied.  The  upper  bound  X'max  is  the  known  limit  of  the 
equivalent-body  flow  field  and  R'max  the  shock  radius  at  that 
X1  =  The  shock  abscissa  X^  is  given  by  the  equation 

9 

x;  (R’)  =  ^  Cm  R'™"1  (91). 

m= 1 

In  the  program  operation,  Equation (91)  represents  an  empirical  fit  of 
the  shock  shape  from  the  method  of  characteristics  solutions. 


If  the  point  (X',  R')is  found  to  lie  within  the  interaction  region,  the 
pressure  assigned  to  the  associated  incremental  surface  area  is  deter¬ 
mined  by  scaling  the  undisturbed  pressure  at(X',  R')bythe  pressure  at 
the  nearest  point  in  the  equivalent  obstacle  pressure  distribution.  With 
the  undisturbed  (jet  off!  pressure  at  the  point  (X1,  R')  on  the  missile 
surface  denotedp(X',  R'l,  the  disturbance  pressure  at(X',R')is  given 
by 

P  (X1,  R')  =  pPc 
where  p 

P  =-=£•=  P  (M,  X',  R\  S) 
c  c 

is  the  pressure  ratio  in  the  equivalent  obstacle  pressure  distribution 
at  the  point  (X1,  R').  The  nearest  point inthe  characteristics  pressure 
distribution  is  found  by  ahunting  procedure  based  on  the  fact  that  the 
characteristics  points  are  arranged  in  order  of  ascending  x-coordinate. 
Hunting  for  the  nearest  point  may  therefore  be  confined  to  a  circle  of 
radius 


d  =  min 


I  R'-R' 

1  min 

|  x’-x; 


around  the  point  (X1,  R').  The  relatively  small  amount  of  computing 
time  required  with  this  streamlined  hunting  procedure  made  surface 
fitting  the  equivalent -body  pressure  field  unnecessary. 

An  empirical  adjustment  has  been  made  in  the  program  that  limits  the 
extent  of  the  interaction  region  in  the  vehicle  cross-sectional  plane. 
The  limitation  is  on  the  radial  angle  from  the  jet,  in  the  cross- 
sectional  plane.  No  pressure  scaling  because  of  the  disturbance  is 
done  beyond  an  angle  of  150°  away  from  the  jet.  ThL.  'djustment  was 
made  because  it  is  known  from  flow  visualization  ir  ws  d-tunnel  tests 
that  the  jet  bow  shock  dissipates  in  the  crosa-t* .’ana. 


79 


4.  5.  3  Boundary-Layer  Separation  Effects 

The  interaction  of  the  jet-induced  bow  shock  with  the  boundary  layer 
is  shown  schematically  in  Figure  40.  The  flow  geometry  and  the 
corresponding  pressure  distribution  are  shown  for  a  longitudinal  plane 
located  some  distance  laterally  from  the  jet  nozzle.  The  geometry  of 
the  lambda -type  shock  pattern  shown  is  quite  speculative  aB  there  is  no 
direct  way  of  actually  observing  these  details  in  a  three-dimensional 
experiment.  The  details  of  the  interaction  phenomena  must,  therefore, 
be  surmised  from  the  measured  pressure  distributions.  The  inter¬ 
action  of  a  swept  planar  shock  with  a  laminar  boundary  layer  is  ana¬ 
lyzed  in  some  detail  in  Reference  40,  It  is  reported  there  that  at  least 
two  plausible  flow  models  can  be  postulated  which  will  produce  the 
observed  results.  Further,  it  is  possible  that  either  type  of  shock 
structure  may  exist,  depending  upon  the  state  of  the  boundary  layer, 
the  Mach  number,  the  shock  strength,  and  other  parameters. 

As  indicated  by  comparing  the  viscous  and  inviscid  pressure  distribu¬ 
tions  illustrated  in  Figure  40,  the  primary  effect  of  the  boundary 
layer  is  to  reduce  the  peak  pressure  and  to  distribute  the  load  over  a 
greater  area.  The  initial  pressure  rise  occurs  as  a  result  of  either  a 
thickening  or  separation  of  the  boundary  layer  ahead  of  the  shock  loca¬ 
tion.  A  fully  separated  boundary  layer  with  reverse  flow  probably 
occurs  only  quite  close  to  the  jet,  where  the  local  pressure  gradients 
are  high.  Over  much  of  the  disturbed  area  only  a  local  thickening  of 
the  boundary  layer  occurs. 

The  maximum  pressure  at  the  wall  probably  occurs  just  downstream 
of  the  main  shock  at  the  location  of  the  rear  leg  of  the  lambda  shock 
pattern.  The  pressure  decay  downstream  of  the  peak  then  follows 
quite  dose  to  the  inviscid  pressure  profile,  because  there  is  no  mech¬ 
anism  in  this  region  to  support  a  large  pressure  gradient  normal  to  the 
wall. 

The  interaction  force  is  dependent  upon  the  distance  that  the  boundary 
layer  is  affected  (X8-Xl  in  Figure  40)as  well  as  the  peak  pressure. 
Numerous  attompts  have  been  made  to  correlate  this  or  related  dis¬ 
tances  with  the  pressure  rise  for  two-dimensional  separated  boundary 
layers.  For  example,  the  results  of  Needhan  and  Stollery  in  Reference 
41  are  correlated  by  Equations  13  of  Reference  19.  In  another  ana¬ 
lysis  more  applicable  to  the  present  situation,  Hakkinen,  et  al.  ,  in 
Reference  42  correlate  the  extent  of  boundary-layer  separation  pro¬ 
duced  by  two-dimensional  incident  shocks.  They  conclude  that  the 
extent  of  the  separation  should  correlate  with  the  "driving  pressure" 
which  induces  the  separation.  The  driving  pressure  is  taken  to  be  the 
difference  between  the  "final"  pressure  and  the  pressure  required  for 
incipient  separation.  In  the  two-dimensional  separation  analysis  of 
Reference  42,  the  final  pressure  Is  equal  to  the  inviscid  pressure 
behind  the  shock.  Though  their  correlation  is  not  directly  applicable 
to  the  present  analysis,  it  is  reasonable  to  expect  that  somewhat 
similar  trends  may  exist.  In  their  analysis  (which  they  verified  by 
experiment)  they  showed  that  the  separated  length  varied  almost  lin¬ 
early  with  the  difference  between  the  final  or  inviscid  pressure  and  the 


80 


plateau  pressure.  Assuming  a  similar  behavior  may  exist  for  the 

more  complex  three-dimensional  case  of  interest  here,  the  distance 

(Xg-Xji  in  Figure  40  should  increase  as  the  pressure  difference 

(p.  -  p  )  increases, 

rinv  rexp 

An  empirical  hyperbolic  curve  fit  for  the  shape  of  the  bow  shock  caused 
by  jets  exhausting  transverse  to  a  flat  plate  is  presented  in  Reference 
19.  The  equation  for  the  shock  radius  (R)  normalised  with  respect  to 
the  jet  height  (h)  is 


6.25 


(92) 


where 


and  x  is  the  distance  downstream  of  the  shock  apex.  Differentiating 
the  above  equation  provides  the  local  shock  angle  (6),  in  the  form 


P2tan20  -i+i^iSJt 

(R/hr 


(93) 


The  shock  angle,  combined  with  the  expression  describing  the  pressure 
rise  across  an  oblique  shock  in  air,  yields  an  equation  for  the  pressure 
(pinv)  immediately  behind  the  shock,  in  the  form 


(94) 


where  p*.  is  the  pressure  behind  a  normal  shock.  This  expression 
provides  the  pressure  which  would  exist  immediately  behind  the  shock 
(Pinv  in  Figure  40)  in  the  absence  of  any  viscous  effects  and  is  the 
pressure  immediately  behind  the  main  incident  shock  outside  of  the 
shear  layer  illustrated  in  Figure  40. 

The  measured  surface  peak  pressure  must  be  closely  related  to  the 
pressure  which  exists  behind  the  shock  outside  the  shear  layer. 
Therefore,  one  would  expect  that  experimental  data  should  correlate 
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Figure  40.  Bow  Shock/Boundary  Layer  Interaction  Approximation 
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with  the  same  parameter  which  governs  the  inviscid  pressure  rise. 
Experimental  values  of  the  peak  pressure  are  presented  in  Figure  41 
•  in  terms  of  the  parameters  suggested  by  the  form  of  Equation  94.  The 

inviscid  pressure  computed  from  Equation  94  is  also  presented  for 
comparison.  As  shown,  the  experimental  pressures  are  small  com¬ 
pared  to  the  predicted  inviscid  values.  The  characteristic  plateau 
,  pressure  for  a  two-dimensional  turbulent  separation  is  also  shown  on 

Figure  41  for  the  Mach  numbers  of  the  test  data  from  References  16 
?nd  43.  It  is  seen  that  the  maximum  pressure  can  exceed  the  two- 
dimensional  plateau  pressure  near  the  jet;  these  points  correspond 
to  the  "second  peak"  pressure  which  occurs  immediately  upstream  of 
the  jet.  In  general,  however,  the  pressures  are  substantially  less 
than  the  two-dimensional  values.  A  curve  fit  of  the  experimental  data 
is  provided  by 


pexp * pl\  _  0>4e-0.8(R/h)Mp'1/2 

PN'pl  / 


In  the  equivalent  obstacle  analogy  computer  program,  Equation  (95)  is 
used  to  account  for  the  effect  of  boundary  layer  separation  on  the 
interaction  control  moment.  The  inviscid  flow  pressure  distribution 
downstream  of  the  bow  shock  is  distributed  over  a  distance  extending 
from  Xj  upstream  of  the  shock  to  X2  downstream  of  the  shock,  as 
shown  in  Figure  40. 

The  distance  Xj  is  determined  by  truncating  the  inviscid  spike  in  the 
pressure  profile  at  the  experimentally  observed  value  determined 
from  Equation  (96)  and  redistributing  the  remaining  interaction  force 
into  a  region  of  constant  pressure  gradient  upstream  of  the  shock 
impingement  line.  The  value  of  the  constant  pressure  gradient,  and, 
hence,  the  distance  Xj  upstream  of  the  shock,  is  fixed  by  the  require¬ 
ment  that  the  interaction  force  due  to  the  spike 


A 


♦ 


AF,  =■£  {pinv'pexp){X2  _  Xs) 


(96) 


equal  the  area  under  the  upstream  triangle.  This  requirement  gives 
the  nondimensional  ratio 


-1 

fifty  I  pt  _! 

D  -  p  • 
rexp 


(97) 
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where  Pinv  is  the  calculated  inviscid  flow  pressure  and  pcx~  is  the 
experimentally  determined  value  as  given  by  Equation  (95)?*^  By  defini¬ 
tion  then,  any  point  with  coordinates  (X'j,  R1)  such  that 


X'  <  X'  <  X' 
1  u  s 


where  R,j  <  R^ax  ,  has  a  pressure  ratio  Pu  =  pu/pi  associated  with  it. 
The  pressure  ratio,  Pu,  is  given  by 


x»  -X* 

P  =  P  +  —2 - iL  _  P  \ 

U  eXp  X'  -X'  '  ®*P/ 

s  1 


(98) 


For  any  point  (X^,  R^,)  such  that 


x;  <  xb s  xz 


and  <  R^^,  the  associated  pressure  is  given  by  Equation  95  as 


Figure  42.  AMICOM  Ruction  Jot  Foret  Modtl 
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For  points  in  the  interaction  region  downstream  of  X2,  the  inviscid 
pressure  scaling  method  discussed  in  Subsection  4.  5.  2  is  employed. 

4.  6  EQUIVALENT  OBSTACLE  PROGRAM  RESULTS 

The  equivalent  obstacle  analogy  program  has  been  used  to  predict 
amplification  factors  for  a  variety  of  configurations  and  flow  conditions 
for  which  wind  tunnel  data  are  available.  The  only  check  on  the  validity 
of  the  analogy  is  the  accuracy  of  such  force  and  moment  data  compari¬ 
sons,  since  details  of  the  flow  are  not  simulated.  In  Figures  43  to  50, 
some  data  comparisons  are  shown  for  the  AMICOM  wind  tunnel  model 
illustrated  in  Figure  42  (for  details  of  the  model  see  Reference  31). 
These  comparisons  are  representative  of  the  general  accuracy  level 
achieved  with  the  computer  program  applied  to  other  flow  conditions 
and  missile  geometry.  Results  are  generally  better  for  aft  located  jets 
and  higher  freestream  Mach  numbers.  For  all  the  following  data  com¬ 
parisons,  the  jet  is  at  the  center  location  indicated  in  Figure  42.  All 
moments  are  referred  to  the  nose  and  the  jet  is  on  the  leeside  of  the 
model  at  positive  angle  of  attack. 

The  basic  normal  force  and  pitching  moment  coefficients  are  shown  ver¬ 
sus  angle  of  attack  for  M^  3.  0  in  Figures  43  and  44.  For  these  coeffi¬ 
cients,  the  variation  with  pressure  ratio  is  due  principally  to  increasing 
jet  thrust.  The  basic  accuracy  level  of  the  jet- off  aerodynamics  pre¬ 
dictions  at  low  Mach  numbers  is  shown  in  these  figures.  Better  accur¬ 
acy  is  achieved  at  higher  Mach  numbers. 

In  Figures  45  and  46,  force  and  moment  amplification  factors  are  shown 
as  functions  of  angle  of  attack  and  Mach  number.  As  indicated,  the 
accuracy  of  the  prediction  method  is  worse  for  lower  pressure  ratios, 
particularly  when  interaction  forces  are  negative.  This  is  believed  due 
to  lack  of  proper  compensation  for  low  pressures  aft  of  the  jet  on  the 
missile  surface.  The  angle  of  attack  effects  are  difficult  to  generalize 
because  of  limited  extent  of  the  data. 

Prediction  of  force  and  moment  amplification  factor  sensitivities  to  jet 
thrust  are  shown  in  Figures  47  and  48.  Again,  the  low  accuracy  level 
at  low  thrust  can  be  observed. 

At  a  constant  pressure  ratio,  the  sensitivities  of  amplification  factors 
to  Mach  number  are  shown  in  Figures  49  and  50.  The  geneval  accuracy 
level  indicated  for  this  relatively  high  pressure  ratio  is  believed  to  be 
as  good  as  can  be  expected  from  the  analogy,  for  this  type  of  configura¬ 
tion  (1.  e.  ,  forward  jets). 


It  is  believed  by  the  authors  that  progress  must  be  tmde  in  basic 
understanding  of  three  dimensional  effects  due  to  J1  bow  shock  dissipa¬ 
tion  and  viscous  effects  downstream  of  the  jet  plume  before  significant 
increases  in  the  accuracy  level  of  amplification  factor  predictions  can 
be  made.  In  the  interim,  the  equivalent  obstacle  analogy  method 
appears  to  offer  as  accurate  a  prediction  scheme  as  is  available. 
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Section  5 

JET-FIN  INTERFERENCE  EFFECTS 


Transverse  jet-induced  interference  effects  in  subsonic  or  supersonic 
mainstreams  have  been  divided  into  the  two  general  categories,  near¬ 
field  effects  and  fin-interference  effects.  Near-field  effects  are  con¬ 
fined  to  the  neighborhood  of  the  nozzle  and  represent  the  direct  JI 
effect  in  amplifying  or  degrading  the  force  and  moment  which  the  jet 
thrust  would  produce  in  still  air.  In  the  category  of  fin-interference 
effects  are  phenomena  assumed  to  occur  many  nozzle  diameters  down¬ 
stream,  where  the  deflected  jet  may  interact  with  aerodynamic  control 
fins.  The  discussions  in  Sections  2  through  4  pertain  tp  the' near-field 
category  of  JI  effects. 

The  calculation  of  fin  interference  effects  has  been  restricted  to 
vehicle  configurations  where  the  fins  lie  at  a  considerable  distance 
downstream  of  the  transverse  jet  nozzle  exit.  It  is  assumed  that  at 
these  downstream  distances,  the  jet  is  almost  aligned  with  the  free 
stream.  It  is  further  assumed  that  the  jet  has  been  reduced  to  two 
counterrotating  vortices  and  a  region  of  turbulence  whose  average  axial 
velocity  is  almost  equal  to  the  free  stream  velocity.  Models  have  been 
developed  that  predict  the  variation  of  vortex  strength  with  distance 
from  the  nozzle  exit  and  other  jet  and  free  stream  parameters, for  both 
subsonic  and  supersonic  free-stream  Mach  numbers.  With  the  strengths 
and  location  of  the  jet-induced  vortices  known,  interference  forces  and 
moments  due  to  far  downstream  interactions  between  the  jet  and 
fins  can  be  calculated.  The  methods  developed  in  this  study  are  appli¬ 
cable  to  cruciform  missiles  at  arbitrary  flight  Mach  numbers  and 
attitudes. 

A  complete  solution  of  the  interference  problem  requires  a  knowledge  of 
fin-jet  interference  effects  as  well  as  of  jet-fin  interference  effects. 

That  is,  the  effects  of  the  body  and  fins  upon  the  jet  trajectory,  vortex 
strengths,  etc.  should  be  estimated.  However,  it  is  assumed  here  that 
fin-jet  interference  effects  are  small,  and  that  the  jet  behaves  at  all 
times  as  if  it  were  exhausting  into  a  uniform,  infinite  crossflow.  With 
this  restriction,  the  calculation  of  jet-fin  interference  effects  is  broken 
into  two  parts.  First;  a  semi-empirical  model  of  the  jet  in  an  infinite 
crossflow,  valid  at  large  distances  from  the  orifice,  is  developed. 
Second,  the  interference  forces  induced  by  the  jet  are  calculated. 

5.  1  GENERAL  CONSIDERATIONS 

For  subsonic  Mach  numbers,  data  such  as  those  shown  in  Figures  7 
and  8  indicate  that  the  major  interference  pressure  disturbances  are 
confined  to  a  distance  of  five  or  six  Mach  disk  heights  from  the 


centerline  of  the  nozzle.  For  supersonic  flight  Mach  numbers,  data 
such  as  those  shown  in  Figure  51  indicate  the  major  interference  nor¬ 
mal  force  occurs  less  than  ten  penetration  heights  downstream  of  the 
jet  exit.  In  Figure  51,  the  normal  force  increment  is  defined  as 


AC  =(C„)  -  (C7) 

fin  on  fin  off 


with  the  jet  and  mainstream  flow  conditions  the  same  with  fin  on  and 
fin  off.  The  configuration  is  described  in  Reference  19  or  31.  The 
conclusion  to  be  drawn  from  the  limited  extent  of  the  major  interaction 
disturbances  in  either  subsonic  or  supersonic  mainstreams  is  that  fins 
located  downstream  of  this  region  will  encounter  relatively  small  pres¬ 
sure  disturbances.  However,  the  resultant  interference  control 
moment  due  to  fin  interference  may  still  be  large. 

Some  of  the  data  from  the  AM1C0M-CAL  tests  (References  27 
and  28)  were  obtained  with  instrumented,  cruciform  rectangular 
fins  on  the  model  to  measure  interference  forces  on  the  fins. 

The  configuration  shown  in  Figure  4  corresponds  to  the  tests 
in  Reference  27,  and  for  this  case  the  sensitivity  of  the  fin 
force  and  moment  balance  was  apparently  too  small  to  detect 
the  interference  forces.  For  the  tests  described  in  Reference  28, 
however,  a  more  sensitive  balance  was  used,  and  significant 
interference  forces  were  measured. 

The  data  of  Burt  and  Dahlke  in  Reference  44  show  that,  for  a  configura¬ 
tion  with  opposed  transverse  jets,  the  strongest  fin  interference  effects 
occurred  when  the  fins  nearest  the  jet  plumes  were  placed  in  a  slightly 
asymmetric  position  relative  to  the  plume.  The  most  significant  fin 
interference  effects  appear  to  be  caused  by  the  two  counterrotating 
vortices  created  by  the  interaction  of  the  jet  and  the  cross  flow.  In 
References  45  and  46,  Dahlke  has  measured  the  strength  of  these 
vortices  at  one  station  downstream  of  the  nozzle,  based  on  flow  field 
surveys  conducted  with  a  special  probe. 

In  the  present  study,  several  assumptions  have  been  made,  based  on 
the  available  experimental  data,  to  derive  a  semi -empirical  mathe¬ 
matical  model  for  jet-fin  interference .  It  is  assumed  that,  in  the 
region  where  the  fins  lie,  the  jet  is  almost  aligned  with  the  free  stream. 
Indeed,  it  will  be  assumed  that  the  fins  lie  within  the  "vortex  zone"  of 
the  jet,  which  has  been  discussed  at  the  end  of  Subsection  2.  1. 
Restricting  the  analysis  to  small  missile  angles  of  attack  then  permits 
bringing  the  entire  fin  interference  problem  within  the  context  of  the 
slender-body  approximation  (i.e.,  of  crossflow  velocities  that  are  much 
smaller  than  the  free-stream  velocity).  The  above  restrictions  are 
satisfied  by  configurations  of  practical  interest,  as  examination  of 
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Figure  52  as  well  as  the  configurations  tested  in  References  31  and  44 
will  reveal.  Some  of  the  data  leading  to  the  assumptions  will  be  illus¬ 
trated  below. 

5.2  JET  PROPERTIES  AT  LARGE  DISTANCES  FROM  THE  NOZZLE 

It  is  known  that  a  jet  in  a  subsonic  or  supersonic  crossflow  contains  two 
counterrotating  vortices.  This  is  true  of  subsonic  (References  21,  22, 
and  23),  as  well  as  sonic  or  supersonic,  highly  underexpanded  jets 
(References  16  and  46).  For  example,  Figure  52,  which  has  been  taken 
from  Dahlke's  report  (Reference  46)  clearly  shows  the  two  vortex 
regions.  The  vectors  in  the  figure  represent  the  Mach  number  compo¬ 
nent  in  a  plane  perpendicular  to  the  body  axis.  The  circulatory  nature 
of  the  flow  is  clearly  visible.  For  subsonic  jets,  Pratte  and  Baines  in 
Reference  22  indicate  that  at  large  distances  from  the  nozzle  the  axial 
velocity  in  the  jet  is  almost  equal  to  the  free  3tream  velocity.  Further, 
the  vortices  are  effectively  convected  at  this  velocity,  while  their 
strength  decays  because  of  viscous  dissipation.  It  is  assumed  that  highly 
underexpanded  sonic  or  supersonic  jevs  exhausting  into  subsonic  or 
supersonic  streams  exhibit  similar  behavior  at  large  distances  from 
the  nozzle. 

5,  2.  1  Vortex  Strengths  in  a  Subsonic  Jet 

A  semi-empirical  model  to  predict  the  variation  in  vortex  strength  with 
distance  iB  postulated  in  this  section.  It  is  first  assumed  that  the  jet  is 
everywhere  subsonic,  but  similarities  between  subsonic  jets  and  highly 
underexpanded  jets  in  subsonic  or  supersonic  crossflows  are  formu¬ 
lated,  which  allow  the  results  obtained  for  subsonic  jets  to  be  extended 
to  the  latter  cases. 


The  vortices  are  assumed  to  be  convected  downstream  at  Ugo  (Fig¬ 
ure  53),  and  the  flow  is  analyzed  as  an  unsteady  flow  in  the  y-z  plane. 
This  is  consistent  with  the  assumption  of  crossflow  velocities  which 
are  small  compared  to  Ugo-  In  the  y-z  plane,  the  jet  is  represented  by 
two  counterrotating  vortices  located  at  ( - yQ ,  zQ)  and  (yQ,  z0)*,  and  con¬ 
nected  by  a  vortex  sheet  of  vanishing  strength  as  illustrated  in  Fig¬ 
ure  54.  The  vortex  positions  and  strengths  are  assumed  to  depend  on 
time.  From  one  instant  of  time  to  another,  impulsive  pressures  of 
different  magnitude  would  have  to  be  applied  across  the  vortex  sheet  to 
generate  the  fluid  motion  instantaneously  from  rest.  The  resulting 
impulse  T_  may  be  calculated  from  the  relation 


p  $  n  ds 


(100) 


« 


* 


t 


* 


*As  before,  these  coordinates  have  been  normalized  by  the  nozzle  exit 
diameter. 
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Figure  52.  Jot  Wako  in  Fret  Stream  Flow  Fiald  at  Mach  OS  and  Zara  Anglo  of  Attack 


where  is  the  potential  of  the  fluid  motion,  and  p  is  the  fluid  density. 
The  integral  in  Equation  (100)  is  evaluated  on  a  path  C  enclosing  both 
vortices  and  their  connecting  sheet  as  discussed  in  Reference  37.  For 
the  vortices  shown  in  Figure  54,  the  potential  may  be  written  as: 


Substituting  Equation  (1 01)  into  Equation  (1  00)  and  performing  the  inte¬ 
gration,  yields: 


X  =  k<2pryo) 


(102) 


where  k  is  a  unit  vector  in  the  z  direction.  The  rate  of  change  of 
J  with  respect  to  time  is  equal  to  the  net  force  which  must  be  apolied 
to  the  vortices  and  connecting  sheet  system  to  generate  the  fluid 
motion  instantaneously  from  rest, as  discussed  in  Reference  37.  This 
force  is  given  by 


dT  , 

IT  =  *  (2p}  I  <rV 


(103) 


In  the  present  development,  it  will  be  assumed  that  vortex  strength 
and  separation  must  vary  in  such  a  way  that  the  net  force  on  the  sys¬ 
tem  of  vortex  sheet  and  vortices  is  zero.  This  assumption  has  been 
made  by  Bryson  in  Reference  47,  in  computing  lift  forces  on  slender 
bodies  at  high  angles  of  attack.  Then  Equation  (10’')  ~equires  that 


d 

dt 


(ryQ)  =  0 


so  that 


(104) 


where  K'  is  constant  for  fixed  Ugo/Uj  and  dj.  Equator  ;i04)  implies 
that  if  the  vortices  draw  apart  their  strength  must  decrease.  If 
viscous  dissipation  is  included,  it  seems  reasonable  that  the  vortex 
strengths  should  decay;  this  model,  however,  is  ir.viscid,  so  it  is 
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difficult  to  explain  what  happens  to  the  vorticity  caused  by  the  decrease 
in  vortex  strength.  In  Bryson's  model  of  separated  flow  about  a  body 
of  revolution  at  high  incidence,  vorticity  generated  in  the  boundary 
layer  about  the  body  is  assumed  to  be  "fed"  to  the  separation  vortices 
through  the  intervening  connecting  sheets.  No  such  mechanism  may  be 
used  to  account  for  vorticity  lost  in  the  present  model,  since  no  external 
boundaries  are  present.  The  only  explanation  to  account  for  the  "lost 
vorticity"  is  that  as  the  vortex  on  one  side  decays  in  strength,  a  t.mall 
amount  of  vorticity  is  carried  via  the  connecting  sheet  to  the  plane  of 
symmetry,  to  be  cancelled  there  by  the  vorticity  of  opposite  sign 
arriving  from  the  vortex  on  the  other  side. 

The  vortices  shown  in  Figure  54  will  convect  upward  at  the  velocity 
induced  by  one  vortex  at  the  location  of  the  other.  This  velocity  is  given 
by 


dz 
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r 


dt  "  4  Tf  y  d. 

yo  j 


and,  since  the  vortices  were  assumed  to  convect  downstream  at  U,*,, 


dz 

_ c 

dx 


■=  tan  p  = 


4rrUooy  d. 

o  j 


where  the  angle  p  is  defined  in  Figure  53. 
If  p  is  assumed  to  be  small  so  that 


tan  p  -  p  -  sin  p 


then 
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4  rrU 


(105) 


Given  the  dependence  of  yQ  on  £,  Equations  (1 04)and  (1 05)may  be 
used  to  calculate  the  change  in  P  and  the  vortex  trajectory. 

In  Reference  22,  Pratte  and  Baines  find  that  they  can  correlate  their 
data  for  jet  trajectory  and  thickness  by  the  use  of  variables  scaled  by 
the  jet  to  free-stream  velocity  ratio  such  as 

X  =  o'  5 
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Figure  54.  Jst-lnductd  Vorticn  in  the  Crowflow  Plane 


where  £  is  the  c  >rdinate  along  the  jet  trajectory  shown  in  Figure  53 
and 
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Pratte  and  Baines  also  find  that,  in  the  vortex  region,  the  jet  cross- 
section  grows  as  X^/3,  apparently  in  self-similar  fashion  (Refer¬ 
ence  22).  This  is  the  same  behavior  exhibited  by  turbulent  jets  in 
coaxial  external  streams,  when  the  difference  between  jet  and  free 
stream  velocities  is  small  (Reference  48).  In  the  present  case  it 
will  be  assumed  that  the  vortices  also  spread  as  X*'3,  so  that 

1 

Y  =  Y  >  ^  (106) 

o  v 


whe  re 


Y 
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a 


lOt* 


and  Yv  is  a  constant  which  has  the  same  value  for  all  velocity  ratios 
and  jet  diameters. 

Substituting  Equation  (104)  into  Equation  (105)  yields 


Written  in  terms  of  similarity  variables  X,  Y  ,  and 
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this  becomes 
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and  substitution  of  Equation  (1  06)  then  yields 
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(107) 


The  data  of  Prattc  and  Baines  show  that  ZQ  is  a  universal  function  of 
X,  and  consequently  it  follows  from  Equation  (1  07)  that 


vi  2 
K  Q- 

4ttU  d. 
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(108) 


where  K  is  a  universal  constant  which  has  the  same  value  for  all 
velocity  ratios.  Furthermore,  it  is  possible  to  integrate  Equation  (107) 
and  obtain 


(109) 
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The  behavior  predicted  by  Equation  (1 09)  for  the  jet  trajectory  is 
verified  by  the  data  of  Pratte  and  Baines,  who  measured: 


4 


z 

o 


(const.) 
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as  shown  in  Figure  5  of  Reference  22.  Finally,  writing  Equation  (104) 
in  terma  of  similarity  variables  leads  to  the  result 
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v 
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(110) 


where 


(111) 


Equations  (1 1 0)  and  (1 1 1)  state  that  the  product  of  normalized  vortex 
strengths  (T MrrUoodj)  and  velocity  ratio  is  a  universal  function  of. 
the  similarity  variable  X. 

In  summary,  the  postulated  variation  in  vortex  spacing  Y0,  leads  to  a 
vortex  strength  variation  which  can  be  used  to  predict  the  correct 
formforthe  jet  trajectory.  Indirectly,  at  least,  this  appears  to  verify 
the  proposed  relation  between  vortex  strength  and  spacing.  The  results 
derived  are  expected  to  hold  in  the  vortex  zone,  whieft  has  been  found 
to  lie  downstream,  of  the  value  X  =  5,  as  shown  in  Reference  22.  Two 
empirical  constants  have  been  introduced,  K  and  Yv.  These  will  be 
calculated  using  Dahlke's  measurements  of  vortex  strength  and  posi¬ 
tions.  Before  this  can  be  done,  it  is  necessary  to  postulate  the  equiv¬ 
alence  between  the  subsonic  and  a  sonic,  highly  underexpanded  jet. 

5.  2.  2  Equivalence  Between  Subsonic  and  Sonic  Underexpanded  Jots 

The  equivalence  between  the  subsonic  and  sonic  jets  is  formulated  in 

terms  of  a  scaling  length.  For  underexpanded  jets,  the  scale  chosen  » 

is  proportional  to  the  Mach  disk  height  or  jet  penetration  height, 

depending  on  whether  the  free  stream  is  subsonic  or  supersonic. 
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Equivalent  similarity  variables  are  defined  as  follows: 


» 


(112) 


where: 
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d^j  equivalent  subsonic  jet  diameter 
PeUe:  equivalent  subsonic  jet  mass  flux  per  unit  area 


It  is  assumed  that  the  equivalent  subsonic  jot  and  the  actual  jet  have 
the  same  mass  flux,  so  that: 


U  d 


2  T,  .  2 
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2  J  J  J 


Thus, 


(113) 


If  it  is  assumed  that 


T  =  T  , 
°.  o 
J  ® 


the  mass  flux  per  unit  area  ratio  may  be  written  in  the  form 
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and  for  a  sonic  nozzle,  this  becomes: 


U 

r«  OB 


PjUj 


ry  +  1 

rx  i  / 

p.  i 

2 

I  T 


M„ 


(114) 


All  variables  in  Equation  (114)  are  known,  but  in  Equation  (1 13)  it  is 
necessary  to  have  a  relation  between  the  equivalent  and  actual  jet 
diameters.  A  plausible  characteristic  scale  is  suggested  by  the 
behavior  of  the  internal  shock  system  in  a  highly  underexpanded  plume. 

Considering  first  the  case  of  subsonic  mainstream  Mach  numbers,  it 
is  assumed  that  the  plume  behaves  as  if  it  were  exhausting  into 
still  ah.  Then  Reference  30  shows  that  for  high  values  of  the  pres¬ 
sure  ratio  (P  =  pQ./pa,),  the  location  of  the  terminal  shock  or  Mach 
disk  is  proportional  to  Pl/2  multiplied  by  the  jet  exit  diameter. 
Reference  30  also  shows  that  the  diameter  of  the  Mach  disk  varies 
approximately  as  pl/2  multiplied  by  the  nozzle  exit  diameter.  Since 
the  equivalent  subsonic  jet  diameter  should  depend  on  the  subsonic 
conditions  which  exist  in  the  jet  downstream  of  the  Mach  disk,  it  will 
be  assumed  that 


=  b/p 


(115) 


Then,  substituting  Equations  (1 1 4)  and  (1 1 5)  into  Equation  (1 1 3)  leads  to  the 
final  relation 


3) 


(116) 


where  the  quantity  b  is  another  empirical  constant  to  be  obtained  from 
data. 


Considering  now  the  case  of  supersonic  free -stream  Mach  numbers, 
the  relation  for  the  characteristic  scale  must  be  changed.  As  dis¬ 
cussed  in  Reference  19,  the  flow  field  to  some  extent  scales  with  the 
jet  penetration  height  hs.  An  expression  for  this  quantity  is  derived 
in  Reference  19.  It  may  be  written  in  the  form 
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where  C„  denotes  the  drag  coefficient  of  the  equivalent  obstacle.  If, 
as  in  Reference  19,  it  is  assumed  that;. 
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then  the  above  becomes: 
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It  is  also  shown  in  Reference  19  that  the  Mach  disk  height  is  directly 
proportional  to  the  penetration  height  hg.  For  supersonic  free-stream 
Mach  numbers,  then,  Equation  (1 17)  suggests  a  relationship  for  the 
equivalent  jet  diameter  in  the  form 
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(118) 
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where  t  is  an  empirical  constant.  For  supersonic  fret:  streams, 
Equation  (113)  then  becomes 


(119) 


5.2.3  Evaluation  of  Empirical  Constants  and  Comparisons  with  Data 

Reference  46  contains  data  on  vortex  strengths  and  positions  for  a 
sonic,  highly  underexpanded  jet  exhausting  from  an  ogive-cylinder  for 
Ma  =0.9  and  1.2.  The  data  were  obtained  by  surveying  the  flow  field 
at  a  fixed  station  downstream  of  the  jet  nozzle  and  varying  (Po./Pg,)  for 
each  value  of  Mg,.  In  using  these  results,  the  effects  of  the  body  will 
be  neglected  and  it  will  be  assumed  that  the  jet  behaves  as  if  it  were 
exhausting  into  an  infinite  stream. 

To  calculate  the  values  of  6  and  £  defined  in  Subsection  5.2.2,  some  of 
the  results  obtained  by  Pratte  and  Baines  rn  Reference  22  will  be  used. 

In  particular,  if  it  is  assumed  that  the  vortices  lie  in  the  same  plane 
as  the  jet  centerline,  then  the  empirical  relation  in  Figure  4  of  Ref¬ 
erence  22  yields  the  vortex  height  as  a  function  of  the  distance  X  from 
the  nozzle  centerline  (Figure  53)  in  the  form 


Z 

o 


=  (1.76) 


x(0.28) 


(120) 


Using  Equation  (112),  Equation  (120)  may  be  written  in  the  form 


(*) 


0.3  89 


pa1- 389 


Then  substituting  for  {re/de)  from  Equation  (116)  in  the  subsonic  case, 
and  from  Equation  (119)  in  the  supersonic  case: 


6i  -  (0. 575) 

M  =0.9 

OO 


,-,0.389 

(x) 

,1.389 

(zo) 


(121) 


106 


=  (0.492) 


(122) 


'M  =1.2 

OO 


(s)°-389 


<jp 


The  above  relations  have  been  written  for  the  specific  cases  M*  =0.9 
and  =1.2,  for  which  Dahlke  surveyed  the  flow  field.  Since  sur¬ 
veys  were  conducted  at 


X  =  47.  465  in. 


downstream  of  the  nozzle,  substitution  of  this  number  in  Equations  (121) 
and  (122)  yields: 


6lM  =0.  9  = 

oc 


(2.  582) 


it m 


‘Im  =1.2 

CO 


(2.  207) 


yff  -.<r- 

(z  )1,3*9 


Figure  55  shows  graphs  of  the  above  relations  as  calculated  from 
Dahlke's  data  in  Reference  46.  Although  there  is  considerable  scat¬ 
ter,  for  pressure  ratios  greater  than  10  the  points  do  seem  to  lie  on  a 
constant  line.  Based  upcn  Figure  55  for  F  >10,  the  average  values 


6  =  1.40  (123a) 


«  =  1.04  (123b) 

have  been  chosen. 

The  theoretical  model  for  the  jet  vortices  implies  that  certain  combin¬ 
ations  of  parameters  should  be  independent  of  the  pressure  ratio.  In 
particular,  the  product  of  normalized  vortex  strength  and  separation 
should  be  constant.  From  Equation  (1 10) 
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'AMTS  FOR  U  <1  (5)  AND  M  >1  (0 


U  sing  Equations  (1 1 1)  and  (1 1 2) 


(124) 


For  the  subsonic  case.  Equation (1 16) is  substituted  in  Equation  (124) 
above.  This  results  in  the  expression: 


(125) 


I 


I 
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Figure  56.  Scata  Constants  for  Equivatant  Subsonic  Jst  Oiametsre 


For  supersonic  free  streams,  substitution  of  Equation  (119)  into 
Equation  (124)  yields 


'M  >1 


Y  ~r  1 

j_  /y+A  Y‘1  j_ 

4ir  y  2  J  M, 


(126) 


Since  the  terms  in  brackets  in  Equations  ( 1 25)  and  (126)  are  constant 
for  fixed  Mac,  Dahlke's  data  should  indicate  that  the  term  (ry0/a«eP) 
is  independent  of  the  pressure  ratio  at  each  value  of  M,*.  This  seems 
to  be  the  case,  as  shown  in  Figure  56;  at  'east  to  within  approximately 
11%  of  the  averages  indicated  by  the  solid  and  dotted  lines  in  Fig¬ 
ure  56,  The  values  chosen  are: 


0. 0615 


(127a) 


M  <1 

oo 


IV) 

a"P/  MOC>1 


0. 0512 


(127b) 
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Figure  56.  Constant  Product  of  Vortsx  Strength  and  Sspsntion 
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These  values,  along  with  6  and  t  as  given  in  Equations  (1  23a)  and  ( 1 23b), 
are  substituted  into  Equations  (125)  and  (126),  respectively.  The 
appropriate  Mach  numbers  and  y  =  1.  4  are  also  substituted.  The 
results  are 


(K)m  0  =  0.155 

OO 

(128a) 

(k>m  >1  =  °-073 

(128b) 

OO 


Equations  (128)point  to  an  inconsistency  in  the  postulated  equivalence 
between  subsonic  and  underexpanded  jets.  Since  the  constant  K  is 
presumed  to  be  a  universal  constant  characteristic  of  the  subsonic  jet, 
the  values  of  K  calculated  starting  from  Mw  <  1  data  or  >  1  data 
should  coincide.'  Since  they  do  not,  it  appears  that  to  some  extent  the 
postulated  equivalences  are  not  valid.  Taken  individually,  however, 
both  the  Ma)<  1  and  >  1  data  indicate  that  K  is  a  constant,  and  using 
the  appropriate  value  of  K  for  each  case  leads  to  good  agreement 
between  predicted  and  measured  vortex  strengths,  as  is  shown  below. 
The  best  way  to  resolve  the  above  difficulty  would  be  to  evaluate  K  by 
using  experimentally  determined  vortex  strengths  for  a  subsonic  jet 
directly,  but  such  data  are  unfortunately  not  available. 

Referring  to  Equations  (106) and  (109),  the  model  predicts  that  the  follow¬ 
ing  ratio  should  be  constant. 


This  ratio  is  shown  in  Figure  57,  for  both  the  Moc<l  and  the  M^»l 
cases.  The  data  again  indicate  that  for  P>10,  the  ratio  is  approxi¬ 
mately  independent  of  P,  as  predicted.  The  average  values: 
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Figure  57.  Constant  Ratio  of  Vortex  Separation  to  Height  Above  the  Nozzle  Plane 


have  boon  chosen.  Using  Equations  (128a)  and  (128b)  these  lead  to 

(Y  )  =  0.  653  (12'M) 

1  v  M  <  1 

or 


(Y  )  =  0.385 

M  >1 

or 

Again,  some  inconsistency  is  evident  since  these  are  not  the  same. 

The  above  analysis  indicates  that  vortex  strength  varies  as  a  function 
of  the  similarity  variable,  Xi  which  is  dependent  upon  distance  along  the 
jet  axis,  pressure  ratio  P,  and  the  tree-stream  Mach  number  M, , 
Dahlke's  measurements  in  Reference  46  were  taken  at  a  fixed  location 
downstream  of  the  nozzle,  but  since  the  pressure  ratio  was  varied,  the 
similarity  variable  x  has  been  varied.  The  model  should  therefore  be 
able  to  predict  the  variation  of  T  with  pressure  ratio.  Combining 
Equations  (109),  (110),  and  (120),  the  following  relation  is  obtained: 


r* 


1 


3K  _ 

Y  3  (1.  76)  X  °* 28 

V 
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3uuBi.iLui.ing  for  X  anti  I  v  from  Equations  (111)  and  (I  IZ),  and  sub¬ 
stituting  the  appropriate  numerical  constants  it  is  possible  to  arrive 
at  the  following  expressions: 


=  0. 065  P0,  64 


=  0.114 


p0.64 


(inche  s) 


(inches) 


(130) 


(131) 


The  above  results  hold  for  x  =  47.  465  in.,  which  is  the  station  down¬ 
stream  of  the  orifice  at  which  flow  field  surveys  were  conducted  in 
Reference  46.  The  above  equations  are  compared  to  the  data  in  Fig¬ 
ure  58,  and  it  is  evident  that  the  agreement  is  quite  good.  The  formulas 
developed  above  will  be  used  for  predicting  vortex  strengths  and  posi¬ 
tions  at  the  aft  fin  location,  and  thus  for  calculating  jet-fin  interference 
forces  and  moments. 
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Figure  58.  Comparison  of  Predicted  and  Exparimantal  Vorttx  Strengths 
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5.  Z  JET-FIN  INTERFERENCE  FORCES  AND  MOMENTS 

The  interference  effects  calculated  are  the  incremental  forces  and 
moments  acting  on  the  fins  alone,  which  would  be  produced  by  turning 
on  a  control  jet  with  the  missile  at  given  angles  of  incidence  and  bank. 
The  general  approach  is  to  first  obtain  the  additional  upwash  which  the 
jet  vortices  produce  at  the  fin  location,  and  then  use  simple  two- 
dimensional  strip  theory  for  calculating  forces  and  moments  on  the 
fins. 

Reference  axes  and  angles  which  are  used  in  the  calculation  of  forces 
and  moments  are  shown  in  Figure  59.  The  missile  attitude  relative 
to  the  free  stream  coordinates  (x*,  y*.  z*)  is  defined  by  the  angles  of 
incidence, or ,  and  of  bank,  $  ,  as  in  Reference  49.  As  noted  in  the 
figure,  however,  the  missile  is  pitched  and  banked  about  axes  centered 
at  the  nozzle  station,  instead  of  about  the  nose.  Forces  and  moments 
are  defined  with  respect  to  the  body-oriented  coordinate  system  |x', 
y 1 ,  z').  It  will  be  assumed  throughout  that  only  one  control  jet  is 
turned  on;  that  the  jet  is  sonic;  and  that  the  nozzle  is  aligned  with  one 
of  four  cruciform  fins,  as  shown  in  Figure  59. 

5.3.1  Vortex  Strengths  and  Vortex  Locations  in  the  Body-Oriented 
Coordinate  System 

Let  t  denote  the  distance  along  the  body  axis  between  the  nozzle  and 
the  midpoint  of  the  fin's  geometric  mean  chord  (assumed  to  be  approx¬ 
imately  the  fin  center  of  pressure).  For  purposes  of  calculating 
interference  effects,  it  is  necessary  to  know  the  strength  and  location 
of  the  jet  vortices.  Consequently,  a  jet-oriented  coordinate  system  is 
introduced,  with  suitable  transformations  defined  below.  As  the 
missile  pitches  and  banks,  it  is  assumed  that  the  jet  remains  aligned 
with  the  free  stream ,  although  the  axis  of  the  nozzle  banks  with  the 
missile.  When  the  angle  of  incidence  is  different  from  zero,  the  jet 
will  no  longer  he  normal  to  the  free  stream  at  the  nozzle,  but  this 
effect  is  neglected  because  it  is  small  for  small  angles  of  attack. 

Coordinate  transformations  between  the  (x* ,  y*.  z,!t)  axes  and  the  body- 
oriented  (x1,  y1,  z')  axes  system  are  given  in  Reference  49.  They  are: 


x'  =  x1!'  cos®  -  z*  sina 


(132a) 


y'  =  -x*  (sina  sin<j>)  +  y#  cos$  -  z#  (cos a  sin$) 


(132b) 


z'  =  x#  (sina  cos<)i)  +  y*  oin$  +  z*  (cosa  cos<£) 


(132c) 
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New  reference  axes  are  now  introduced,  denoted  by  (x",  y",  z"). 


These  are  also  centered  along  the  body  centerline,  at  the  nozzle 


station.  The  x"  coordinate  is  aligned  with  the  free-stream  direction, 


and  one  plane  of  the  system  coincides  with  the  plane  defined  by  the 
free-stream  {or  x*)  direction  and  the  nozzle  centerline  (or  z')  direc¬ 
tion,  It  is  therefore  possible  to  define  the  unit  base  vectors  (i",  j", 
k")  for  this  new  coordinate  system  as  follows: 


i"  =  i* 


(133a) 


(k'  x  i#) 


l  -  lii'  x  j*i 


(133b) 


k"  =  j"  x  j"  =  i*  x 


k'  x_j# 


|k'  x_j#| 


g  .  j*  (i*  •  £') 


|k'  x  i*| 


(133c) 


The  j,  and  k  vectors  correspond  to  the  x,  y,  and  z  directions, 
respectively,  in  whatever  coordinate  system  is  indicated  by  the  super 
script.  Equation (13 3b) defines  the  i"  vector  as  being  normal  to  the 
plane  defined  by  the  free-stream  direction  and  the  nozzle  centerline, 
and  Equation  (133c) defines  the  lj"  vector  as  being  normal  to  the  i"  and 


I* 


,j"  vectors.  Using  the  transformation  in  Equations  (i 32),  it  is  possible 
to  derive  the  relations  between  the  (x*,  y*,  z#)  axis  base  vectors  and 
the  (x1,  y'  z')  axis  base  vectors.  These  are: 


=  i1  (cosa)  -  j1  (sina  sin$)  +  k1  (sina  cos«|>) 

j*  =  j'  (cos4>)  +  k'  (sin4>) 

Jj*  =  -Utsina)  -  j'  (cosa  sin$)  +  k1  (cosa  cosd>) 
Substitution  of  Equations  (134)  into  Equations  (133) then  gives: 
iT  =  (cosa)  -  j1  (sina  sin$)  +  It.1  (sina  cos<J>) 


r  -  r 
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sina  sin4> 

+£ 

cosa 
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.  2  2. 

-  sin  a  cos  9 

\fl  -  sin2o  cos2<|> 

»  ■ 

k"  =  -i' 

sina  cosa  cos^ 

+  3' 

2 

sin  a  sin<|>  cos4> 

1  2  2 
\l  1  -  flin  a  cos  4> 

L  .2  2 . 

yjl  -  sm  a  cos  9 

y  j 

L  ”  J 

+JslV 

.  .  2  2. 

1  -  sin  a  cos  9 

(134a) 


(134b) 


(134c) 


(135  a) 

(135b) 


(135c) 


With  the  base  vector  transformation  given  by  Equations  (135),  it  is  easy 
to  arrive  at  the  coordinate  transformations: 


x1  =  x"  [cosa]  +  y" 


sina  s inb 
LVl  -  sin^a  cos^fj 


-z 


sina  cosa  cos$ 


4 


.  .  2  2., 

1  -  sin  a  cos  4> 


(136a) 


115 


Finally,  the  jet-oriented  (x,  y,  z)  coordinate  system  is  defined.  It  is 
the  same  as  the  (x",  y",  z")  axis  system,  but  its  origin  is  located  at 
the  nozzle  exit  instead  of  at  the  axis  of  the  body.  Letting  R  denote  the 
radius  of  the  missile  at  the  nozzle  location,  the  transformation 
between  jefc -oriented  and  body-oriented  coordinates  iBj 


x'  =  x  [coso]  +  y 


sing  sin<j> 

r  71  27 

1  -  sin  a  cos  9 


sing  cosg  cos9 
J  1  -  sin^or  cos2<|> 


(137a) 


The  equations  for  vortex  strengths  and  positions  derived  in  Sub¬ 
section  5.  1  are  based  upon  the  (x,  y,  z)  coordinate  system.  For  pur¬ 
poses  of  computing  fin  interference  forces  and  moments,  it  is  desired 
to  know  the  location  and  strengths  of  the  jet  vortices  at  the  body  sta¬ 
tion  located  at  distance  I  downstream  cf  the  nozzle  station,  and  in  a 
plane  perpendicular  to  the  body  axis.  The  coordinates  of  the  fin 
station  in  rhe  body-fixed  system  are: 

x'  =  f 

y1  =  Q 
z1  ---  0 
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Letting  the  subscript  f  denote  the  coordinates  of  this  same  point  in 
the  jet-oriented  system.  Equations  (i  37) may  be  inverted  to  obtain  the 
result: 

x^  =  feosev  -  R  sina  cos$  (138a) 


Yf 


l  sina  sin<j> 


yjl  -  sin2a 


() 38b) 


cos*"4> 


=  ■< 


sinacosacosji  .  R  Vl  -  sin2Q  cos> 

r  n  t. 

1  -  sina  cos  9 


(138c) 


The  equation  for  a  plane  normal  to  the  missile  centerline  at  I  has  the 
following  vector  form 


(139) 


Equation (1 35)  mav  be  inverted  to  obtain  in  terms  of  the  base  vectors 
(i".  k"),  which  are  identical  to  the  jet-oriented  base  vectors  (I,  J, 

k),  with  the  result:  "*  ** 

i1  -  1  fcosal  +  f  sing  sin$  .  £  sina  cosa  cos<t> 

’  L  J  -  / - — J - 1  I - 2 - 2~ 

Jl  -  sina  cos  $  Jl  -  sina  cos  4 

When  this  is  substituted  into  Equation  (1 39),  the  expression  obtained  for 
the  plane  in  terms  of  jet-oriented  coordinates  is 


[cosa]  (x  - 


i-iJ+  f  _  I  (y.yj 

[Vl  -  sin2(J  cos2* J 


sina  cosa  cos$ 
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,  .  2  2. 

1  -  sin  a  cos  9 


lz  -  z£)  -■  0 


(140) 
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The  intersection  of  this  piano  with  'lie  jet  trajectory  must  now  be  found, 
in  order  to  determine  the  point  at  which  jet  properties  should  be  calcu¬ 
lated.  This  is  done  by  ncdng  that  the  jet  trajectory  may  be  defined  by 
the  relations: 


trajectory,  so  that  the  distance  at  which  the  jet  vortices  intersect  the 
plane  will  be  somewhat  different  from  that  given  by  Equation (144). 

This  effect,  however,  is  neglected. 

The  formulas  obtained  in  Subsection  5.  2  may  be  used  to  calculate  the 
location  and  strength  of  the  jet  vortices.  The  transformation  Equa¬ 
tions  (137)  will  then  give  their  location  in  the  hody-oriented  coordinate 
system,  at  a  station  located  a  distance  /  dov/nstream  of  the  nozzle, 
and  in  a  plane  normal  to  the  body  axis 

5.  3.  2  Jet  Fin  Interference  Forces  and  Momenta 

As  previously  mentioned,  it  is  assumed  that  crossflow  velocities  are 
much  smaller  than  the  free-stream  velocity.  Further,  it  is  assumed 
that  the  body  cross-section  does  not  change  at  the  fin  location,  (  ,  and 
that  the  vortex  strength  does  not  change  very  much  over  the  space  of 
a  body  radius.  Consequently,  the  flow  induced  by  the  jet  vortices 
about  the  missile  body  is  equivalent  to  the  incompressible  flow  induced 
about  an  infinite  circular  cylinder  by  two  infinite  counterrotating 
vortices,  as  illustrated  in  Figure  60.  The  strength  and  location  of 
the  vortices  are  taken  to  be  those  at  a  distance  t  downstream  of  the 
nozzle.  They  are  calculated  by  methods  described  in  Subsections  5.2 
and  5.  3.  1.  The  angles  of  attack  and  side  slip  induced  by  the  vortices 
at  the  fin  locations  are  calculated.  It  is  assumed  that  the  fins  are  flat 
plates  of  high  aspect  ratio,  so  that  simple  strip  theory  may  be  used  to 
calculate  the  induced  forces  and  moments,  as  described  in  Reference 
50,  Figure  61  shows  schematically  the  upwash  and  sidewash  induced 


by  a  positive  vortex  ot  strength  T,  at  the  plan<  s  z'  =  0  and  y*  =  0, 
respectively.  The  complete  vortex  system  contributing  to  interference 
effects  consists  of  the  two  jet  vortices,  Ti,  a;  d  rrc,  as  well  as  their 
images  inside  the  cylinder.  Since  all  vortices  have  the  same  strength, 
a  positive  strength  T  is  introduced  such  that 


where  the  subscript  (i)  denotes  the  image  vortex  inside  the  cylinder. 
The  vortex  system  is  depicted  in  Figure  62. 

Following  Appendix  B  of  Reference  50,  the  upwash  produced  by  these 
vortices  at  the  y'  axis  is: 


(145a) 


(145b) 
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(148a) 


1  "  2ir 


'2  =  27 


(148b) 


vi.  =  ra 


fz'  -  z'j) 


y'i  2  +  <*'  -  z\  >  2 
*i 


(148c) 
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v2.  =  ‘27 


(*'  -  *'2  ) 


y'z  +  («'  -  ) 

i  “i 


(148d) 


The  fin  geometry  is  depicted  in  Figure  63.  The  variation  of  chord  c 
with  distance  from  the  axis  may  be  written  in  the  general  form 


c  =  c 


c  -  c, 
r  t 

s  -  R 


<y'  -  R) 


(149) 


or,  normalized  by  the  cylinder  radius  R,  in  the  equivalent  from 

(Y1  -1)1  (150) 


C  =  C 


r  ['  -  (fn)  <*’  -  » 


where: 


C  =  normalized  root  chord  (c  / R) 
r  '  r 

\  =  taper  ratio  { / c r ) 

S  =  normalized  semi-span  (s /R) 
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Y1  =  normalized  coordinate  (y'  /R) 


To  compute  the  induced  angles  of  attack  and  sideslip,  the  induced 
velocities  given  by  Equations  (145)and  (149)  are  to  be  divided  by  the  com 
ponent  of  free-stream  velocity  along  the  x'  axis.  For  small  angles  of 
attack,  a,  this  component  is  approximately  equal  to  Uw,  The  total 
induced  angle  of  attack  along  the  y1  axis  is  therefore 


and  the  sideslip  angle  along  the  z1  axis  is 


[V1  +v2  +  vl,tv2.] 


« 

(152) 


P 
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Using  strip  theory,  then,  the  normal  force  on  the  horizontal  fins  iB 
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J  Cjj  (tl)c(n)dT) 
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-R 

f  Ct  (n)c(h)dri 


(153) 


the  side  force  on  the  vertical  fins  is 
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and  the  rolling  moment  on  the  fins  is 
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Pitching  and  yawing  moments  are  obtained  by  multiplying  Fg  and 
Fy  respectively,  by  f.  In  the  above  equations,  represents  the 
two-dimensional  lift  coefficient  for  a  flat  plate  at  an  angle  of  inci¬ 
dence.  It  should  be  noted  that  the  above  equations  assume  that  the 
fins  are  independent  of  each  other  (i.  e.  ,  fin-fin  interference  effects 
are  neglected).  Also,  any  forces  induced  by  the  vortices  on  the 
cylindrical  portion  of  the  body  are  not  included. 
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The  relation  betwc  -u  and  the  induced  angles  of  attack  or  sideslip 
depends  on  whether  the  free  stream  is  subsonic  or  supersonic.  In 
the  subsonic  case,  Reference  51  gives  the  following  relation: 


and  for  supersonic  free  streams 


Pi 


where 


P=  -M*| 


The  derivation  of  final  results  for  the  induced  normal  force  in  the 
subsonic  case  will  now  he  carried  out  in  detail.  With 


M  >  1 
00 


4 

p 


and  Equation (145),  and  Equation  {15 3} may  be  written  in  the  form: 
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where  C j,  is  the  normal  force  coefficient  defined  by 
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If  all  lengths  in  Equation  (158)  are  normalized  by  the  cylinder  radius,  R, 
and  the  symmetry  of  the  fins  is  utilized,  so  that: 


C(*l)  =  C  (-9) 
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Then  it  is  possible  to  write  Equation (158) in  the  following  form: 


where  C(’l)  is  given  by  Equation  (150)l 

The  integrals  in  Equation (I60)are  evaluated  by  standard  techniques, 
and  the  final  result  is: 
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where: 
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The  aide  force  and  rolling  moment  coefficients  are  defined: 
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Integrations  similar  to  those  carried  out  for  the  normal  force  lead  to 
the  following  results  for  the  subsonic  side  force  coefficient 
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where,  for  example, 
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(i.e.,  G  is  obtained  by  interchanging  Y'j  and  Z'  j  in  Equation  (162), 
as  might  have  been  expected  from  the  symmetry  of  the  situation). 

Finally,  carrying  out  the  integrations  indicated  in  Equation (1  55), 
leads  to  the  following  results  for  the  Bubsonic  rolling  moment 
coefficient. 
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where: 


Comparison  of  Equations  (1 56)  and (1  57)  indicates  that  the  interference 
coefficients  for  supersonic  free  streams  may  be  obtained  from  subsonic 
results  by  simply  multiplying  by  (2/ir). 


5.  4  RESULTS  OF  JET-FIN  INTERFERENCE  CALCULATIONS 

The  formulas  derived  in  Subsection  5.  3,  in  ronjunction  with  the  B 

methods  for  predicting  vortex  strengths  and  positions  described  in 
Subsection  5.2,  have  been  used  in  two  computer  programs  for  calcu¬ 
lating  jet-fin  interference  forces  and  moments  on  a  cruciform 
missile.  One  of  these  programs  is  valid  for  subsonic  free-stream  Mach 
numbers,  and  the  other  for  supersonic  free-stream  Mach  numbers.  As  * 

previously  mentioned,  the  interference  effects  calculated  are  the 
incremental  force  and  moment  coefficients  which  are  induced  on  the 
fins  alone  by  the  presence  of  the  jet  vortices. 

A  brief  description  of  the  programs  is  now  given,  with  reference  to 
the  flow  chart  and  tables  of  Figure  64. 

Given  the  geometric  and  aerodynamic  input  parameters  listed  in 

Figure  64,  tho  program  first  calculates  the  appropriate  equivalent 

subsonic  jet  scale,  using  Equation  (lib)  in  the  subsonic  case,  ani 

Equation  (1 19)  in  the  supersonic  case.  For  specific  angles  of  attach 

a  and  bank  4>  ,  the  program  then  calculates  the  value  of  X  at  which 

vortex  properties  are  to  be  computed  by  solving  Equation (144)  .« 

numerically. 

At  this  point,  a  check  is  made  to  ensure  that  the  resulting  X  corresponds 

to  a  value  of  X  which  is  greater  than  5.  This  is  done  to  ensure  that  the  ■ 

fins  lie  within  the  "vortex  region"  of  the  jet  as  defined  in  Reference  22,  » 

The  restriction  is  necessary  because  the  model  for  the  jet  vortices 

is  only  valid  in  this  region.  The  program  is  terminated  if  the  condition 

X  >  5  is  not  met. 

If  the  above  test  is  passed,  the  formulas  derived  in  Subsection  5.  2  ' 

are  used  for  computing  the  strength  of  the  vortices  and  their  position  I 

relative  to  the  jet-oriented  (x,  y,  T)  coordinate  system.  Using  the 

transformation  Equation  (137),  the  vortex  positions  relative  to  the  body 

fixed  (x1 ,  y1 ,  z  ')  coordinate  system  are  finally  calculated.  After 

suitable  normalization  of  the  vortex  coordinates,  the  formulas  derived  | 

in  Subsection  5.3.2  are  used  for  calculating  interference  forces 

and  moments.  At  this  stage,  it  is  again  necessary  to  discriminate 

between  sonic  and  supersonic  froe-stream  Mach  numbers,  as 

described  in  Subsection  5.  3.  2.  For  each  value  of  the  free  stream 

Mach  number,  ,  and  pressure  ratio,  P,  the  program  prints  all 

the  interference  coefficients  listed  in  Figure  64  as  functions  of  the 

bank  angle  <j>,  for  each  value  of  the  angle  of  incidence  #,  The  results  ' 

of  some  sample  computations  are  shown  in  Figures  65  through  67.  As  , 

indicated  in  the  figures,  the  free-stream  Mach  number  is  0.  8  and  the  ; 

angle  of  incidence  is  2° .  » 

i 

i 
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GEOMETRIC 

PARAMETERS 

AERODYNAMIC 

PARAMETERS 


GEOMETRIC  PARAMETERS 

JET  DIAMEtER  <*j 

DISTANCE  BETWEEN  NOZZLE 
CENTERLINE  AND  FIN  HALF, 

CHORD  LOCATION  f 

:  FIN  SEMISPAN  l 

BODY  CROSS-SECTIONAL 
RADIUS  R 

FIN  ROOT  CHORD  6r 

FIN  TIP  CHORD  e, 

AERODYNAMIC  PARAMETERS 

FREE-8TREAM  MACH  NO.  Mw 

ANGLE  OF  ATTACK  a 

JET-CHAMBER-TO- 
FREE-STREAM  PRESSURE 
RATIO  P 

INTERFERENCE  COEFFICIENTS 

NORMAL  FORCE  Cz 

SIDE  FOHCE  Cy 

ROLLING  MOMENT  C|. 

PITCHING  MOMENT  Cm 

YAWING  MOMENT  CN 


NOTE: 

JET  EXIT  MACH  NUMBER  IS 
ASSUMED  TO  BE  UNITY. 
GEOMETRIC  PARAMETERS  MUST 
HAVE  CONSISTENT  DIMENSIONS. 


Figure  64.  Flowchart  for  Fin  IntirfarancB  Fore*  and  Momtnt  Computations 


INTERFERENCE  ROLLING  MOMENT  COEFFICIENT. 


The  missile  ar.d  fin  geometries  correspond  roughly  to  the  configuration 
tested  in  Reference  27  and  depicted  in.  Figure  4.  The  specific  values 
used  for  the  geometric  pararheters  are:  ‘  • 


l  =  57.  405  in. 
s  a  6.  855  in. 

R  =  2.  75  in. 
cr  =  ct =  5.5  in. 

The  nozzle  is  assumed  to  be  sonic,  as  the  program  requires,  and  !■» 
have  a  diameter 

d.  -  0.  22  in. 

J 

Pitching  and  yawing  moment  coefficients  about  the  nozzle  station 
may  be  obtained  from  Figures  65  and  66  by  simply  multiplying 
and  by  the  normalized  moment  arm  (f/R). 

The  behavior  of  the  curves  in  Figures  65  through  67  may  be  explained 
as  f.illows.  The  equations  of  Subsection  5,  2  as  well  as  the  data  of 
Reference  46  indicate  that  the  vortex  spacing  yQ,  vortex  height 
above  the  exit  plane 2Q,  and  vortex  strength  T  increase  with 
increasing  pressure  ratio  P,  for  a  fixed  Ma  and  a  fixed  station 
downstream  of  the  nozzle.  Generally,  F  grows  with  P  at  a  faster 
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Figure  68.  Schematic  Diagram  of  Vortax  Locations  and  Intarfaranca  Forets  for  Various  Roil  Anglos 
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’•ate  than  either  yQ  or  I0  ,  In  a  sense;  then,  an  increase  in  P  will  j 

produce  two  counteracting  effects  on  the  fin  interference  forces  and  I 

moments.  As  the  vortices  move  farther  away  with  increasing  P,  j 

interference  effects  should  decrease;  on  the  other  hand,  an  increase 
in  P  also  tends  to  increase  interference  effects  by  increasing  T . 

Unless  the  vortices  are  quite  close  to  the  fin,  the  latter  effect 

predominates,  since  T  increases  at  a  faster  rate  with  P  than  either 

yQ  or  7.0  .  However,  forces  induced  an  fins  located  near  the  vortices 

are  more  sensitive  to  vortex  position,  and  an  increase  in  P  may  well 

decrease  interference  effects.  These  trends  are  reflected  In  Figures  j 

65  through  67.  Figure  65  indicates  an  increase  in  C-  with  pressure  \ 

ratio.  (It  is  recalled  that  the  nosile  centerline  is  aligned  with  the 

Z-axis. )  On  the  (-♦her  hand,  side  force  and  rolling  moment  j 

coefficients,  which  are  primarily  governed  by  the  upwash  induced  by 

the  vortices  on  fins  that  lie  near  the  vortices,  may  either  decrease  or  j 

increase  with  P,  depending  on  the  bank'  angle  +  .  The  trends  obtained 

when  $  is  varied  may  be  explained  with  the  aid  of  Figure  68.  Because 

of  symmetry,  side  force  and  rolling  moment  coefficients  must  be 

zero  at  $  =  0°  and  at  4  3  180*.  The  normal  force  coefficient  must 

be  greater  at  «ji  =>  180*  than  at  ♦  ■  0%  since  for  a  positive  a  the 

vortices  lie  closer  to  the  fin*' at  the  former  roll  angle.  From 

Figure  68b  it  is  evident  that  the  side  force  should  be  negative  at 

t»  90°,  and  the  side  force  may  further  be  expected  to  peak  near 
-  90° .  Figure  68  also  indicates  that  at  $  ■  90®  the  contribution  of  the 
horizontal  fins  to  the  rolling  moment, is  positive  (counterclockwise), 
while  that  of  the  vertical  fins  is  negative.  Consequently,  it  is  not 
surprising  that  under  certain  conditions  the  interference  rolling 
moment  coefficient  changes  sign  as  ^  is  varied.  This  effect  1b 
evident  in  Figure  67,  for  the  case  P  =  40.  As  previously  mentioned,  fin 
interference  forces  for  subsonic  Ma  have  been  measured  during  the 
AMICOM-CAL  tests  described  in  Reference  28.  A  sample  case  for  , 
the  configuration  tested  yielded  large  errors  in  the  magnitude  of  the 
theoretical  interference  forces,  although  the  measured  trends  were 
predicted  correctly.  The  configuration  tested  in  Reference  28, 
however,  had  a  very  low  fin  panel  aspect  ratio--one  half  of  that  shown 
in  Figure  4.  Consequently  the  simple  strip  theory  used  in  calculating 
interference  forces  would  not  be  expected  to  apply  to  this  configuration. 

Experimental  data  are  available  for  supersonic  free  stream  Mach 
numbers,  and  a  comparison  has  been  made  using  the  supersonic  fin 
interference  program.  The  geometry  chosen  corresponds  to  the 
configuration  tested  in  Reference  31  for  the  center  nozzle  location. 

The  geometric  parameters  have  the  following  values 

d.  ■  0.  1 1  in. 
i 

<  =5.  688  in. 

sal.  375  in. 
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INTERFERENCE  NORMAL  FORCE  COEFFICIENT.  I  *  INTERFERENCE  NORMAL  FORCE  COEFFICIENT.  C* 


R  *  0.  687  in. 


cr  a  1.  375  in. 
ct  a  1.  375  in. 

The  relevant  experimental  interference  normal  force  coefficient  has 
been  obtained  from  data  by  the  operation. 

Cz  m  (Cz)  -  (Cz) 

fins  on  fins  off 


Theoretical  and  experimental  results,  are  compared  in  Figure*  69  and 
70,  for  pressure  ratios  of  60  and  100,  respectively.  The  agreement 
appears  to  be  satisfactory. 


* 


* 
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Section  6 
CONCLUSIONS 


After  extensive  study  of  the  problem  of  JI  with  a  subsonic  mainstream, 
it  appears  evident  that  adequate  analytical  models  can  be  formulated 
from  potential  flow  theory;  however,  more  accurate  models  will  haye 
to  be  based  on  viscous  flow  analysis.  In  general,  the  adequacy  of  any 
analytical  methods  developed  will  remain  unknown  until  more  detailed 
experimental  data  are  available. 

Very  few  measurements  of  the  behavior  of  an  under  expanded  jet 
plume  exhausting  into  a  subsonic  cross  flow  have  been  conducted. 
Measurements  of  the  jet  vortex  strengths  have  so  far  beer,  limited 
to  the  results  described  in  References  45  and  46,  and  these  were 
only  obtained  at  one  station  downstream  of  the  nozzle.  Measure¬ 
ments  of  the  flowfield  in  the  viscous  wake-like  region  on  the  lee¬ 
ward  side  of  an  underexpanded  jet  have  not  been  made,  although  it 
appears  that  this  region  influences  the  interference  pressure  distri¬ 
bution  very  strongly.  Tests  of  an  underexpanded  jet  exhausting  from 
a  flat  plate  have  been  conducted  by  AMICOM  concurrently  with  this 
study.  Data  resulting  from  these  experiments  will  provide  a  basis  for 
evaluating  the  analytical  models  developed  in  this  study  and  others 
reported  in  the  VTOL- related  literature. 

Based  upon  data  concerning  subsonic  jets  in  subsonic  mainstreams 
and  the  limited  underexpanded  jet  data  available,  the  behavior  of  the 
JI  interference  pressure  distribution  appears  amenable  to  empirical 
description.  The  Fourier  series  empirical  fit  method  developed  in 
this  study  is  expected  to  provide  a  relatively  convenient  and  accurate 
empirical  description  of  the  interference  pressure  distribution.  It 
has  been  shown  in  this  study  that  a  characteristic  dimension  of  the 
underexpanded  jet  plume  will  scale  the  interference  pressure  distri¬ 
bution  as  the  jet  exit  dimension  does  for  subsonic  jets.  Consequently, 
it  is  expected  that  data  from  tests  involving  an  underexpanded  jet  will 
be  easily  fit  by  the  Fourier  series  method,  with  jet  exit  Mach  number 
replacing  the  velocity  ratio  (Uoo/Ujjas  a  parameter. 

Several  semi-ernpirical  models  of  the  interference  pressure  distri¬ 
bution  due  to  a  jet  in  a  subsonic  mainstream  have  been  developed  in 
this  study.  Generally,  the  approach  taken  has  been  to  postulate 
equivalent  flowfields  which  appear  plausible  either  on  physical  grounds, 
or  from  a  qualitative  knowledge  of  the  behavior  of  the  induced  pressure 
on  a  flat  plate.  These  flow  models  contain  empirical  constants  which 
are  determined  by  matching  to  experimental  pressure  distributions  in 
some  region.  The  validity  of  a  particular  model  is  then  judged  by  how 
well  the  pressure  distribution  is  represented  in  other  regions.  The 
model  flows  postulated  have  been  assumed  to  be  inviscid.  Reason¬ 
ably  good  representation  of  the  pressure  coefficient  distribution 
has  been  achieved  with  some  of  these  models.  Close  agreement 
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with  data  is  usually  restricted  to  a  particular  range  in  the  azimuthal 
angle.  Unfortunately,  the  models  which  give  best  overall  agreement 
with  data  are  those  which  3eem  most  unrealistic  on  physical  grounds. 
This  difficulty  is  tied  to  the  assumption  of  a  model  flow  which  is 
inviscid.  On  the  lee  side  of  the  jet,  a  realistic  model  must  include 
viscous  effects. 

The  equivalent  solid  obstacle  analogy  provides  a  basis  for  calculating 
approximate  JI  control  effectiveness  in  supersonic  flight.  In  general, 
the  behavior  of  amplification  factor  with  angle  of  attack  and  flight 
Mach  number,  as  well  as  with  jet  thrust,  can  be  predicted.  Account¬ 
ing  for  the  effects  of  boundary  layer  separation  by  redistributing  the 
pressures  due  to  the  inviscid  flow  about  the  equivalent  obstacle  does 
not  appear  to  increase  the  accuracy  of  the  prediction  method.  The 
effects  of  equivalent  obstacle  shock  reflection  from  the  vehicle  surface 
should  be  accounted  for  to  improve  the  accuracy  and  general  validity 
of  the  prediction  method.  Equivalent  obstacle  flowfield  analyses  which 
admit  nonaxisymmetric  shocks  are  required  in  order  to  acr'-.uit  for 
this  shock  interaction. 

Methods  for  calculating  jet-fin  interference  effects  for  subsonic 
and  supersonic  mainstreams  have  been  developed.  They  are  based 
upon  a  simple,  semi-empirical  model  of  the  jet-induced  vortices 
which  is  valid  at  large  distances  from  the  nozzle.  An  equivalence 
between  Bonic,  highly  underexpanded  jets  and  subsonic  jets  has  been 
postulated,  and  data  for  both  situations  are  used  in  evaluating  univer¬ 
sal  empirical  constants.  The  resultant  vortex  strength  variation  has 
been  shown  to  agree  quite  well  with  the  limited  data  available.  For 
given  vortex  strengths  and  locations,  simple  schemes  for  computing 
the  induced  load  on  control  fins  placed  well  aft  of  the  jet  nozzle  have 
been  developed.  One  of  these  is  valid  for  subsonic  free  stream  Mach 
numbers,  and  the  other  for  supersonic  freestream.  Results  of  sample 
calculations  show  that  for  a  fixed  Mach  number,  increasing  the  jet 
chamber  pressure  may  increase  or  decrease  fin  interference  effr.ots, 
depending  upon  the  relative  location  of  the  jet- induced  vortices  and 
the  fins  which  contribute  to  the  interference  force  or  moment.  Limited 
data  comparisons  for  supersonic  mainstreams  have  been  made,  and 
the  agreement  between  calculations  and  experiment  ha9  been  found  to 
be  satisfactory. 
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Appendix  A 

TERMINAL  SHOCK  LOCATION 
FOR  SUPERSONIC  NOZZLES 


In  this  appendix,  a  formula  is  derived  for  calculating  the. location  of  the 
terminal  shock  or  Mach  disk  in  a  highly  under  expanded  rocket  plume 
exhausting  into  still  air.  The  approach  is  the  analytical  equivalent  of  a 
graphical  method  developed  by  Adamson  and  Nicholls  in  Reference  A- 1. 
The  expression  for  the  terminal  shock  location  for  s  nozzle  with  exit 
Mach  number,  Mj,  greater  than  unity,  is  based  upon  a  simple,  semi- 
empirical  representation  of  the  Mach  number  distribution  along  the 
centerline  of  the  plume  of  a  sonic  nozzle. 

A.  1  SONIC  NOZZLE 

For  a  sonic  nozzle  exhausting  into  still  air,  it  has  been  shown  in 
Reference  A- 2  that  the  Mach  number  distribution  along  the  centerline 
of  the  plume  may  be  represented  by  the  flow  from  a  compressible  source 
whose  sonic  sphere  radius  is  given  by 

r*  =  (0.61)  d*  (A- 1 ) 

(Starred  quantities  will  refer  to  conditions  for  the  sonic  nozzle.  ) 

Flow  continuity  for  a  compressible  source  requires  that 

*  #  *  #  #2 
Pl  U1  r  =  P  u  z 


where  (  )j  denotes  conditions  at  the  sonic  radius  M*  =  1 

z*  is  a  coordinate  with  origin  at  the  nozzle  exit  plane  (Sco 
Figure  A- lb) 

The  above  relation  can  be  written  in  terms  of  stagnation  cpnditions  as 


(A- 2) 


Using  isentropic  flow  relations  and  Equation  (A-l),  Equation  (A- 2)  can 
be  written  in  the  form 


(A- 3) 
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b.  SONIC  NOZZLE 


Figure  A-1.  Mach  Disk  Haight  for  •  Suptnnnlc  Nozzla  Bawd  on  Sonic  Nozzla  Rasultt 
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Figure  A-2  compares  Equation  (A-3)  for  y  =  1.4  with  the  results  of  a 
method  of  characteristics  solution  reported  in  Reference  A-3.  Evidently 
the  agreement  is  excellent,  except  in  the  immediate  neighborhood  of  the 
orifice. 


With  the  subscript  f  denoting  flow  conditions  just  upstream  of  the  terminal 
shock,  Equation  (A-l)  and  Equation  (A-2)  can  be  combined  to  yield 


4  & 

The  final  pressure  and  temperature,  p#  and  Tj,  are  determined  by  the 
condition  that  the  pressure  across  the  Mach  disk  rise  to  the  external 
pressure,  pm.  Using  normal  shock  relations,  then,  the  above  equation 
may  be  written  as 


The  final  Mach  number,  M’jf,  may  also  be  determined  by  the  condition  that 
the  pressure  immediately  downstream  of  the  terminal  shock  be  equal  to  p^ 
and  normal  shock  relations  then  yield  the  equation 


V 


r  *2 

2YM*  -  (Y-l) 

9  • 

2 

*  *' 

Y+l 

9  9 

>ji2 

(Y-l)M*  +2 

(A- 5) 


It  is  not  possible  to  invert  Equation  (A-3)  to  obtain  M*  for  a  given 
pressure  ratio,  but  it  is  possible  to  obtain  this  quantity  from  Figure  A-3, 
which  is  a  plot  of  Equation  (A- 5)  for  Y  =  1.4,  Figure  A-3  Indicates 
that  for  moderately  high  pressure  ratios  Mf  is  large,  and  it  is  therefore 
possible  to  write  Equation  (A-4)  in  the  approximate  form 


(A-6) 


If  it  is  further  assumed  that  Ta  ■  Tq  ,  and  that  y  =  1.4,  Equation  (A-6) 
becomes  j 


h 
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Figure  A-2.  Mach  Numbtr  Dotrifeution  for  Sonic  NozzliPlunw 
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Figun  A-3. 


This  may  be  compared  with  the  empirical  formula  derived  by  Crist, 
Sherman,  and  Glass  in  Reference  A-4 


-S—  -  (0.645) 
d! 

J 


(A- SI 


Equation  (A-4)  predicts  that,  the  effects  of  varying  y  and  the  temper¬ 
ature  ratio  (T^/T^.)  will  be  small.  This  general  trend  has  been 
observed  experimentally  in  Reference  A-4. 

A.  2  SUPERSONIC  NOZZLE 


Following  the  method  of  Adamson  and  Nicholls  in  Reference  A-l,  the 
terminal  shock  location  for  a  supersonic  nozzle  will  now  be  derived 
based  on  the  sonic  nozzle  results  of  the  previous  section.  Referring 
to  Figure  A- la,  it  is  evident  that  the  flow  along  the  centerline  of  the 
supersonic  nozzle  plume  will  remain  undisturbed  until  the  first  expan¬ 
sion  fan  from  the  nozzle  lip  strikes  the  axis.  Let  the  distance  from 
the  exit  to  the  point  at  which  this  happens  be  denoted  by  Zf,  In 
Figure  A- lb,  the  distance  z&ij  denotes  the  point  at  which  the  expanding 
flow  in  the  subsonic  nozzle  plume  reaches  a  Mach  number  equal  fa 
the  exit  Mach  number  of  the  supersonic  nozzle.  Adamson  and  Nicholls 
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have  assumed  that  the  Mach  number  distribution  between  h  and  •/.(  is 
the  same  as  the  Mach  number  distribution  between  h*  and  zj£j.  . 
Assuming  once  again  that  the  flow  then  crosses  the  Mach  disk^to 
achieve  a  pressure  equal  to  p^,  it  follows  that 

h  =  z,  +h*- 

J 


Normalizing  by  the  exit  diameter  of  the  aupe r sonic  nozzle  yields 


It  will  now  be  required  that  the  sonic  and  supersonic  nozzles  have 
the  Bame  stagnation  conditions 

P*  -  P 
*o.  ro. 

J  .1 

T#  =  T 
o.  o. 

1  J 

and  then  isentropic  streamtube  ;  /ield 


Furthermore,  it  is  evident  from  Figure  A- la  that  if  the  initial  charac¬ 
teristic  is  assumed  to  be  straight,  the  ratio  (z//  d.  )  is  just  equal  to  one 
half  the  cotangent  of  the  Mach  angle.  ^ 

Thus, 


♦ 


* 


ZL 

d. 

J 
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(A- 1 2 ) 


A.  3  COMPARISON  WITH  DATA 


Figures  A- 4,  A-5,  and  A-6  compare  the  above  formulae  to  the  experi¬ 
mental  results  of  Love,  et  al.  (Reference  A-5).  For  M-  =  1,  some 
discrepancy  is  evident  in  Figure  A-4.  Figure  A-4  also  illustrates  the 
difference  between  taking  the  limit  of  Equation  (A- 13)  as  Mj  -» 1  and 
using  the  sonic  nozzle  results  of  Equation  (A-4).  In  both  instances, 
the  approximate  formulas  obtained  by  assuming  M|'  >>  1  give  results 
which  are  in  close  agreement  with  those  obtained  by  using  the  actual 
value  of  Mf  given  by  Figure  A-3.  The  discrepancies  for  Mj  =  1 
due  to,  the  fact  that  supersonic  source  flow  does  not  accurately  repre¬ 
sent  the  Mach  number  distribution  in  the  immediate  neighborhood  of 
the  nozzle  plane,  as  shown  in  Figure  A-2.  Since  this  region  of  dis¬ 
crepancy  is  removed  by  the  method  of  calculating  h  for  supersonic 
nozzles,  the  agreement  with  data  improves  at  higher  nozzle  Mach 
numbers.;;  Figure  compares  theory  and  experiment  for  Mj  =  2, 
and  thVagfeement  is  very  good.  The  more  complicated  Equation  ^A-13) 
gives  slightly  better  agreement  with  data  but  the  difference  with  the 
results  obtained  by  assuming  M?  »  i  does  not  appear  to  be  significant. 
Figure  A-6  cbmpiares  theory  and  experiment  for  M;  =  3.  It  should  be 
"noted  that  for  this  Mach  number,  the  terminal  shock  does  not  exist  for 
pressure  ratios  of  less  than  90.  The  points  shewn  in  the  figure  below 
this  value  correspond  to  the  distance  from  the  exit  to  the  first  inter¬ 
section  of.  a  diamond  shock  pattern.  The  unshaded  data  points  on  the 
figure  a ru, true  Mach  disk  locations,  and  the  theory  again  predicts  them 
cfulte  accurately. 
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Figure  A4.  Conpvim  of  Tlwofy  and  Data  for  M:  = 
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Appendix  B 

SINGULARITY  STRENGTHS  FOR  VORTEX  MODEL 
B.  1  DERIVATION  OF  FOURIER  COEFFICIENTS 

As  discussed  in  Section  3.2.2,  the  vortexmodel  leada^PithefoUowiiig 
expression  for  the  pressure  coefficients  ’ 

•\  ,  -k': i:  J: «■" 

.2  M  iX 


V  A2 


.i'I-  v, 


e+p]  -  H 

,.\W  I-?;1' 

i;'1 1  -".''i  •'  :',v 


■n-VN 

VA. 


-  ‘A3  <‘o  *  f»' 


*2+aJ-az  ’ 


4  A3  %  '  I  2 


(B-l) 

t  ■  * , . •  i. *  ,  ‘  i  ,  /  . 

'  '■  -V  !■•  I  ’  .11  I}.’  •...  . 

The  objective  is  to  obtain  the  first  three  terms  in  the  Fourier  series' 
representation  for  Equation  (B-l)  by  means  of  the  relations: 


V')4/  Cp(r.8)  d9 


(B-2) 
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(n  =  1,2) 


(The  symmetry  of  the  model  about  the  plane  0  =  C  insures  that  all  sine 
terms  in  the  series  will  be  identically  zero. )  Since  pressure  coeffi¬ 
cients  predicted  by  the  model  and  determined  from  experiment  are  to 
be  matched  at  r 1,  if  is  isirhpler  to  start  out  with  Equation  (B  - 1) 
written  for  .r  =  1.  v  Written  in  real  variables,  Equation  (B-l)  takes  the 
form; 


.v 


|2A4  C<”  28  2‘A  {['  -  A2*o]  c°‘  20  + 

+[Aia+r2)+(;0+5o)IA2-l)]co.8  + 


+  [»o  •  A.  do  *  to'  -  Aj)  (B-4» 

where 

F(6)  “  a  +  a,  co8  0  +  a,  cos  20  (B-5) 

O  1  M 

and:  , 

*0  =  aUoto)2H4+io 

«,  =  -2(1+  t0t0>  (to  +  t„>  (B-6b) 
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(B-6e) 


a2  =  2  lolo 


Then  the  integral  in  Equation  (B-2)  yields 


•'„<»  =  <a5  +  a2>-a*(c0  +  if  4s  + 


2*Aj(;  - 


£  2iA3( 


(B-7) 


where 


.  i*  'i 


I  *  r_- _ dL. 

o,  J  a_  +  a.cosO  +  a. 


2  o  o  “1 


co«20 


(B-8) 


and 


b  +  b,  cosfi  +  b,  cos  26 
O  1 _ c _ 


a_  +  a,  cos 0  4-  a,  cos  20 
o  1  c 


d0 


(B-9) 


and 


bo  “  toto  ’  A1  ^*»o  +  *  A2 


(B-9a) 


V  AI<l+W  +  <lo  +  V  <A2-1> 


(B-9b) 


b2  =  1  -  a2 


(B-9c) 
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Evaluating  Iq 


2 


first, 


it  is  advantageous  to  make  the  substituation: 


X  =  cos.  8 


♦ 


So  that  the  integral  becomes,  with  some  rearranging 


or,  factoring  the  quadratic  in  the  denominator 


+  1 


■-ib  / 


dx 


■1 


-  x2  (x  - 


(B- 10) 


c)  <x  -  c) 

where  it  may  be  shown  (using  Equation  B-6),  that 


—2 

t  +  1 

*o 


(B-lla) 


+  1 


(B-llb) 


The  integral  in  Equation  (B-10)  may  ve  evaluated  as  follows:  Consider 
a  complex  z  plane,  slit  along  the  real  axis  from  •  l  to  +1.  In  this 
plane,  consider  the  contour  integral 


l 


dz 


1  (z-c)  (z-c) 


where  C  is  the  contour  shown  in  Figure  B-l 


(B-  12) 
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Figuia  B-1.  Integration  Contour  for  Equation  (B-12) 


It  is  assumed  that  the  poles  c  and  c  do  not  lie  on  the  contour,  and  that 
the  ci-cle  of  radius  R  is  large  enough  to  enclose  both  singularities. 

In  that  case,  the  use  of  the  Residue  Theorem  a.3  discussed  in  Ref¬ 
erence  B-1  ^ielub 


9  = 


2wi 


(B-13) 


Writing  the  different  parts  of  the  integral  along  C,  and  then  allowing 


it  is  possible  to  show  that 


dx 


<x-c)  (x  -  c) 


(B-14) 


Then  Equations  ( B - 13)  and  (B-14)  lead  to  the  final  result 


rr 

2a j  (c_c) 


c 

The  integral 


(B-15) 


! 


4 


in  Equation  (B-9)  may  be  evaluated  in  an  analogous  fashion,  as  may  the 
other  integrals  which  arise  from  the  application  of  Equation  (B-3).  The 
final  equations  obtained  for  the  three  Fourier  coefficients  at  r  =  1  are: 


AT  + 


A2  ‘  a33 


a23  A2A3 


+  co(l)  =  0 


(B- 16) 


b33A3  +  b23A2A3  +  b12AlA2  "  bl3AlA3  “  ^iAl  "  cl^)  =  0  ^B"17) 


"33  3 


A,  +  c,,A,A,  -  c,,A,Aj  +Y2A-.  - 


'23  2  3 


'13  1 


2  '  *3A3 


A-  -  cJl) 


(B-18) 


where: 


a33  r 


«o  -  tb> 
UJo-1) 


(B-19a) 


a23  =  2i  "  to> 


(B- 19b) 


B.  2  SOLUTION  OF  THE  EQUATIONS 


Substituting  numerical  values  for  c0(l),  cj(l),  and  C2 <  1 )  evaluated 
using  Vogler'e  data  from  Reference  (B-2)  will  yield  the  five  simul¬ 
taneous  equations  for  the  five  unknowns  Aj,  A2.  A3,  r„.  and  0O.  These 
are  Equations  (31)  and  (32),  in  Section  3.2.2,  (B-16),  (B-i7),  and 
(B-18),  with  the  supplemental  Equations  (B-)9a)  through  (B-21e),  in 
which  it  is  recalled  that 


ie  -i®_ 

to  =  roe  °’  t0sro* 


It  is  possible  to  combine  the  above  set  of  equations  into  another  set 
which  is  more  suitable  for  numerical  calculation.  After  a  consider¬ 
able  account  of  algeibra,  the  results  are: 


Q4A*  +  QjAj  +  Q2A^  +  QjA3  +  Q0  =  0  (B-22) 


where 


Q4  =  b33^3  + 


4  (x?  -  1)  sin  20 

0  0 

C33P3  + 

4  cos0Q 

r 

0 

c‘3  (B-23a) 


Q3  =  ^33^3  + 


8(r  -  1 )  sin20^ 

'  O  O 


j  [c33  +ZP3r 


sin0 
o  o 


COS0 


+  3 ^ °  ^  +  ^c33  ®*n 


(B-23b) 
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Q?  =  4b33  -  [4  <r*-l)  ein2flo]  {P3C^{!)  -  8rosineo]  + 


4 


4 


8cos8 

+( — _ — a 


1  t  c33  +2P3  (ro5in0o)]  -0,(1)  p‘  - 


8cos8^  ,  , 

•I—? - ]  tc33c2<1)  "  8ro  ain  6o] 
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Qj  =  “  (8  (rJ+l)8in2eo]  c2<l)  -  J  ^3^(1)  -8r 
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o-  o 


-403^(1) 


(B-23d) 


4cos0  1 
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‘ 8cos0 
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r 

0 

0 

c2(l) 


(B-23e) 


an  expression  for  boj  is  given  by  Equation  (B»20a).  Written  in  real 
variables  this  has  the  form 
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3  2 

16r  sin  8  cos  0^ 
O _ 0 _ o 


33 


(rJ-D 


r^  -  2r^  coe28  +  1 
oo  o 


The  variable  P3  is  given  by 


(B-24) 


P3  =  -  (4r0sin8o  ) 


(rj"  -  1  )  <  1  +  2  cos  28  )  -  1 
0  o 


(B-25) 
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The  variable  C33  has  been  defined  in  Equation  (B-12a),  and  in  terms 
of  real  variables  it  has  the  form 


•> 

c,,  =  "4  sin  0 

JJ  C 
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<*» 


2 

o 


+  1)  l+2cos20  +  1 
_ _  o 

2r^  cos  20  +1 

o  o 


(B-26) 


Note  that  if  values  of  r_  and  0O  are  assumed,  Equation  (B-22)  becomes 
a  fourth  order  polynomial  for  A3.  The  other  equations  in  the  set  are: 


A  c2(1)  -  4A3(r08inV  ■  c33A3 
2  =  P3A3  +  2 


(B-27) 
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and 
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{4A2A3sineoK  +  {A1  +  A2  +  Co(1)}  rc 
■  j[4A2A38ineo]  [l+2co826oj  r  *- 


A1  +  A2  +  co^)J  [  1  +  2eos20  I  + 


+  4A2  sin2  6  1  r*  +  j  [4A,A,«in&J  |l  +  2  cos2ej|  r3  + 


axu  V 

J  O 


+  I  [  A1  +  A2  +  coU)|  (1  +  2cos26o)  - 


'»+|l 


2n3““’0 


-4A2sin^o  r2  -  4A,A,sine_|  r_ 


2  3  o  o 


A1  +  A2  +  c  (1)  =  0 


(B-30) 


Eqiiations  (B-22)  through  (B-30)  have  been  used  in  a  numerical  scheme 
to  obtain  solutions.  Starting  with  assumed  values  for  rQ  and  Ba,  the 
coefficients  of  the  fourth-order  polynomial  in  Equation  (B-22)  are 
calculated  using  Equations  (B-23a)  through  (B-26).  All  roots  of  this 
polynomial  are  then  obtained  numerically,  and  any  complex  roots  dis¬ 
carded.  The  remaining  values  of  A3  are  used  to  calculate  correspond- 
mg  values  A2  and  Al  by  means  of  Equations  (B-27)  and  (B-28).  The 
cubic  in  Equation  (B-29)  is  then  solved  numerically  to  obtain  new  values 
for0pj  and,  finally,  the  seventh -order  polymonial  in  Equation  (B-30) 
is  solved  for  rQ.  Any  complex  roots  are  again  discarded. 

Note  that  if  all  roots  in  the  above  scheme  are  real,  for  each  initial 
value  of  r0  and  60  there  exist  the  following  possibilities: 

4  possible  Aj,  A2>  A3's 

12  possible  0  's 
o 

84  possible  r  's 
o 
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In  practice,  however,  very  few  roots  turn  out  to  be  real.  In  addition 
to  requiring  that  all  variables  be  real,  two  more  constraints  are 
imposed  on  the  results.  It  is  known  physically  that  the  vortices  must 
lie  in  the  leeward  quadrants,  and  consequently  the  requirement 


0*  <  9  <  90* 

o 


(B-31) 


is  imposed.  Since  it  is  also  known  that  the  vortices  lie  near  the  jet, 
it  is  also  required  that 


0  <  rQ  <  1  (B-  32) 

B.  3  NUMERICAL  RESULTS 

The  validity  of  the  vortex  model  as  tested  by  carrying  through  a 
complete  case  for  Uw  /Uj  =  0.  4,  based  upon  the  data  of  Vogler  in 
Reference  B-2,  A  direct  iterative  scheme  as  suggested  by  the  discus¬ 
sion  above  was  unsuccessful  because  in  many  cases  an  input  (rD,  0O) 
pair  which  satisfied  the  constraints  of  Equations  (B-31)  and  (B-32),  did 
not  lead  to  any  calculated  (rQ,  0O)  values  which  satisfied  this  constraint. 
Therefore,  it  was  decided  to  search  the  matrix  to  acceptable  rQ  and  0O 
inputs  for  those  values  which  led  to  new  values  which  also  satisfied 
constraints  of  Equations  (B-31)  and  (B-32).  A  flow  chart  for  thiB 
scheme  is  shown  in  Figure  B-l.  Briefly,  the  program  takes  input 
initial  values  for  the  components  of  the  vortex  position  vector  (r^,  0O) 
which  satisfy  the  constraints  of  Equationo  (B-31)  and  (B-32),  and 
searches  for  solutions  to  the  set  of  five  nonlinear  equations  which  also 
satisfy  the  above  constraints.  The  most  relevant  input-output  diagram 
is  depicted  in  Figure  B-2,  The  diagram  indicates  how  the  acceptable 
output  varies  as  the  input  0O  is  varied  for  a  fixed  value  of  the  input  rQ. 
The  input  values  vary  over  the  entire  acceptable  range  (solid  lines), 
but  only  those  input  values  which  lie  within  the  regions  shown  by  dotted 
boundaries  lead  to  output  r„  and  0O  which  satisfy  the  constraints  of 
Equations  (B-31)  and  (B-32),  For  the  cases  shown  in  Figure  B-3,  the 
first  input  values  (solid  points) 


r  =  0.  65 
0 


8=5° 

o 
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lead  to  an  output  (open  point)  of 


77.  5* 


As  the  input  6  is  increased,  the  output  points  approach  the  input 
points,  until  at  input  point  rQ=  0,6 5,  6  Q  =  35*;  the  corresponding 
output  point  is  rQ  =  0.64,  0Q  =  36.  5*.  since  the  system  of  equations 
is  nonlinear,  more  than  one  output  may  correspond  to  a  single  input 
point.  This  fact  is  illustrated  in  Figure  B-2  by  the  points  labeled 
"second  branch".  As  the  input  &  is  increased  beyond  35*,  the  output 
values  no  longer  satisfy  the  constraints  (B-31)  and  (B-32),  »nd  are 
therefore  ignored.  At  rQ  =  0.65,  80  =  75*,  satisfactory  r  -put  values 
are  again  obtained,  and  they  are  shown  in  the  lower  right-hand  corner 
of  Figure  B-  1.  Exact  equality  of  input  and  output  would  signify  that  a 
solution  to  the  set  of  equations  has  been  found.  This  has  obviously  not 
yet  been  achieved  by  the  above  results,  and  iterations  in  the  neighborhood 
of  the  points  labeled  "approximate  covergence"  would  be  necessary  to 
obtain  a  more  accurate  answer.  However,  it  was  felt  that  the  above 
results  were  accurate  enough  to  allow  an  overall  evaluation  of  the 
model.  Examination  of  diagrams  analogous  to  Figure  B-2  has  shewn 
that  convergence  1b  not  approached  at  any  other  point.  The  final 
results  for  this  case  are  taken  as: 


A  j  =  -0. 1207 
A2  *  -0,0478 
A3  =  0.668 

*o  -  36.31° 


=  0.64 


Substitution  of  these  values  into  the  set  of  Equations  (B-22)  through 
(B-30)  verifies  that  they  are  approximately  satisfied.  These  are  the 
values  which  have  beon  used  for  plotting  the  model  three  term  Fourier 
series  representation  shown  in  Figures  20  to  26, 
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Appendix  C 

EQUIVALENT  SOLID  OBSTACLE  ANALOGY  COMPUTER  PROGRAM 


C.  1  MAIN  PROGRAM 
C.  1.  1  Input  Tape 

The  first  task  performed  by  the  main  program  is  that  of  reading  the 
input  tape.  The  input  tape  is  an  even-parity  BCD  tape,  density  800 
bits  per  inch,  which  contains  the  images  of  punched  output  from  a 
method  of  characteristics  analysis  of  the  flow  fields  around  a  unit 
hemisphere-cylinder  at  various  Mach  numbers  {Mf},  The  input  tape 
contains  the  images  of  the  following  cards: 

1.  The  first  card  contains  JMACH,  the  number  of  Mach  numbers 
in  the  sequence  {Mi}.  (For  operation  on  the  IBM  7094, 
JMACH-6. ) 

2.  The  next  JLIACH  cards  contain 
j,  nj,  Mj,  Xjmin,  Xjmax 

where  j  -  1,2,....  JMACH;  nj  is  the  number  of  net  points  in 
the  jth  flow  field;  Mj  is  the  Mach  number, 

Xj  .  and  Xj 
Jmin  Jmax 


denote  the  limits  of  the  stored  flow  field  data. 

(M>  Ximl„'  XW,) 

are  stored  in  the  arrays  AMLOC(6),  XMINC(6),  and 
XMAXC(6)  respectively. 

3.  The  next  cards  contain  the  actual  net  point  data 


j,  i,  Xji,  Rji,  Pji 


j  =  1,2,,...  JMACH 
i  =  l  *  2 . n  j 


where  (Xji.  R  jil  Pji)  are  logitudinal  and  radial 
coordinates,  and  the  associated  pressure  ratio,  p/pj. 
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4.  The  next  9  x  (JMACH)  cards  contain  the  coefficients  Cjm  of 
the  set  of  equations 

9 

X=  £  CjmRm-l;  j  =  1,  2 . JMACH 

m=l 


describing  the  shock  shape  associated  with  the  Mach  number 
Mj.  These  coefficients  are  stored  in  the  array  COEF(6,9). 

All  data  read  from  the  input  tape  is  stored  in  COMMON/BLK2/ 
so  it  is  available  to  the  integration  subroutine  (INTEG).  The 
input  tape  is  read  first,  and  none  of  the  information  in  BLK2 
is  altered  during  the  execution  of  the  program.  Thus  it  is 
necessary  to  read  this  tape  only  once,  regardless  of  the 
number  of  cases  to  be  run.  COMMON/BLK2/  accounts  for 
11,  773  storage  locations  in  the  program  when  JMACH  =  6. 

C.  1.  2  Input  Cards 

After  the  input  tape  has  been  read  and  COMMON/BLK2/  filled,  the 
input  cards  for  the  first  case  are  read.  The  data  on  the  cards  fill  the 
input  array  RR  (100),  which  contains  jet  and  free-stream  data 
described  in  Table  C-l  in  locations  1-19  and  vehicle  geometry 
specifications  in  locations  20-99.  Each  (K  =  1-100)  input  card 
contains  five  combinations  (K,  RR  (K))  in  a  format  specification 
(IX,  5  (13,  E9.4))  where  only  nonzero  values  RR  (K)  need  be  input. 

The  input  cards  for  each  case  must  be  preceded  by  a  card  containing 
the  number  of  cards,  NCARD,  to  be  read  for  that  particular  case 
punched  in  an  (110)  format. 

C.  1.  3  Program  Logic 

The  logic  involved  in  the  main  program  is  summarized  by  listing 
the  four  main  subroutines  in  the  following  manner: 

1.  SUBROUTINE  KWKBOD(NC,  RR) 

Inputs  via  calling  arguments:  Component  type  flags, 

Component  end  points 


Inputs  via  COMMON:  None 
Outputs  via  calling  arguments:  None 

Outputs  via  COMMON:  All  integration  mesh  data  trans¬ 
ferred  to  integration  subroutine  through  COMMON/MAIN/. 
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Table  C- 1 

INPUT  LOCATIONS  FOR  VEHICLE  GEOMETRY  SPECIFICATIONS 


Location 

Quantity 

Item 

1 

Free  stream  Mach  number 

2 

LBLSEP 

r>i  |0  No 

FU*  It  Yes 

3 

NCOMP 

Number  of  vehicle  components 

4 

db 

Vehicle  diameter 

5 

Lb 

Vehicle  length 

6 

6 

Nose  half  angle  (degrees) 

7 

xCg 

Vehicle  c*  g.  location 

8 

yeg 

Vehicle  c.  g.  location 

9 

zCg 

Vehicle  c.  g.  location 

10 

xj 

Jet  location 

11 

z  j 

Jet  location 

12 

CDIS 

Nozzle  discharge  coefficient 

13 

NPJ 

Jet  patch  number 

14 

dt 

Jet  diameter 

15 

4> 

Nozzle  cant  angle 

16 

Me 

Jet  exit  Mach  number 

17 

Jet  specific  heat  ratio 

18 

p°j/p„ 

NEW  TPM 

Jet  pressure  ratio 

19 

1(0  Newtonian  theory  on  lee  side 
[(1  Prandtl-Meyer  theory  on  lee  side 

(Locations  20  -  100  contain  the  body  section  inputs  as  described  in 
Section  C.  2. ) 


External  references:  Several  other  subroutines  are  called 
by  KWKBOD  during  vehicle  geometry  calculations.  Since 
these  require  no  special  inputs  or  handling  by  the  user,  a 
detailed  description  of  them  is  omitted  here. 

2.  SUBROUTINE  LOCFLO  (AMLOC,  JMACH,  ALPHA,  NPJ, 
GAMINF,  AMINF,  CONANG) 


Inputs  via  calling  arguments:  Number  of  local  Mach  number 
JMACH.  Sequence  of  local  Mach  numbers  [Mj]  .  Jet  patch 
number,  NPJ.  Nose  half  angle,  6  . 

Inputs  via  COMMON:  None 

Outputs  via  calling  arguments:  Sequence  of  angles  of  attack 


3. 


Outputs  via  COMMON:  Pressure  ratios 


[pi]£ 


(£i) 

VpJ  i 


at  jet  location. 


Local  flow  direction  vectors  [ux.]  [uyi ] '  [uZi] 
Angle  of  attack  flags  i_(JFLAG)i] 

Unit  outer  normal  at  jet  nj 


SUBROUTINE  JETHIT  (AMLOC,  GAMINF,  GAMJET, 
POJPIN,  DT,  JMACH,  THRUST,  PHI,  AME) 


Inputs  via  calling  arguments:  Jet  data 

Yj,  ,  dt,  ♦  ,  Me,  CDIS 
Sequence  of  Mach  numbers  [Mi] 

Inputs  via  COMMON:  Pressure  ratios  [Pj] 

Angle  of  attack  flags  [(JFLAG)^] 

Outputs  via  calling  arguments:  Jet  thrusts  [t^] 
and  normal,  sonic,  vacuum  thrust  [T3]  . 

Outputs  via  COMMON:  Sequence  of  equivalent  body 
radii  [Si]  . 

4,  SUBROUTINE  1NTEG  (JMACH,  CGX,  NPJ,  AMINF, 
GAMINF,  ALPHA,  AREF,  DREF,  TC) 


Inputs  via  calling  arguments:  Number  of  local  Mach 
Numbers,  JMACH. 

Free  stream  M,*,, 

Vehicle  c.  g.  location  xc,,,  yCg,  zcg 
Reference  area  and  length  Arcf,  dref 
Vehicle  angles  of  attack  [04] 

Jet  patch  number,  NPJ 
Jet  thrust  [Ti] 
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Inputs  via  COMMON:  Angle  of  attack  flags  [(JFLAG)ij 

Outputs  via  COMMON:  Vehicle  force  and  moment  coefficients 

External  references:  Another  subroutine,  PRESS,  is  called 
by  INTEG  during  the  integration  procedure. 

5.  SUBROUTINE  WRITR  (ALPHA,  DREF,  AREF,  AMINF,  N) 

Inputs  via  calling  arguments:  Vehicle  angles  of  attack  [ojJ 
Reference  length  and  area,  DREF,  AREF 
Free  stream  Moo 

Number  of  local  Mach  numbers  (N) 

Inputs  via  COMMON:  Jet-off  and  jet-on  aerodynamic 

coefficients 

Jet  thrusts  [Ti] 

Jet  pressure  ratio  (p0j/ Poo  ) 

Unit  outer  normal  at  jet  (nj) 

Outputs:  All  aerodynamic  coefficient  and  JI  amplification 
factor  outputs  are  printed  by  this  subroutine. 

C.  2  VEHICLE  GEOMETRY 

The  vehicle  is  assumed  to  be  fixed  in  the  coordinate  system  xyz  with 
its  nose  at  the  origin  and  with  the  positive  x-axis  as  the  axis  of 
symmetry.  The  cross  section  of  the  vehicle  can  be  made  up  of 
N  components,  where  N  lies  in  the  range  1SN5  8.  Thus,  the  input 
scheme  requires  that  each  component  be  described  by  a  separate 
curve  in  the  x-z  plane.  The  curves  may  be  one  of  three  types: 
straight  line,  circular  arc,  or  arbitrary  curve.  The  input  data 
for  the  first  curve  is  entered  in  locations  20-29  of  the  input  array, 

RR,  the  second  curve  in  locations  30-39,  etc,  A  summary  of  the 
inputs  required  for  each  type  of  curve  is  shown  in  Figure  C-l. 

It  should  be  noted  that  the  flag  denoting  the  type  of  curve  (1,  2,  or  3) 
must  be  entered  in  locations  20,  30,  40,  etc.,  and  the  coordinates 
entered  in  the  following  four  to  eight  locations. 

After  the  vehicle  shape  is  read  in,  the  integration  mesh  is  set  up. 

Each  component  is  divided  into  four  patches,  each  subtending  a  45° 
angle  on  the  surface.  As  is  also  indicated  in  Figure  4.  1,  the  free 
stream  velocity  vector  is  assumed  to  lie  in  the  x-z  plane,  so 
consideration  of  the  half- space  y  £  0  is  sufficient  for  vehicle  geometry 
considerations.  A  16  rectangle  per  patch  mesh  is  constructed  on 
all  the  patches  upstream  of  the  jet.  Aft  of  and  including  the  patch 
NPJ  on  which  the  jet  is  located,  the  mesh  fineness  is  chosen 
to  be  64  rectangles  per  patch.  The  patch  number  of  the  jet  is  an 
input  quantity  and  can  be  determined  readily  from  the  numbering 
system  indicated  in  Figure  4.  1.  If  an  even  finer  mesh  is  desired 
in  the  vicinity  of  the  jet,  the  reg.on  may  be  constructed  of  several 
short  components. 


X,  XJ 


a.  STRAIGHT-LINE  SEGMENT 


*2  *3  “4 


e.  ARBITRARY  CURVE 


INPUT 

LOCATION 

FLAG  -  1  .... 

.  10  ♦  ION 

*1 . 

.  11  ♦  ION 

*1 . 

.  12  ♦  ION 

*2 . 

.  13  +  ION 

*2 . 

.  14  +  ION 

FLAG  ■  2 . 10  ♦  ION 

*1 . 11  +  ION 

Hi . 12  ♦  ION 

*2 . 13  ♦  ION 

>2 . 14  ♦  ION 

$3  (DEGREES) . 16  ♦  ION 


FLAG  -  3  .... 

. 10  ♦ 

ION 

»1 . 

. 11  ♦ 

ION 

*1 . 

. 12  4 

ION 

“2 . 

. 13  ♦ 

ION 

‘2 . 

. 14  ♦ 

ION 

*3 . 

. 16  ♦ 

ION 

'3 . 

. 16  ♦ 

ION 

*4 . 

. 17  ♦ 

ION 

u . 

. 18  ♦ 

ION 

Figure  C-1.  Vehicle  Profile  Oeecription  Options 
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The  coordinate  system  described  above  is  a  standard  aerodynamic 
coordinate  system  for  the  user's  benefit.  Inside  the  program, 
however,  a  transformation  is  made  to  another  coordinate  system 
in  which: 

X  (INPUT/OUTPUT)  =  X(3)  or  operating  z  axis 
Y  (INPUT/OUTPUT)  =  -X(2)  or  operating  negative  y  axis 
Z  (INPUT/OUTPUT)  =  X ( 1 )  or  operating  x  axis 

for  the  purpose  of  internal  program  operation.  The  user  Bhould 
bear  this  in  mind  when  consulting  the  flow  charts. 

C.  3  LOCAL  FLOW  SUBROUTINE  (LOCFLO) 

In  this  subroutine,  where  the  angles  of  attack  Ja.l  corresponding  to 
local  Mach  numbers  {M^}  at  the  jet  location  are  Calculated,  several 
algebraic  equations  are  solved  numerically.  The  schemes  used 
(binary  chop  and  regula  falsi)  are  straightforward  and  should  present 
no  problem  to  the  user.  The  convergence  criteria  and  basic  logic  are 
readily  apparent  from  the  flow  charts,  so  no  further  discussion  is 
necessary  here. 

Although  the  program  has  the  capability  of  using  either  Newtonian  or 
Prandtl- Meyer  lee  side  aerodynamic  theories,  the  Prandtl-Meyer 
»  theory  is  used  for  determining  the  angles  of  attack  {cr^}. 

For  every  angle  of  attack,  o.  ,  corresponding  to  a  local  Mach  number, 
M^,  at  the  jet  location,  the  flag  (JFLAG)^  is  set  equal  to  zero.  If  no 
afc  is  found  corresponding  to  a  particular  Mach  number,  (JFLAG)k 
is  set  equal  to  1  indicating  to  the  integration  subroutine  that  no 
pressure  Integration  is  to  be  carried  out  at  this  angle  of  attack. 

If  the  input  nose  half-angle  is  less  than  85°,  the  attached  shock  version 
of  LOCFLO  is  executed.  If  the  angle  of  attack  corresponding  to  a 
particular  local  Mach  number  is  greater  than  0.  3  radian,  an  appropri¬ 
ate  message  is  printed.  If  data  is  desired  for  larger  angles  of  attack 
for  a  vehicle  with  a  sharp  nose,  it  is  recommended  that  the  program 
be  rerun  for  a  slightly  blunted  vehicle  of  the  same  overall  dimensions. 
In  this  case  the  input  nose  half-angle  will  be  90*  and  the  detached 
shock  option  will  be  executed. 

C.  4  PRESSURE  INTEGRATION  SUBROUTINE 

The  pressure  integiation  over  the  vehicle  surface  is  carried  out  by 
,  subroutine  INTEG.  This  subroutine  uses  the  integration  me«h 

constructed  by  subroutine  KWKBOD.  The  integration  is  performed  for 
the  vehicle  at  those  angles  of  attack,  a^,  of  the  sequence,  [a^]  ,  for 
which  the  flag,  (JFLAG)^,  of  the  sequence,  [  { JFLAC»)jJ  ,  was  set 
equal  to  zero  in  subroutine  LOCFLO. 

At  each  mesh  element,  before  calculating  the  surface  pressure, 
another  subroutine,  PRESS,  is  called.  This  subroutine  tests  to  see  if 
the  mesh  element  lies  within  the  region  of  influence  of  the  jet.  If  it 
does,  the  pressure  used  in  the  integration  scheme  is  calculated  by  the 
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method  described  in  section  C,  4.  1.  If  the  mesh  element  lies  outside 
the  region  of  influence  of  the  jet,  control  is  returned  to  subroutine 
INTEG  where  the  pressure  is  computed  without  regard  for  the  jet. 

In  subroutine  INTEG,  the  pressure  on  each  mesh  element  outside  the 
region  of  influence  of  the  jet  is  determined  by  a  local  inclination 
pressure  law.  On  the  winuward  side,  tangcnt-cone  theory  is  used.  On 
the  lee  side,  depending  on  the  value  of  the  input  flag  NEW  TPM,  either 
Newtonian  or  Prandtl-Meyer  hypersonic  small  disturbance  theory  is 
used. 


SUBROUTINE  PRESS  (X,  Y,  Z,  J.  KF LAG,  PEPI,  KPRESS)  is  given 
below. 

Inputs  and  Outputs 

Inputs  through  calling  arguments: 

X,  Y,  Z  are  the  coordinates  of  the  mesh  element  in  question. 
J  is  the  index  denoting  an  angle  of  attack  of  the  sequence 
{oi}. 

KPRESS(6)  is  an  array  zeroed  in  INTEG  and  used  tn  count 
the  times  PRESS  is  called. 

Outputs  through  calling  arguments: 

KFLAG  is  set  equal  to  one  if  the  point  (X,  Y,  Z)  is  inside 
the  region  of  influence  of  the  jet,  and  set  equal  to  zero 
otherwise. 

Inputs  through  COMMON/BLK1/: 

PIPINF  (6)  contains  the  pressure  [  P^  ]  at  the  jet  location. 

XYZJ(3)  contains  the  coordinates  (X.,  Y.,  Z.)  of  tho  jet 
location.  J  } 

UNJ(3)  contains  the  components  (nx,  n  ,  nz)  of  the  unit  outer 
normal  n  at  the  jet  location.  y 

UX(6),  UY(6),  UZ(6)  contain  the  components  of  the  series 

of  unit  vectors  [  ]  aligned  with  the  local  flow  at  the  jet 

location. 

SCALE(6)  contains  the  scale  factors  (  ]  which  are 

calculated  as  described  in  Section  4. 

NEWTPM  is  the  input  lee  side  aerodynamic  theory  flag. 

(0  -  Newtonian,  1  -  Prandtl-Meyer). 

IBLSEP  is  the  input  boundary  layer  separation  flag, 

(0  -  Inviscid  pressure  profile,  1  -  Modified  pressure  profile). 
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Inputs  through  COMMON/BLK2/: 

NCHAR{6)  contains  the  numbers  [NC:j  of  characteristic  net 
points  in  the  equivalent  body  flow  field  at  local  Mach 
numbers  [Mj]  . 

XCHAR(6,  650),  RCHAR(6,  650),  PCHAR(6,  650)  contain  the 
coordinates  [  XQ_  ]  ,  [  Rc^]  •  °*  the  net  points  and  their 

associated  pressure  ratios  [  Pc_  1  • 

XMINC(IO)  contains  the  abscissas  [X-  .1  of  the  most 

1  cmin  4 


forward  net  point. 

XMAXC(IO)  contains  the  abscissas  [Xc 
aft  net  point. 


]  of  the  farthest 


Beyond  this  point  the  pressure  perturbation  due  to  the 
hemisphere- cylinder  is  negligible.  The  pressure  perturbation 
is  ignored  at  radial  angles  in  body  cross-section  larger  than 
±150  degrees  from  the  jet  centerline. 

COEF(10,  9)  contains  the  coefficients  Cilt  in  the  equations 


shock 


9 

I 

k  =  1 


shock 


j  =  1,  2,  . .  .  JMACH 


for  the  shock  shapes  in  the  equivalent  body  coordinate  system, 

In  this  subroutine,  the  coordinates  of  a  point  (x,  y,  a)  are  transformed 
into  the  equivalent  body  coordinate  system  (X1  ,  R1  )  and  examined 
to  determine  whether  they  lie  within  the  region  of  influence  of  the  jet. 
The  limits  of  the  interaction  region  depend  on  whether  the  boundary 
layer  separation  simulation  option  is  being  exercised.  In  the  inviscid 
case,  the  upstream  limit  is  the  equivalent  body  shock  wave.  In  the 
case  of  the  modified  pressure  profile,  the  problem  is  reduced  to 
determining  whether  a  hypothetical  point  (X*,  R'  )  has  an  Inviscid 
pressure  ratio  associated  with  it  greater  than  the  experimentally 
observed  plateau  pressure. 

The  logic  involved  in  this  subroutine  may  be  considered  as  consisting 
of  three  main  parts.  First,  the  limits  of  the  interaction  region  are 
calculated;  second,  it  is  determined  whether  the  point  in  question  lies 
within  the  interaction  region;  and  finally,  the  appropriate  pressure 
ratio  is  assigned  to  the  point  (X1 ,  R' ). 

An  array  of  counters  KPRESS(6)  zeroed  in  subroutine  INTEG 
indicates  when  PRESS  is  being  called  for  the  first  time.  The  first 
time  it  is  called  for  a  given  angle  of  attack  (i.  o.  ,  a  particular  local 
Mach  number  and  equivalent  body  flow  field)  the  limits  R'min  and 
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R'max  corresponding  to  the  shock  end  points  X'min  and  X  max  read 
in  on  the  input  tape  must  be  determined.  The  numerical  methods 
UBed  to  solve  the  equation 


X'  =  > 

n  =  1 


(n-  1) 


are  unique.  Because  of  the  unknow  i  nature  of  the  above  curve  fit 
outside  the  interval  (X'min,  X'max),  the  iteration  schemes  were 
designed  to  conve-ge  *o  the  end  points  of  the  interval  (R’min.  R'max  )• 
from  the  interior.  To  find  R'max  »  5111  initial  R'  =  lb  is  guessed,  and 
the  iteration  proceeds  outward  in  steps  of  size 


/m2  -  1 


DR  = 


The  resulting  sequence  { R*  +  n  D  R}  approaches  R'max  monotonically 
because  the  Mach  angle  H  *  tan*l  -  l ) - 1  / 2  aiways  smaller 

than  the  shock  angle  at  any  point  X'8  (R' ).  To  find  R'min  »  an  initial 
value  of  R'  a  10  is  guessed  and  the  iteration  proceeds  toward  the 
X'  axis  in  steps  of  size 


<R,>-Xmm 


The  resulting  sequence  {  R1  -  n  DR  }  approaches  R'min 

dX' 


monotonically 


- 1  aA  s 

because  the  shock  angle  P  =  tan  — is  a  monotonic  increasing 

function  of  R'  and  always  greater  than  zero.  The  pattern  of  Iteration 
in  both  the  above  cases  should  be  clear  from  Figure  C-l.  An  examina¬ 
tion  of  the  flow  charts  for  subroutine  PRESS  will  show  that  there  are 
safeguards  In  the  event  of  a  numerical  overshoot  of  an  end  point. 

The  second  task  of  the  subroutine  is  to  determine  whether  a  point  lies 
inside  the  interaction  region.  Clearly,  ifX1  5  X*  max  or  if 
R1  >  R'max.  the  point  in  question  1b  outside  the  region  and  control 
la  returned  to  INTEG  with  KFLAG  set  equal  to  zero,  indicating  that 
the  surface  pressure  at  the  point  (x,  y,  z)  is  unaffected  by  the  presence 
of  the  jet.  Next,  the  shock  abscissas  are  found  for  the  particular 
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value  of  R’  in  question.  If  R' 
station  is  given  by 
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it  is  clear  that  the  nhock 


X’  - 

s 


C  R  (n  *  1} 


n  =  1 
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If  H  v  Rmm  »  ^  normal  shock  Is  ossutncd  (since  the  point  (X  ■  R  } 
is  near  the  sonic  line  on  the  equivalent  body)  and  X's  is  set  equal  to 
X'^p.  The  inviscid  pressure  ratio  Pmv  across  the  shock  i^  then 
calculated  for  the  appropriate  shock  angle  p. 


R  min 


P  =  tan 


■*u  * 


min  S  R’  * 


If  tile  boundary  layer  separation  simulation  is  being  exercised 
r.BLSEP  *  1),  a  new  abscissa 


=  <X.  *  Xi  >  1X70 c;>  +  Xs 


is  computed  for  every 


point  upstream  of  .  The  original  value  X^  is  stored  in  XTEMP 

and  the  flag  LFLAG  is  set  equal  to  one  to  indicate  to  later  sections  of 
the  subroutine  that  the  point  (X1 ,  R1  )  is  an  image  of  an  upstream 
point.  In  the  inviscid  case,  KFLAG  is  sriroed  and  control  returned  to 
INTEG  if  X'<  X*  8. 

The  pressure  perturbation  ratio  lJ  associated  with  a  point  (X1  ,  R'  ) 
is  found  by  locating  the  nearest  characteristics  net  point 
(XCHAR  (J,  Ic),  RCHAR  (J,Ic))  in  the  stored  equivalent  body  flow 
field  and  using  its  associated  pressure  ratio  PCHAR  (J,  Ic).  The 
index  J  indicates  the  local  Mach  number,  and  Ic  the  I-4*1  net  point, 
when  arranged  in  order  of  ascending  X1  coordinates.  Instead  of 
searching  all  NCHAR(J)  points  for  the  nearest  one,  searching  is 
confined  to  a  circle  of  radius 


d  »  min 


R  -  R'  min 
X’  -  X’s 


around  the  point  (X'  ,  R'  ).  This  is  accomplished  by  determining  the 
index  I3  corresponding  to  the  net  point  with  the  largest  XCHAR(J,  I) 

<  X'  -  d  and  the  index  I*  corresponding  to  the  point  with  the  smallest 
XCHAR(J,  I)  >X'  +d.  Then,  the  search  for  the  nearest  point  need  take  place 


REFERENCES 


1.  Vinson,  P.  W.  ,  Amick,  ,1.  L.  ,  anil  Liepmsnn,  H.  P.  Interaction 
Effects  Producted  by  Jet  Exhausting  Laterally  Near  Base  of 
Ogive- Cylinder  Model  in  Supersonic  Mainstream.  NASA 

Memo  12-5-58W,  February  1959. 

2.  Amick,  J.  L.  and  Hays,  P.  B.  Interaction  Effects  of  Side  Jets 
Issuing  from  Flat  Plates  and  Cylinders  Aligned  with  a  Super¬ 
sonic  Stream.  WADI)  TR-60-329,  June  I960. 

3.  Phinney,  R.  E.  ,  Werle,  M.  J.  ,  Knott,  J.  ,  and  Volz,  W.  C. 

Slot  Jet  Interaction  Studies  of  an  Ogive- Cylinder  at  M^  =  4 
and  5.  NOLTR-68- 143 ,  September  1968. 

4.  Spring,  D.  J.,  Street,  T.  A.,  and  Amick,  J.  L.  Transverse 
Jet  Experiments  and  Theories-- A  Survey  of  the  Literature, 

Part  I.  U.  S.  Army  Missile  Command,  Redstone  Arsenal, 

Report  RD-TR-67-4,  June  1967. 

5.  Barnes,  J.  W.  ,  Davis,  J.  G.  ,  and  Tang,  H.  H.  Control 
Effectiveness  of  Transverse  Jets  Interacting  With  a  High-Speed 
Free  Stream.  AFFDL-TR-67-90,  Vol.  I,  July  1967. 

6.  Carvalho,  G.  F'.  and  Hays,  P.  B.  Jet  Interference  Experiments 
Employing  Body-Alone  and  Body-Fin  Configurations  at  Supersonic 
Speeds.  University  of  Michigan  Technical  Report  03942- 7- T , 
GM-979.  December  1970. 

7.  Swanson,  R.  S.  Parametric  and  Preliminary  Design  Comparison 
of  Transverse  Jet  Controls  and  Flap- Type  Aerodynamic  Controls 
for  Re-entry  Vehicles  and  High  Speed  Cruise  Aircraft. 
AFFDL-TR-68-98,  December  1968. 

8.  Spring,  D.  .1.  and  Street,  T.  A.  Experimental  Investigation  of 
the  Pressure  Distributions  Induced  Upon  a  Body  of  Revolution  by 
Several  Transverse  Jets  at  Low  Speeds.  U.  S.  Army  Missile 
Command  Report  No.  RD-TM-68-7,  August  1968. 

9.  Kuiper,  R.  A.  Control  Jet  Effectiveness  in  the  Subsonic  and 
Transonic  Flight  Regimes.  Philco  Aeronutronic  Division 
Publication  No.  U-  2932,  December  1964. 


183 


10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


j 

i 

i 

-5 

} 

j 

Heyser,  A.  and  Maurer,  F.  Experimental  Investigations  on  ! 

Solid  Spoilers  and  Jet  Spoilers  at  Mach  Numbers  of  0.6  to  2.8,  ; 

Translation  No.  J2,  Jet  Propulsion  Laboratory,  21  February  1964.  j 

*  i 

Vinson,  P.  R.  Prediction  of  Reaction  Control  Effectiveness  at  j 

Supersonic  and  Hypersonic  Speeds.  Martin- Orlando  Report  j 

OR- 6487 ,  March  1965.  1 

Snaiu,  F.  W .  and  ZcVr»«ki,  E.  E.  A  Study  of  the  Interaction  of  ' 

Gaseous  Jets  from  Tran.3Vf.rse  Slots  with  Supersonic  External 
Flow.  AIAA  Journal,  Vol.  6,  No.  2,  February  1968. 


Werle,  M.  J.  A  Critical  Review  of  Analytical  Methods  for  ij 

Estimating  Control  Forces  Produced  by  Secondary  Injection.  5 

NOLTR  68-5,  January  1968. 

Kaufman,  II,  L.  G.  Hypersonic  Flows  Past  Transverse  Jets. 

Presented  at  the  AIAA  Fifth  Aerospace  Sciences  Meeting, 

Paper  No.  6-7-190,  New  York,  New  York,  23  to  26  January  1967. 

Strike,  W.  T.  Analysis  of  the  Aerodynamic  Disturbances 
Generated  on  a  Flat  Plate  Containing  Lateral  Jet  Nozzles 
Located  in  a  Hypersonic  Stream.  AEDC-TR-67- 1 58, 

January  1968. 

Spaid,  F.  W.  ,  Zukoski,  E.  E.  ,  and  Rosen,  R.  A  Study  of  I  j) 

Secondary  Injection  of  Gases  Into  a  Supersonic  Flow.  j. 

Jet  Propulsion  Laboratory  TR  32-834,  1  August  I960.  j 

Street,  T.  A.  and  Spring,  D.  J.  Experimental  Investigation  j 

of  a  Transverse  Jet  Injecting  From  a  Flat  Plate  info  a  Mach  | 

Number  5.0  Free  Stream.  U.  S.  Army  Missile  Command  j 

Report  No.  RD-TM-68-6,  August  1968.  Y 


Lexko,  W.  Loads  Induced  on  a  Flat  Plate  at  a  Mach  Number 
of  4.5  with  a  Sonic  or  Supersonic  Jet  Exhausting  N«>i  ml  to 
the  Surface.  NASA  TN-D-  193  5,  July  1963. 

Cassel,  L.  A.,  Davis,  J.  G.  ,  and  Engh,  D.  P.  Lateral  Jet 
Control  Effectiveness  Prediction  for  Axisymmetric  Missile 
Configurations.  U.  S.  Army  Missile  Command  Report 
No.  RD-TR-68-  5,  June  1968. 

Lee,  Jr.  ,  E.  E.  and  Willis,  C.  M.  Interaction  Effects  of  a 
Control  Jet  Exhausting  Radially  from  the  Nose  of  an  Ogive- 
Cylinder  Body  at  1  ransonic  Speeds.  NASA  TN-D-3752, 
January  1967. 


l 


i 


I 


184 


21.  Jorclinson,  R.  Flow  in  a  Jet  Direct  Normal  to  the  Wind 
t,i.  Brit.  Aero.  Res.  Council.  R&M  No.  3074,  1958. 

22.  Pratte,  B.  D.  ,  and  Baines,  W.  D.  Profiles  of  the  Round 
Turbulent  Jet  in  a  Crossflow.  Jour,  of  the  Hydraulic  Div.  , 

Proc.  of  the  ASCE,  November  1  967. 

2  5.  Knffer,  J.  F.  and  Baines,  W.  D.  The  Round  Turbulent 
Jel  in  a  Crosswind.  Jour,  of  Fluid  Mcch.  ,  Vol.  15, 
pp.  481-496.  1963. 

24.  Abranovich,  G.  N.  Theory  of  Turbulent  Jets.  Reproduced 
by  the  Armed  Services  Technical  Information  Agency, 

Arlington,  Va. 

25.  Reichenau,  D.  E.  A.  Interference  Effects  Produced  by  a 
Cold  Jet  Issuing  Normal  to,  the  Airstream  from  a  Flat  Plate 
at  Transonic  Mach  Numbers.,  AEDC- TR-67-220,  Arnold 
Air  Force  Sta.  ,  Tenn.  ,  October  1967. 

26.  Reichenau,  D.  E.  A.  Interference  Effects  of  Cold  and  Hot 
Rocket  Exhaust  Issuing  Normal  to  the  Airstream  from  a  Flat 
Plate  at  Free  Stream  Mach  Numbers  from  0.  6  to  1.4. 

AEDC- TR- 66- 127,  Arnbld  Air  Force  Sta.  ,  Tenn.  ,  June  1966. 

27.  Reid,  C.  F. ,  Jr.  The  Effect  of  Several  Forward-Mounted 
Control  Jet  Nozzles  on  a  Typical  Missile  Configuration  at 
Transonic  Speeds.  Cornell  Aero.  Lab.  ,  Inc.  Report 

No.  AA-2234-W- 1  (Vols.  1  and  II).  ,  Buffalo,  N.  Y.  , 

January  1967. 

28.  Reid,  C.  F.  ,  Jr.  The  Effect  of  Several  Forward-Mounted 
Control  Jet  Nozzles  on  a  Typical  Missile  Configuration  at 
Transonic  Speeds.  Cornell  Aero.  Lab.  ,  Inc.  Report 

No.  AA-2267-W-3,  July  1967. 

29.  Love,  E.  S.  ,  Grigsby,  G.  E.  ,  Lee,  L.  P.  ,  and  Woodling,  M.  J 
Experimental  and  Theoretical  Studies  of  Axisymmetric  Free 
Jets.  NASA  TR-R-6,  1959. 

30.  Crist,  S.  ,  Sherman,  P.  M.  ,  and  Glass,  D.  R.  ,  Study  of  the 
Highly  Underexpanded  Sonic  Jet.  AIAA  Journal,  roi.  4, 

No.  1,  January  1966. 

31.  Spring,  D.  J.  An  Experimental  Investigation  of  the  Interference 
Effects  Due  to  a  Lateral  Jet  Issuing  from  a  Body  of  Revolution 
Over  the  Mach  Number  Range  0.8  to  4.5.  U.  S.  Army  Missile 
Command,  Redstone  Arsenal,  Report  RD-TR-68- 10 

August  1968. 


185 


i 


32.  Vogler,  R.  D.  Surface  Pressure  Distributions  Induced  on  a 
Flat  Plate  by  a  Cold  Air  Jet  Issuing  Perpendicularly  Irum  the 
Plate  and  Normal  to  a  Low-Speed  Free-Strearn  Flow, 

NASA  TN-D-1629.  March  1963. 

33.  Bradbury,  L.  J.  S.  and  Wood,  M.  N.  The  Static  Pressure 

Distribution  Around  a  Circular  Jet  Exhausting  Normally  from  a  t 

Plane  Wall  into  an  Airstream.  Gt.  Brit.  Aero.  Res.  Council, 

C.  P.  No.  822,  1965. 

34.  Martin,  W.  and  Gelb,  G.  An  Experimental  Investigation  of 
the  Flow  Field  About  a  Subsonic  Jet  Exhausting  into  a  Low 
Velocity  Airstream.  Northrop  Corp. ,  Norair  Div.  NOR-65-229, 

August  1965. 

35.  Peake,  D.  J.  The  Pressures  on  a  Surface  Surrounding  a  Jet 
Issuing  Normal  to  a  Mainstream.  Natl.  Res.  Council  of 
Canada.  Aero.  Rept.  LR-410.  Ottawa,  November  1964. 

36.  Wooler,  P.  T.  ,  Burghart,  G.  H.  ,  and  Gallagher,  J.  T. 

Pressure  Distribution  on  a  Rectangular  Wing  With  a  Jet 
Exhausting  Normally  into  an  Airstream.  Jour,  of  Aircraft 
Vol.  4.  ,  No.  6.  November  -  December  1967. 

37.  .  Milne -Thompson,  L.  M.  Theoretical  Hydrodynamics, 

3rd  Edition.  The  MacMillan  Co.  ,  New  York,  1955.  ’ 

38.  Timmer,  H.  G.  and  Stokes,  T.  R.  Rapid  Prediction  of  Hyper¬ 
sonic  Aerodynamic  Characteristics  -  Computer  Program  F218. 

McDonnell  Douglas  Astronautics  Company--Western  Division, 

DAC-62432,  December  1968. 

39.  Simon,  W.  E.  and  Walter,  L.  A.  Approximations  for  Super¬ 
sonic  Flow  Over  Cones.  AIAA  Journal,  Vol.  No.  7,  July  1963. 

40.  Martellucci,  A.  and  Libby,  P.  A.  Supersonic  Flow  About 
General  Three-Dimensional  Blunt  Bodies  (Volume  II).  Aero¬ 
nautical  Systems  Division,  Air  Force  Systems  Command 
Report  No.  ASD-TR-61-727  (Vol.  II).  October  1962. 

41.  Needham,  D.  A.  and  Stollery,  J.  L.  Boundary  Layer  Separation 
in  Hypersonic  Flow.  AIAA  Paper  No.  66-455,  June  1966. 

42.  Hakkinen,  R.  J.  ,  Greber,  I.  ,  Trilling,  L.  and  Abarbanel,  S.  S. 

The  Interaction  of  an  Oblique  Shock  Wave  With  a  Laminar 

Boundary  Layer.  NASA  Memo.  2-18-59W,  March  1959.  1 


186 


I 


« 


* 


f 


X 


43.  Street,  D.  R.  Effects  of  Injection  Nozzle  Configuration  on 
Secondary  Injection  Into  Supersonic  Flow.  California  Institute 
of  Technology  Mechanical  Engineering  Degree  Thesis,  1966. 

44.  Burt,  J.  R. ,  Jr.  ,  and  Dahlke,  C.  W.  Effects  of  Spin  and  Vent 
Tube  Jet  Flow  on  Lance  Missile  Aerodynamic  Coefficients  from 
Analysis  of  Wind  Tunnel  Test,  LTV  Test  246  (U).  U.  S.  Army 
Missile  Command  Report  No.  RF  TM-67-4,  Redstone  Arsenal, 
Ala.  June  1967,  (Confidential). 

45.  Dahlke,  C.  W.  An  Experimental  Investigation  of  Downstream 
Flow-Field  Properties  Behind  a  Forward  Located  Sonic  Jet 
Injected  into  Transonic  Freestream  from  a  Body  of  Revolution. 

U.  S.  Army  Missile  Command  Report  No.  RD-TM-68-2, 

Redstone  Arsenal,  Ala.,  January  1968. 

46.  Dahlke,  C.  W.  An  Experimental  Investigation  of  Downstream 
Flow-Field  Properties  Behind  a  Sonic  Jet  Injected  into  Trans¬ 
sonic  Free  Stream  from  a  Body  of  Revolution  (Series  II). 

U.  S.  Army  Missile  Command  Report  No.  RD-TM-69-2, 

Redstone  Arsenal,  Ala.  ,  February  1969. 

47.  Bryson,  A.  E.  Symmetric  Vortex  Separation  on  Circular 
Cylinders  and  Cones,  journal  of  Applied  Mechanics,  Trans, 
of  ASME,  December  1959. 

48.  Hinze,  O.  ,  Turbulence,  McGraw-Hill  Book  Co.  ,  Inc.  , 

New  York,  1959. 

49.  Nielsen,  J.  N.  Missile  Aerodynamics,  McGraw-Hill  Book  Co.  ,. 
Inc.  ,  New  York,  I960. 

50.  Pitts,  W.  D.  ,  Nielsen,  J.  N.  and  Kaatari,  G.  E.  Lift  and 
Center  of  Pressure  of  Wing- Body- Tail  Combinations  at  Subsonic, 
Transonic,  and  Super  sonic  Speeds.  NASA  Report  1307,  July  1953. 

51.  Landahl,  M.  and  Ashley,  H.  Aerodynamics  of  Wings  and  Bodies. 
Addison  Wesley,  1964. 

A-l.  Adamson,  T.  C.  and  Nicholls,  J.  A.  On  the  Structure  of  Jets 
from  Highly  Underexpanded  Nozzles  into  Still  Air.  Journal  of 
the  Aero.  Sciences,  Vol.  26,  1959,  pp  16-24. 

A-2.  Ashkenas,  H.  and  Sherman,  F.  S.  The  Structure  and  Utiliza¬ 
tion  of  the  Supersonic  Free  Jet  in  Low  Density  Wind  Tunnels. 
Rarefied  Gas  Dynamics,  J.  H.  DeLeeuw  (ed.  '  Academic 
PresB,  Vol.  2,  1965. 


» 


18/ 


A -3.  Owen,  P.  L.  and  Thornhill,  C,  K,  The  Flow  in  an 

Axially-Symmetric  Supersonic  Jet  From  a  Nearly-Sonic 
Orifice  Into  a  Vacuum.  Aero.  Res.  Council  R&M  26l6, 

Gt.  Brit.,  1952. 

A-4.  Crist,  S. ,  Sherman,  P.  T.  and  Glass,  D.  R.  Study  of  the 
Highly  Underexpanded  Sonic  Jet.  AIAA  Journal,  Vol.  4, 

No.  1,  January  1966. 

A-5.  Love,  E.  S.  ,  et  al.  Experimental  and  Theoretical  Studies  of 
Axisymmetric  Free  Jets.  NASA  TR-R-6,  1959. 

B-l.  Morse,  P.  M.  and  Feahbach,  H.  Methods  of  Theoretical 
Physics,  Volume  I.  McGraw-Hill  Book  Company,  1953. 

B-2.  Voglev,  R.  D.  Surface  Pressure  Distributions  Induced  on  A 
Flat'j  Plate  by  a  Cold  Air  Jet  Issuing  Perpendicular  From  the 
Place  and  Normal  to  a  Low  Speed  Free-Stream  Flow. 

NASA  TN-D-1629,  March  1963. 


188 


UNCLASSIFIED 

Security  Classification 


nnrnufekJT  rnutofti  r» » *r  *  d  •  jj 

^•entity  cl*»  ailfrmfion  of  tl*la  *»nrf)'  **  •hatract  and  fn^axin^  annotation  mu*t  be  entered  when  the  ovrrot^e^of^^tcio^aiiled) 


1-  originating  activity  (Co rporata  author) 

McDonnell  Douglas  Astronautics  Company 
Western  Division 

3000  Ocean  Park  Bl/d,  ,  Santa  Monica,  California 


2fl.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


3b.  ONOUP 


N/A 


3  REPORT  TITLE 

Jet  Interaction  Control  Effectiveness  for  Subsonic  and  Supersonic  Flight 


4.  d£*cripiivc  not£«  (ivpi  oi  MpMi  wj  /nr/niiK  iiaKu  Final  Report-  -  Lateral  Jet  Interaction  Study- - 
Phase  IT,  1  June  1068  -  30  September  19&9 _ —  _ 

n  au  THORisi  (Ural  name,  middie  initial,  loot  name) 

Louis  A.  Cassel,  Norbert  A.  Durando,  Clark  W.  Bullard,  James  M.  Kelso 


C-  REPORT  DATE 

September  1969 

7a.  TOTAL  NO,  Or  F'AGEfl 

l  82 

76.  NO.  OF  ntFt 

58 

•A.  CONTRACT  OR  GRANT  NO. 

DAAH01-68-C- 1919 

»«.  ORIGINATOR'*  REPORT  NUM  BE  R(*l  .  . 

MDAC-G1 164 

t.  PROJEC  T  NO. 

DA  Project  Number  1M2623XXA206 

0. 

AMC  Management  Structure  Code 
*  No.  522C.  11.148 

*6.  other  REPORT  NdU  fAny  other  numbeta  that  may  be. maligned 

npa,l>  tl.  S;  Army  Missile  Command 
Report  RD-TR-69-21 

10.  OIITRISUTION  ITATIMKNT 


This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  C.  G, , 
U.  S.  Army  Missile  Command,  ATTN:  AMSMI-RD 

'll.  SUPPLEMENTARY  NOTES 


None 


19.  »PON»OAINO  MILITARY  ACTIVITY 

Advanced  Systems  Laboratory 
Research  and  Engineering  Directorate 
U.  S.  Army  Missile  Command 

RnrUtnna  Araanal Alahnma 


TS7  ABSTRACT 


Interference  effects  between  a  highly  under  expanded,  sonic  or  supersonic  jet  in  a 
subsonic  or  supersonic  crossflow,  and  the  surface  from  which  the  jet  exhausts  are 
examined.  For  subsonic  freestream  Mach  numbers,  existing  data  ^  examined  and 
correlated.  Various  semi- empirical  models  to  represent  the  interference  pressure 
distribution  on  flat  plates  are  then  developed.  For  supersonic  freestream  Mach 
numbers,  a  computer  program  for  calculating  jet  interference  effects  on  axiaym- 
metric  bodies  at  angle  of  attack  is  described.  Interference  effects  between  the  jet 
plume  and  control  fins  on  a  cruciform  missile  are  analyzed.  A  semi- empirical 
model  of  the  jet  in  a  crossflow,  valid  at  large  distances  from  the  nozzle  is  developed. 
The  results  of  this  model  are  then  used  to  compute  interference  forces  and  moments 
on  fins  located  aft  of  the  nozzle,  both  for  subsonic  and  for  supersonic  freestream 
Mach  numbers, 


DD  :r..1473 


191  UNCLASSIFIED 

Security  ClusMIication 


SSIFIE 


Security  Clafiiiflcetlon 


KEY  WORDS 

Aerodynamics 
Fluid  Mechanics 
Gas  Dynamics 
Fluid  Injection 
Gas  Injection 
Secondary  Injection 
Gas -Gas  Interaction 
Surface  Interaction 
Reaction  Control 

Reaction  Jet  Attitude  Control  Techniques 

Control  Device 

Jet  Augmented  Wing  Flap 

Jet  Control 

Jet  Flap 

Jet  Mixing 

Reaction  Jet 


nOH  1  WT  I  ROLE  I  WT 


UNCLASSIFIED 

Security  CUi'lflcition 


